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Vorwort 


Die mit Band 1 begonnene Reihe Anwendung programmierbarer Taschenrechner wird mit dem 
vorliegenden Band 2 zu den Fachgebieten Allgemeine Elektrotechnik, Nachrichtentechnik und 
Impulstechnik fortgesetzt. 

Dieser Band wendet sich insbesondere an Ingenieure der Elektrotechnik aus Entwicklungs-, Labor- 
und Betriebsbereichen, die an einer rationellen und systematischen Bearbeitung sich wiederholender 
Aufgabenstellungen aus den genannten Fachgebieten interessiert sind. Studenten an Technischen 
Hochschulen und Fachhochschulen erhalten einen Einblick in die Möglichkeiten und Grenzen der 
Anwendung programmierbarer Taschenrechner. 

Dem Leser wird ein Nachschlagefundus mit Beispielen für programmierte Aufgabenlösungen prak- 
tischer Anwendungsfälle in die Hand gegeben. Die hierbei gegebenen Programmierhinweise haben 
den Zweck, die Ausarbeitung spezieller, auf die eigene Problemstellung zugeschnittene Programme 
zu erleichtern. 

Zu den Aufgabenstellungen werden die mathematischen Grundlagen nur soweit dargestellt, wie es 
zur Programmierung des Lösungsalgorithmus erforderlich ist. Zum tieferen Eindringen in die jewei- 
lige Thematik wird auf die vorhandene Fachliteratur verwiesen. 

Die angegebenen Programme sind auf die UPN-Technik HP-67 und HP-97 mit externem Magnet- 
karten-Speicher zugeschnitten. Mit Hilfe der angegebenen Lösungsgleichungen und den Hinweisen 
zu der jeweiligen Programmstruktur sind die Programme jedoch auch auf andere Rechnersysteme 
übertragbar. Der Verfasser geht davon aus, daß die den Rechnern aller Systeme zugehörigen Be- 
dienungshandbücher und die hierzu vorliegende Sekundärliteratur alle Fragen zur eigentlichen Be- 
dienung und praktischen Einübung umfassend beantworten. 

Der vorliegende Band soll als Lehr- und Arbeitsbuch eine Brücke zwischen den objektiven Möglich- 
keiten der Technik (Hardware) und den subjektiven Fähigkeiten zu einer rechnergerechten Problem- 
analyse und rationellen Problemlösung (Software) bauen. 

Dank gebührt auch für diesen Band Herrn Dipl.-Ing. J. Hildebrandt für viele fruchtbringende Fach- 
diskussionen und Anregungen. Weiter habe ich wiederum zu danken Frau E. Lindenlauf für die 
Übernahme der Schreibarbeiten und Frau G. Stüsser für die Ausführung der Zeichenarbeiten. Dan- 
ken möchte ich auch Herrn H.-J. Niclas, Lektor im Verlag Vieweg, für die fortdauernd gute Betreu- 
ung der Arbeit und dem Verlag für die gute Aufnahme. Nicht zuletzt gilt besonderer Dank meiner 
Frau und meinen Kindern, durch deren Verständnis mir auch die Bearbeitung dieses Bandes wieder- 


um möglich wurde. 


Aachen-—Brand, Juli 1979 
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1 Einführung 


1.1 Allgemeines 


Der Band 1 dieser Reihe enthält ein Angebot von Programmen zur Lösung allgemein-mathematischer 
Aufgabenstellungen mit den programmierbaren Rechnern HP 67/97/41C (UPN-Technik). In diesem 
Band wird das Angebot auf den Fachbereich Elektrotechnik für die Fachgebiete Allgemeine Elektro- 
technik, Nachrichtentechnik und Impulstechnik erweitert. Zu den aufgegriffenen Problemstellungen 
werden alle zur Beschreibung einer programmierten Lösung erforderlichen Beziehungen angegeben. 
Mit Hilfe der zu jedem Thema angegebenen Fachliteratur kann der Leser auch ein vertieftes Verständ- 
nis der Lösungsalgorithmen erlangen. Die Arbeitsweise der Programme sollte zunächst mit Hilfe der 
angegebenen Test- oder Anwendungsbeispiele erprobt werden. In Verbindung mit der Programmbe- 
schreibung und der Anweisungsliste sind individuelle Programmänderungen oder Programmergänzun- 
gen leicht möglich. Die zu einigen Themen aus den Lösungswerten erstellten Diagramme sollen auf- 
zeigen, wie durch vielfältige Variantenrechnungen ein umfassendes Verständnis der physikalischen 
Gegebenheiten erzielt werden kann. 

Da bei der Programmierung stets von dem exakten Lösungsansatz ausgegangen werden kann, besteht 
bezüglich der Größenordnung der Eingabeparameter keine Anwendungsbeschränkung. Dies wird be- 
sonders im Bereich der Leitungsgleichungen deutlich, die sowohl im Niederfrequenzbereich wie 

auch im Hochfrequenzbereich unverändert anwendbar sind. Die Herleitung bereichsweise zuge- 
schnittener Näherungsgleichungen ist damit nicht mehr erforderlich. Um das Anwendungsspektrum 
der einzelnen Programme möglichst weit zu fassen, werden die Ansätze allgemein formuliert. Dies 
hat natürlich zur Folge, daß im konkreten Anwendungsfall einzelne Parameter theoretisch gleich O 
oder © gesetzt werden müßten. Durch null-setzen einzelner Parameter können aus dem allgemeinen 
Ansatz heraus rechentechnisch Divisionen durch null entstehen, die zu einem Fehlerabbruch des 
Programms führen. Dieser Mangel kann durch die Wahl extrem kleiner Eingabeparameter (z.B. 10"'° 
statt O) beseitigt werden. An Stelle des Wertes © kann z.B. 10'° gesetzt werden. Hierbei ist zu be- 
achten, daß bei Register-Arithmetik-Operationen das Ergebnis kleiner als 9,99 - 10° sein muß, da 
sonst die Programmbearbeitung mit ERROR-Anzeige unterbrochen wird. 

Außerdem ist bei der Rechnung mit um Größenordnungen > 10 verschiedenen Werten zu beachten, 
daß durch Addition oder Subtraktion mit anschließenden Multiplikationen unbrauchbare Rechen- 
ergebnisse entstehen können. In solchen Fällen ist ggf. eine Grenzwertbetrachtung zu empfehlen. 
Die Gutartigkeit des programmierten Lösungsalgorithmus kann durch schrittweises Herantasten an 
den jeweiligen Extremwert untersucht werden. Dabei muß der Graph der Lösungswerte einen 
asymtotischen Verlauf zu dem wahren Lösungswert zeigen. 

Nach einem Fehlerabbruch ist es zweckmäßig zu prüfen, ob vor einem neuen Start des Programms 
der richtige Speicherbereich (Primär- oder Sekundärspeicher) gemäß den Eingabebedingungen aktiv 
ist. Diese Prüfung kann unterbleiben, wenn im Programm eine automatische Ausrichtung der 
Speicherbereiche programmtechnisch eingebaut ist, oder durch konsequent angewandte indirekte 
Adressierung der Speicher auf die Anwendung der Speicherbereichswechselanweisung verzichtet 
wurde. Bei der Dokumentation der Programme ist es zweckmäßig, in den Anweisungslisten die 
Bereiche, in denen der in Bezug zum Programmstart sekundäre Speicherbereich aktiv ist, mit einer 
senkrechten Strichmarkierung zu kennzeichnen. 


1.2 Modularer Programmaufbau 


Bei der Programmierung umfangreicher Aufgabenstellungen ist es zweckmäßig, das Programm in 
Form einzelner, eigenständig ausgetesteter Module aufzubauen. Die Fehlersuche kann dann auf die 
Richtigkeit der Übergabeparameter an den Einbindungsstellen der Programmodule beschränkt 
werden. 

Insbesondere erweist sich hierbei die Verwendung ausgetesteter Unterprogramme als sehr vorteilhaft. 


Da die Rechner nur mit absoluten Adressen arbeiten, ist der Einbau übernommener Programmodule 
immer dann unproblematisch, wenn in diesen Modulen alle Operationen ausschließlich im Stack- 
Register ablaufen und keine Sprungadressen, Flag-Abfragen oder Speicheraufrufe enthalten sind. 
Anderenfalls muß auf das Vorhandensein doppelter Adressen besonders geachtet werden. 


1.3 Standard-Unterprogramme 


Standard-Unterprogramme sind Programmodule mit einer definierten Startadresse und einer RTN- 
Endanweisung. Sie dienen zur Realisierung allgemeiner Funktionen der Übergabeparameter. Hierzu 
gehören z.B. komplexe Grundoperationen und Hyperbelfunktionen mit reellem und komplexem 
Argument oder auch Hilfsprogramme zur Erzeugung besonderer Markierungen. Für den Benutzer 
sind vollständige Angaben über die Form der Übergabeparameter und der Ergebnisparameter sowie 
der intern benötigten Adressen, Speicher, Flags zu machen. Bei den Adressenangaben ist außerdem 
von Bedeutung, ob diese innerhalb des Unterprogramms in Vorwärtsrichtung oder in Rückwärts- 
richtung zu erreichen sind. Die letzteren dürfen weder im Hauptprogramm noch in weiteren Unter- 
programmen auf keinen Fall mehrfach verwendet werden. 


1.3.1 Komplexe Grundoperationen 


Die Operanden zur Ausführung der komplexen Grundoperation (Addition, Subtraktion, Multiplika- 
tion, Division, Inversion) können in Polarkoordination oder in kartesischen Koordinaten vorgegeben 
werden. Bei der Subtraktion und der Division ist es wichtig, daß der Minuend bzw. der Divisor mit 
Betrag und Phase oder Real- und Imaginärteil jeweils im X- und Y-Register übergeben werden. Die 
komplexe Ergebnisgröße wird mit Imaginärteil bzw. Winkel im Y-Register und mit Realteil bzw. 
Betrag im X-Register an das rufende Programm übergeben. Wegen der besonderen Bedeutung dieser 
Unterprogramme für die Elektrotechnik wird die schon in Band 1 dieser Reihe enthaltene Anwei- 
sungsliste zur Realisierung der komplexen Grundoperationen für die Varianten mit ausschließlich 
Stack-Register-Operationen in Tabelle 1.3.1 nochmals angegeben. 
Die einzelnen Grundoperationen ; n'Palar- in kartesischen 
sind über folgende Adressen Operation Köardinsten erlinaien 
aufzurufen: 

Ei 


a 
+ 


Iw 


Tabelle 1.3.1 Anweisungsliste komplexe Grundoperationen 


Polarkoordinaten Kartesische Koordinaten 
062 SR 852 Fr37 
083 R+ 853 Rr 
884 RE a54 + 
805 >R 855 R+ 
Ba6 xeY e56 + 
ae7 [23 057 kt 
656 + 85E RTN 
65 Ri 
810 + 868 er 
811 Rt 861 Rs 
al2 »P 662 = 
613 RIN 863 R+ 

064 KEY 
815 R 665 = 
el Ri 866 kt 
617 Ru 867 RTN 
8 @R 
619 Kay 863 EY2 
828 Kr 878 R+ 
e2i Key 871 Ri 
622 = 872 +F 
823 Rs 873 er 
824 = 874 R+ 
825 Rt 75 x 
86 >P ar6 Rt 
827 RTN 877 + 

878 kt 
829 xEY 879 -R 
838 RE 888 RTN 
831 x 
832 R+ 882 F 
633 + 883 R+ 
834 Rt 084 Rt 
435 RTN 885 >P 

&86 er 
837 x 887 Rı 
835 Rs 688 ker 
833 = 869 = 
934 RE 638 R+ 
841 KEY es - 
d42 - 892 Rt 
643 kt 893 »R 
644 RTN 494 RTN 
846 FR 896 »F 
847 R 897  6SBE 
846 CHS 898 >R 
849 Rt 093 RTN 


658 RTN 


1.3.2 Hyperbelfunktionen 


a) Hyperbelfunktionen mit reellem Argument 


Die Hyperbel- oder hyperbolischen Funktionen sind folgendermaßen definiert [1,2]: 


sinhx = Lei -e”) 11.3.1) 
coshx = I eite”) (1.3.2) 
tanhx = snz (1.3.3) 
cothx = en (1.3.4) 


Für die Umkehrfunktionen (Areafunktionen) gilt: 


arsinnx = In(x + Yx? +1) (1.3.5) 

arcoshx = In{x+ Yx? - 1) für x2 1 (1.3.6) 

artanhx = m 2% für -I1sxsi1 (1.3.7) 
2 1-x 

arcothx = zn für -12x und x21 (1.3.8) 


Der Übergabeparameter x und der Ergebnisparameter der Funktion werden jeweils im X-Register 
übergeben. 


b) Hyperbelfunktionen mit komplexem Argument 


sinh Z = sinh (x + jy) = 5 [te -e’*) cosy +jle*+te’”) siny] (1.3.9) 


cosh Z = cosh (x + jy) - 4 [ie* +e*)cosy+jle“-e”*)siny] (1.3.10) 


sinh x 


Der komplexe Übergabeparameter Z und der komplexe 
Ergebnisparameter der Funktion werden in kartesische 
Koordinatenform mit dem Imaginärteil im Y-Register 
und dem Realteil im X-Register übergeben (der Übergabe- 
parameter y muß beim Grad Mode (DEG) in Winkelgrade 
übergeben werden). 


cosh x 


tanh x 
In der Tabelle 1.3.2 ist ein Programm zur Berechnung der 


Hyperbelfunktionen über folgende Adressen angegeben: coth x 
arsinh x 
arcosh x 
artanh x 


arcoth x 


sinh Z 


DIUEEBIGIBECHGIGHE 


cosh Z 


Tabelle 1.3.2 Anweisungsliste „‚Hyperbelfunktionen‘ 


sinhx ai »LELA 847 ‘X 893 + 
802 ENT 848 - 894 LSTX 
803 e* 849 LSTK 895 Rt 
884 aY es xer 896 + 
8655 CHS a51 R+ 897 Rt 
886 er 852 + 898 Sin 
887 - 853 LN 899 x 
888 5 854 Rt 188 2 
803 * 855 LN 181 # 
eid8__RIN : T 182 Fi 

Coshx el #LELE 183 x8y 
612 ENTF 055 ENT 184 cüs 
813 er 08593 ENT? 185 x 
814 Kar 868 1 186 2 
815 60H 861 + 167 ® 
8l6 er a2 x 183 Rt 
617 + 863 cHs 189 xy 
818 2 864 1 118 _RIN 
619 ® 865 + coshix+jy)_ 111 #LBLe 
626 __RIN 866 & 112 ENT? 

tarhx 821 #LBLtl 867 LN 113 CHS 
822 ENT 865 2 114 er 
823  65BR 863 ® 115 ENT 
024 XKer 6e__RIN 116 RY 
025  65BB arcoth x 871 __x*LBLd 117? Kay 
825 & 872 ENT? 118 er 
827__RIN 873 ENT? 118 + 

cofh x 826_ xLBLD 674 N 128  LSTa 
825  6SBC 75 & 121 Rt 
0838 1X 876 xer 122 C#5 
831 RTN 877 1 123 + 

arsinhx 032 _ #LBLa 878 - 124 Rt 
033 ENT? 873 3 125 SIN 
834 K2 888 LN 126 x 
835 1 861 2 27 2 
836 + 882 5 128 n 
837 ‘x 8832 __RTN 123 R+ 
836 + sinhixsjy)_ 884 #LBLE 138 Kay 
839 LN 885 ENT? 131 cos 
848 ___RTN BBE CH 132 x 

orcoshx 841 _ XLBLE 887 ex 133 2 
842 ENTt es8 ENT? 134 * 
843 ENT? 089 Rr 135 Rt 
844 x? e38  cHs 136 xer 
845 1 821 ar 13? RIN 
846 - 892 er 


1.3.3 Periodische Schaltfunktionen 


Bei der Untersuchung periodischer nicht-sinusförmiger Vorgänge, treten häufig als Eingangsfunk- 
tionen periodische Schaltfunktionen in Rechteck-, Dreieck- oder Sägezahnform auf. Derartige Funk- 
tionen lassen sich vorteilhaft aus den Kreisfunktionen und deren Umkehrfunktionen gewinnen. 

Die Periode von 27 im Bogenmaß, entsprechend 360° im Gradmaß, ist rechnerintern vorgegeben. 
Für den Winkel-Modus dezimale Grad {Normalzustand nach dem Einschalten des Rechners) sind 
nachfolgend einige Schaltfunktionen angegeben: 


1. Rechteckfunktion 


Die Rechteckfunktion, Bild 1.3.1, ergibt sich durch Multiplikation einer Konstanten mit der Vor- 
zeichenfunktion sign der Sinusfunktion 


y = sign (sin x) (1.3.11) | ' 
Bild 1.3.1 Rechteckfunktion -1 


Die Rechteckfunktion wird durch den Programmbaustein mit der Startadresse Label O realisiert: 


LBLO x#0 
SIN + 
ABS RTN 
LAST x 


2. Dreieckfunktion 


Die Dreieckfunktion, Bild 1.3.2, ergibt sich durch Anwendung der inversen Sinusfunktion auf die 
Sinusfunktion selbst: 


‚_—1i “ 
ya; Ar ini (1.3.12) = 
90 


Bild 1.3.2 Dreieckfunktion 


Die Dreieckfunktion wird durch den Programmbaustein mit der Startadresse Label 1 realisiert: 


LBL 1 90 
SIN ja 
sin! RTN 


3. Sägezahnfunktion 


Die Sägezahnfunktion, Bild 1.3.3, ergibt sich durch Multiplikation der inversen Cosinusfunktion 
angewandt auf die Cosinusfunktion selbst mit der Vorzeichenfunktion angewandt auf die Sinus- 
funktion: 


cos ' (cosx) 


y- 180 sign (sin x) (1.3.13) 


Bild 1.3.3 Sägezahnfunktion 


Die Sägezahnfunktion wird durch den Programmbaustein mit der Startadresse Label 2 realisiert: 


LBL2 x2y 
ENTER cos 
SIN cos“! 
ABS x 
LAST x 180 
x#0 = 

+ RTN 


4. Trapezfunktion und lückende Dreieckfunktion 


Die Trapezfunktion und die lückende Dreieckfunktion, Bild 1.3.4, lassen sich durch Kombination 
der Hauptwerte der Cosinusfunktion mit dem Betrag 
der Hauptwerte der Sinusfunktion realisieren: 


u cos"! (cosx) + Isin”' (sin x) | 1.3.19) 
’ 180 ar Ba Ca "zu Bi 
+ Zeichen: Trapezfunktion o) 
- Zeichen: lückende Dreieckfunktion 
H 1 
Bild 1.3.4 i 
a) Trapezfunktion 0° 180° 360° n 
b) lückende Dreieckfunktion b) 


Die Trapezfunktion und die lückende Dreieckfunktion werden durch die Programmbausteine mit 
den Startadressen Label 3 und Label 4 realisiert: 


a) Trapezfunktion b) lückende Dreieckfunktion 
LBL3 sin! LBLA4 sın! 
ENTER ABS ENTER ABS 
cos + cos = 
cos"! 180 cos"! 180 
x2y 7 x2y = 
SIN RTN SIN RTN 


5. Lückende Sägezahnfunktion 


Eine lückende Sägezahnfunktion, Bild 1.3.5, mit 180° Pulsbreite und 180° Pulspause läßt sich durch 
Kombination der Sägezahnfunktion mit einer Dreieckfunktion realisieren: 


-1 h s 
cos  (cosx) + [1 + sign (sin x)] 
= 1.3.15 
r 360 yı 
0 
0° 180° 360° BEBEEN. 
Bild 1.3.5 Lückende Sägezahnfunktion 


Die lückende Sägezahnfunktion wird durch den Programmbaustein mit der Startadresse Label 5 
realisiert: 


LBL5 x#0? 
ENTER + 
cos 1 
cos"! + 
ENTER Rt 
ENTER x 

At 360 
SIN = 
ABS RTN 
LAST x 


Ersetzt man in Gl. 1.3.15 das Plus-Zeichen durch ein Minus-Zeichen, so entsteht eine lückende 
Sägezahnfunktion mit umgedrehtem Verlauf und negativen Werten (senkrechter Anstieg und ge- 
neigter Abfall) im Intervall 180 < x < 360°. Die negativen Amplitudenwerte können im Programm 
durch eine Vorzeichenwechselanweisung vor dem Rücksprung umgekehrt werden. 


6. Periodisch alternierende Nadelfunktion 


Eine periodisch alternierende Nadelfunktion, Bild 1.3.6, läßt sich durch zusätzliche Vorzeichen- 
steuerung aus dem Ergebnis einer „‚Größer null’‘-Abfrage der Cosinusfunktion erzeugen. In dem 

Programmodul wird hierzu als Entscheidungsmerkmal das Flag 2 benutzt. 

Die alternierende Nadelfunktion wird durch den Programmbaustein mit der Startadresse Label 6 
realisiert: 


LBL6 x#y Y 

COS 0 e 180° 

x<o F2? p° ER 

SF2 CHS R 

ABS RTN 

1 Bild 1.3.6 Alternierende Nadelfunktion 


Falls der x-Maßstab für die dargestellten Funktionen das Periodenintervall 1 haben soll, ist nach 

dem jeweiligen Start-Label eine Multiplikation mit 180 bzw. 360 einzufügen. Für beliebige Perioden- 
intervalle | ist das Argument x vor Aufruf des entsprechenden Programmoduls mit = bzw. 280 

zu multiplizieren [KI<xS (k + 1) 1]. 

In der Tabelle 1.3.3 sind die Programmodule der einzelnen Schaltfunktionen zusammengefaßt. 


1.4 Zahlenbereichsüberschreitung 


Unabhängig von dem gewählten Anzeigeformat werden rechenintern 10 signifikante Ziffern einer 
Zahl zuzüglich einem zweistelligen Exponenten verarbeitet. Die darüber hinausgehenden Ziffern 
werden abgeschnitten. Die Addition oder Subtraktion einer Zahl kleiner 0,5 zu einer zehnziffrigen 
Zahl bleibt z.B. auch bei beliebig häufiger Wiederholung wirkungslos. Eine ziffernmäßig wirksame 
Verarbeitung von zwei Zahlen ist nur dann für alle Operationen möglich, wenn sich beide Zahlen 
um weniger als 10'° voneinander unterscheiden. Durch diese Beschränkung der internen Ziffern- 


verarbeitung auf 10 signifikante Ziffern bleibt der extern zulässige Zahlenbereich von + 9,99 109% 
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Tabelle 1.3.3 Anweisungsliste „Schaltfunktionen‘ 


Rechteck-F. Bu xLBL& lück Dreieck 845_ xL6 


22  SIN B48  ENTT 
ee ABS 87 008 
BB LSTK B4E COS“ 
a05  Xr6r 849 Kr 
086 ei esa Sin 
a67_RTN 51 SIn 
Dreieck-F__BBE_#L6L1 eSE ABS 
a8 SIN 653 - 
dia  SIN 054 ı 
61 $ 855 & 
el: e esc ö 
613 * e57 # 
Bl __RTN 65$__RTN 
else ENT? B60  ENTT 
ei?  5SIn 861 cos 
616 ABS 6: Cüs- 
219 LSTK @62  ENTt 
826 x#0° 864  ENTI 
821 = 865 Rt 
BEE Kar @6E  SIN 
823  cüs 867 ABS 
824 COS" 868 LSTÄ 
82 x 869 XrB® 
826 1. 878 5; 
27 8 erl 1 
628 & 72 + 
623 ä 873 Rt 
&30___ RTN a74 x 
Trapez-F. 851 »LBL3 075 Ki 
032  ENTT are 5 
033 C08 arr & 
034 Cüs" e78 5 
835 ey @3___KIN 
e36  SIN 
837 SIn“ a1 cos 
63E ABS BB: Ne 
839 + ass  SFE 
848 1 B4 RBS 
84 8 885 ! 
842 ö BEE Arr” 
843 5 687 ) 
644 RIN E88 Fa 
as chs 

a3a  RTN 


unberührt. Dieser Sachverhalt kann für Rechenoperationen mit Additionen oder Subtraktionen, bei 
denen Operanden vorkommen, die sich um mehr als 10'° voneinander unterscheiden, zu vollkom- 
men falschen Ergebnissen führen. An zwei Beispielen soll das ‚„‚Umkippen‘ der Rechnung durch 
Operationen an der Zahlenbereichsgrenze aufgezeigt werden: 


a) Die Gleichung 


El h)-(-h) von 


hat für alle x die Lösung y = 1. 


Dagegen liefert die schrittweise Lösung gemäß dem Programm mit der Startadresse Label A nach 
Tabelle 1.4.3 in Abhängigkeit von x (STOA) und n (STO B) die in Tabelle 1.4.1 angegebenen 
Lösungswerte. Nur die Werte oberhalb der treppenförmigen Begrenzungslinie sind richtige Lösungen 
der Gleichung. Es ist zu ersehen, daß z.B. die Kombination der Werte x = 100 und n=5 zu einem 
falschen Ergebnis führt, ohne das diese Fehlrechnung durch eine Fehlermeldung unterbrochen wird. 


Tabelle 1.4.1 Lösungsmatrix zu Gl. (1.4.1) 


b) Die Gleichung: 


N 2 ' 
= 1- ]X1- =) 1.4.2 
ME | | x” 
hat für große x die Näherungslösung: 
1\] 
x" 11-1 -5) 1.4.3 
ae 1 | we] 


mit dem Grenzwert y=1 für x"— =. 
Die schrittweite Lösung der Gl. (1.4.2) gemäß dem Programm mit der Startadresse Label B und 
Tabelle 1.4.3 liefert die in Tabelle 1.4.2 angegebenen Lösungswerte. 


Tabelle 1.4.2 Lösungsmatrix zu Gl. (1.4.2) 


| 10? E 10° | 10° 


1,06 1,0 
101 1 
1 
1 
1 


1 


1 
1 


1 
1 
1 
0 
1 0 


Die Lösung zeigt asymtotische Annäherung an den Grenzwert y = 1 für x” + ® bis zu den Operan- 
denkombinationen an der Zahlenbereichsgrenze von 10'°. Hier kippt die Rechnung um, d.h. bei 


der Berechnung des Radikanten 1 - 2 fällt der Subtrahend aus dem Rechenregister nach rechts 


heraus und wird ignoriert. Die hierdurch entstehende geringfügige Ungenauigkeit liefert als Ergebnis 
der Wurzel den Wert 1 und somit das falsche Ergebnis y=0O. 

Anhand dieser beiden Beispiele sollte deutlich werden, daß bei Nichtbeachtung der Zahlenentwick- 
lung innerhalb der Rechnung beliebig falsche Ergebnisse entstehen können. Da bei den üblichen 
technischen Berechnungen der ziffernmäßig verfügbare Zahlenbereich von 10"° nicht verlassen 
wird, treten rechentechnische Probleme dieser Art nur in speziellen Fällen auf, die dann jedoch 
einer sorgfältigen Grenzwertbetrachtung bedürfen. 
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Tabelle 1.4.3 Anweisungsliste „Zahlenbereichsüberschreitung” 


Start 801 _ «LBLA 


682 
803 
884 
885 
886 
087 
888 
889 
818 
011 
812 
813 
814 
815 
el6 
817 
818 
819 
828 
821 
822 
023 
024 
825 
826 


Start 027 

1 e28 
RCLA 823 
PRTX=D |x 838 
RCLB 631 
Yx 832 
$T08 833 
17x 834 
+ 835 
LSTX 836 
CHS 837 
1 838 

+ 833 

> 848 
RCLE 841 
2 842 

* 043 

x 044 
PRTS=D |y 845 
SPC 846 
RCLA 847 
1 646 

6 842 

x 858 
STOA 851 
£TOA 852 


&LBLB 
2 
RCLR 
RCLB 
Y* 
CHS 
I 

+ 
‘x 
CHs 
l 

+ 
RCLA 
PRTX=D |x 
RCLB 
yr 

x 
PRTX => |y 
SPC 
RÜLA 
1 
[ 


x 
STOA 
&TOB 


2 Allgemeine Elektrotechnik 


2.1 Netzwerkberechnung 


Die Berechnung linearer elektrischer Netzwerke mit den Methoden der Maschen- oder Knotenana- 
Iyse führt durch Anwendung der Kirchhoffschen Regeln zu linearen Gleichungssystemen, zu deren 
Lösung bis zur vierten Ordnung in Band 1 ein Programm angegeben wurde. Bei der Maschenanalyse 
werden die Ströme in den Verbindungszweigen als unabhängige Veränderliche aufgefaßt. Jede 
Masche enthält einen Verbindungszweig mit dem unbekannten Strom und sonst nur Baumzweige, 
die von mindestens zwei Maschenströmen durchflossen werden. Die Ströme in den Baumzweigen 
lassen sich mit Hilfe der Knotenregel sofort angeben, wenn die unabhängigen Ströme in den Ver- 
bindungszweigen berechnet sind. 

Ein Netzwerk mit k Knoten und z Zweigen enthält b=k- 1 Baumzweige und n=z-b Verbin- 
dungszweige. Das Gleichungssystem für die Maschenanalyse umfaßt n=z-b=z-k + 1 Gleichun- 
gen für die n unbekannten Maschenströme [3, 4]. 

Bei Anwendung der Knotenanalyse wird das gegebene Netzwerk durch b=k- 1 Gleichungen für 
die Spannungen an den Baumzweigen beschrieben. Zur praktischen Anwendung der Knotenanalyse 
werden nicht die Zweigspannungen als unabhängige Veränderliche eingeführt, sondern die Potentiale 
an b = k - 1 Knotenpunkte gegenüber dem mit Null-Potential gewählten Bezugsknoten. 

Der Rechenaufwand für die Anwendung einer der beiden Methoden ist abhängig vom Aufbau des 
gegebenen Netzwerkes. Bei der Maschenanalyse ist die Anzahl der Gleichungen identisch mit der 
Anzahl der Verbindungszweige. Die Anzahl der Verbindungszweige hängt bei gegebener Knoten- 
zahl davon ab, wie stark das Netzwerk „besetzt” ist. Aus der Netzwerktheorie ist bekannt, daß ein 
Netzwerk mit k Knoten maximal 5 k {k - 1} Zweige enthält. Davon sind maximal: 


Han Dei [=] (2.1.1) 


Verbindungszweige. 

Die Zahl der Gleichungen bei der Maschenanalyse wächst also mit der Knotenzahl quadratisch an. 
Bei den praktisch vorkommenden Netzwerken wird dieser Maximalwert jedoch meist bei weitem 
nicht erreicht. Für die Mindestzahl der Verbindungszweige eines Netzwerkes als untere Grenze für 
den Umfang des Gleichungssystemns gilt: 


S+ 1 für k gerade 
Amin” (2.1.2) 
+ 
k 2 . +1 für k ungerade 


Bei der Mindestzahl ist unterstellt, daß bei gerader Knotenzahl an jedem Knoten drei Zweige und 
bei ungerader Knotenzahl an einem Knoten vier Zweige und an den übrigen Knoten drei Zweige 
angeschlossen sind. 

Die Zahl der für die Maschenanalyse benötigten Gleichungen wächst bei einem minimal besetzten 
Netzwerk linear mit der Knotenzahl an. Bei der Knotenanalyse wächst die Zahl der Gleichungen 
wegen b=k- 1 linear mit der Knotenzahl an und ist unabhängig von der Besetztheit des Systerns 
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durch Zweige. Aus diesem Sachverhalt ergeben sich insbesondere für Netzwerke der Energiever- 
sorgung Vorteile für die Knotenanalyse [5, 6]. 


2.1.1 Netzwerke mit reellen Widerständen 


Zur Berechnung einfacher Netzwerke mit reellen Widerständen läßt sich der Lösungsalgorithmus für 
Netzwerke bis zu drei Maschen und vier Knoten in dem verfügbaren Programmspeicher mit 224 An- 
weisungen unterbringen. Hierzu soll von einem allgemeinen Netzwerk mit sechs Spannungsquellen 
und sechs Widerständen gemäß Bild 2.1.1 ausgegangen werden. 


Bild 2.1.1 
Allgemeines 3-Maschen-Netzwerk 
mit 4 Knoten und 6 Zweigen 


Die Maschengleichungen für die drei Maschen M,, M, und M3 führen zu folgendem Gleichungs- 
system mit den unbekannten Zweigströmen I,,12 und 13: 


(MI): IH(Ry+ Rs+ Re) Io Re - IaRs =U,+W,-U. 
(M2): Re +12 (Ra+ Ra+ Re) — IaRa =U+U.-U, (2.1.3) 
(M3): -14 Rs -13R4 + 1,{Ra+ Ra+ Rs) = Ug+ U,- U, 


Das Gleichungssystem (Gl. 2.1.3) läßt sich in Matrixschreibweise wie folgt angaben: 

R-1I=U* (2.1.4) 
Mit der Lösungsgleichung für den Stromvektor: 

I=R'-U* (2.1.5) 


Zur Lösung muß aus der quadratischen Widerstandsmatrix R die invertierte Matrix A” gebildet 
werden. 
Für die Widerstandsmatrix R gilt: 


Rı+tRstRg Re Rs a) bı cı 
R= Re Ra+ Ry+ Re -Ry =» b & (2.1.6) 
Rs; -R4 Rat Rıt Rs 93 ba Ca 


Für den Spannungsvektor U* gilt: 


U, R U, U. 
U*=|U,+U.-U, (2.1.7) 
U;+U,-% 


Für den Ergebnisvektor 7 gilt: 


I=|1, (2.1.8) 
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Die Berechnung der invertierten Matrix R' kann mit einem geschlossenen Formelsatz durchge- 
führt werden. Voraussetzung dazu ist, daß die Determimante D = |R| der Widerstäandsmatrix 


existiert. 


a ßı Yı 
|. hy 
a3 Ba Ya 


Die Determinante der Matrix |R| ergibt sich zu: 


D =4, bac3 + bı Ca; +tcC} b332 -cCı b2 a3 -c2b34ı -C3240b} 


Für die Elemente der invertierten Matrix R”' gilt: 


2.1.2 Programmstruktur 


(2.1.9) 


(2.1.10) 


a = (b2c3 - b3ca)/D Pı = (bacı -byc3)/D Yırlbic2 - bacı)/D 
a = (a3C2 - 22 c,)/D ß> = (aı C3 - a3cı)/D Y=la2cı - a162)/D (2.1.11) 
a3 = (aa b3 - a3b>)/D Ba = (a3bı - a} b,)/D Y,= (aıb2 - 22b,)/D 


Das Programm soll neben der eigentlichen Netzwerkberechnung auch die Möglichkeit des Aufrufes 
allgemeiner Matrizenoperationen, die ohnehin als Unterprogramme benötigt werden, bieten. Das 
über Label E aufzurufende Netzwerkprogramm enthält daher folgende eigenständige Unterpro- 


gramme: 
Label A: 


Label B: 
Label C: 
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Bild 2.1.2 Flußdiagramm „Reelle Netzwerkberechnung” 


Berechnung der Determinante D einer 3 X 3 Matrix A gemäß GI. (2.1.10) 
Berechnung der inversen Matrix A" zu einer gegebenen Matrix A gemäß Gl. (2.1.11) 


Berechnung des Matrizenproduktes C' einer quadratischen Matrix A mit einem Spalten- 
vektor B durch Zeilen-Spaltenmultiplikation C=A-B. 


Berechnung der 
Widerstandsmatrix R 
Determinante |R]| 
UP: Label A 
Matrix-Inversion R"' 
UP: Label B 


Matrizen-Multiplikation 
UP: Label C: /=R'-U* 


Ausgabe: 
I7,12,13, U4,Pı 


Vor dem Aufruf dieser Unterprogramme müssen die Elemente der quadratischen Matrix A in den 
Speichern STO 1 bis STO 9 in der Reihenfolge a}, a2, 33, bı., D2, Da. Cı., Ca, C3 (spaltenweise) 

und die Elemente des Spaltenvektors B in den Speichern STOA bis STOC für das erste, zweite 
und dritte Element vorliegen. Bei der Berechnung der inversen Matrix A! wird die gegebene 
Matrix durch die berechnete inverse Matrix überschrieben. Die Ergebnisse der Produktberechnung 
werden in die Speicher STO D,E und I abgelegt. 

Das eigentliche Netzwerkprogramm wird über die Startadresse Label E aufgerufen. Zunächst werden 
die Elemente der Widerstandsmatrix R berechnet, um dann durch Aufruf der Unterprogramme 
Label A, Label B und Label C entsprechend dem Flußdiagramm in Bild 2.1.2 den Ergebnisvektor 
mit den Maschen-Strömen I, ,1,, Ia, der Spannung U, und der Leistung P, zu erhalten. 


2.1.3 Programmbeschreibung ‚„Reelle Netzwerkberechnung” 


Für das Programm „Reelle Netzwerkberechnung‘’, siehe Tabelle 2.1.1, mit der Startadresse Label E 
in Zeile 001 sind neun Eingabeparameter vorgesehen. Die Eingabespannungen gemäß dem Span- 
nungsvektor U* in Gl. (2.1.7) sind den Speichern STO A, STO B und STO C zugeordnet. Sind 
Netzwerk nach Bild 2.1.1 alle sechs Spannungen vorhanden, so sind folgende Eingabedaten erfor- 
derlich: 


U,+W,-U. in STOA 
U,+U.-U, in STOB 
Us+U,-U, in STOC 


Die Widerstände sind den Speichern STO 1 bis STO 6 wie folgt zugeordnet: 


R, in STO 1 R, in STO4 
R, in STO2 Rs in STO5 
Ra in STO3 Re in STO6 


Nach Start des Programms über Label E wird zunächst in einer Dekrementschleife von Zeile 004 
bis Zeile 010 eine Umspeicherung der Widerstands-Eingabedaten in die Sekundärregister STO 1’ 
bis STO 6° vorgenommen. Durch diese Maßnahme brauchen die Eingaben bei wiederholtern Pro- 
grammstart nicht erneuert zu werden. Ab Zeile 011 bis Zeile 039 ist der sekundäre Speicherbereich 
aktiv. Mit dem Programmaufruf GSBB in Zeile 037 wird die Berechnung der inversen Matrix ge- 
mäß GI. (2.1.11) eingeleitet. Innerhalb des Unterprogramms Label B ab Zeile 079 wird in Zeile 080 
ein weiteres Unterprogramm Label A zur Berechnung der Determinante gemäß GI. (2.1.10) aufge- 
rufen. Der Wert der Determinante wird in einer kurzen Pause in Zeile 081 zur Anzeige gebracht. 


Mit dem Aufruf des Unterprogramms Label C in Zeile 038 wird der berechnete Stromvektor mit den 
Komponenten I}, 2, !z ausgegeben. Mit den folgenden Anweisungen der Zeilen 040 bis 049 werden 
die Spannungen U, = I4 Ra und die von der Spannungsquelle U, abgegebene Leistung P} = U; ' 1} 
berechnet und ausgegeben. Die als eigenständige Hauptprogramme anwendbaren Unterprogramme 
Label A für die Determinantenberechnung, Label B für die inverse Matrixberechnung und Label C 
für die Produktberechnung enden jeweils mit einer RTN-Anweisung. Die Determinantenberechnung 
endet in Zeile 078 mit dem Ergebnis im X-Register. Die inverse Matrixberechnung endet in Zeile 165 
mit gelöschtem Anzeigeregister und der Ergebnismatrix in den Sekundärregistern STO 1’ bis STO 9’ 
in der Reihenfolge &,,, &2, &3, ßı , ß2. ß3, Yı. Ya, Ya. Die Produktberechnung endet in Zeile 206 
mit dem Ergebnisvektor im X-, Y- und Z-Register sowie zusätzlich im I-, E- und D-Register. Im 
D-Register ist das Produktergebnis erste Zeile mal Spaltenvektor, im E-Register zweite Zeile mal 
Spaltenvektor und im I-Register dritte Zeile mal Spaltenvektor abgelegt. 
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Tabelle 2.1.1 Anweisungsliste „Reelle Netzwerkberechnung” 


Netzwerk 


Det A 
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a6l _xLELE 657  RCLS 
a2 [3 058  RCLI 
ea3  STOI 059 6567 
a04 xLBLE e68 RELS 
065 Rüti ael RÜL? 
BRe F2S 862  65B? 
ea?  5T0i 863 CHS 
Bas FrS 864 REL3 
eos DS 065  RCLE 
eıe ETOE 866  65B? 
811 PES 067 RELI 
e12 | RÜLS 868  RCLI 
813 | 5T03 869  65SB7 
eı4 | RCLS e78  RCL2 
815 | ST+$ 871 __RCL? 
sie | STH UP 872 #LBL’ 
817? | RELE 873 1521 
eis | STOS 674 RLLi 
818 R+ 875 x 

528 CHS are x 

021 RöLE ar? + 

822 sT+ vB RTN 
823 | ST+5 K e73 _xLBLB 
824 CHS 688  65BA 
625 | RCL4 881 

82€ | ST+S 082 1/8 
Ber | ST+4 885  RCLi 
828 CHS 084 RCLS 
822 | STOs 085 x 

ese | STüc @BE  RELS 
831 Ry 887 _ RÜL? 
832 | STUs 0388 65B3 
833 | STOZ 889 STOD 
834 Ri 898 CLx 
835 | STUr 6821  RCL3 
ale | STE 692  RCL4 
837 | 656€ 893 x 

838 | ES5C 834  RÜLI 
835 P2S 895 RÜLE 
848 ÜHs 026 8583 
841  RCLE 897  STÜE 
842 + 895 CLA 
043  RÜLI 8932  RCLZ 
0944 N 188  RCL? 
85  FRTS=D U, 181 x 

#46  RÜLD 182  RELI 
a47: RÜLR 183  RCLE 
046 x 164  65B3 
049 PFTA=D Pı 185  STOI 
85% Rx 186 CLX 
851  «LBLA 187  RCLI 
852 6 188 RCLS 
853  STOI 183 x 

0854  RÜLE 118 RELZ 
6855  RCLE 111 RCL4 
856  65SE7 112  65SB3 


PSE=> Det 


up 


113  STOB 
114 CLX 
115 RCL3 
116  RCLE 
117 x 

118  REL2 
119 RCLI 
126  6SB3 
121 $T01 
122 CLX 
123  RCL2 
124 RCL6 
125 x 

126  RCL3 
12? RELS 
128 6583 
1293  5T03 
138 CLX 
131 RELS 
132 RELS 
133 x 

134 RELE 
135 RCLE 
136 6583 
137 $STO2 
138 CL 
139 RCLE 
148 RCL? 
141 x 

142 RCL4 
143 RCLI 
144 6583 
145 'STOE 
146 CLX 
147 RCL4 
148 RCLE 
149 x 

158  RCLS 
151 RCL? 
152 6SB3 
153  RCLI 
154  RCLE 
155 6SBe 
156  RCL2 
157 REM 
158  RCL3 
159  65Ba 
168  RELE 
161 RÜLD 
162  RCLE 
163 6SBb 
164 CLX 
166 #LBL3 
167 x 

168 - 


169 x 
178 RIN 
UP 171 «#LBla 
172 8 
173 6705 
174 xLBLb 
175 3 
176 6705 
177 #LBLe 
178 6 
179 &LBLS 
188 STOI 
181  6SBE 
162 65B6 
UP 183 sLBL6 
184 Rt 
165 18521 
186  STOi 
187 RTN 
189 SPC 
198 1 
191  STOI 
192 65B1 
193 STOD 
194 2 
195 STOI 
196  6SBi 
19? _STOE 
198 3 
199  $TOl 
206 65Bi 
281 $Tül 
2082 8 
2063 RÜCLD 
264 RCLE 
265 RELI 
2086 RTN 
UP 267 «LBLi 
288 6 
269 RCLA 
218 6584 
211  RCLB 
212  65B4 
213  RELC 
214 6584 II 
215 PRTX=DII 
216 RTN (I 
UP 2i7 »LBL4 
218 RCLG 
213 x 
228 + 
221 1821 
222 1821 
223 IS2I 
224 RTN 


2.1.4 Test- und Anwendungsbeispiele 


a) 


b) 


Zu der Schaltung (Bild 2.1.3) sind mit den angegebenen Daten die Ströme 1;, >, I, die 
Spannung U, und die von der Spannungsquelle U, abgegebene Leistung P, zu bestimmen: 


U, =40 V R7=100 2 R5 
U2 =30V Ra =200n 
U3=50 V R3=30N 
R4=0N 
Rs=500 
Re= 60 {7} 
Eingaben: 40 STOA 10 STO1 r 
30 STOB 20 STO2 
50 STOC 30 STO3 Bild 2.1.3 3-Maschen-Netzwerk 
40 STO4 
50 STO5 
60 STO6 


Start [El 


Ergebnisse: 1. Kurze Pause Determinante D = 564.000 

2. Lange Pause bzw. Ausdruck I, = 2,14 A 
3. Lange Pause bzw. Ausdruck 12 = 1,98 A 
4. Lange Pause bzw. Ausdruck 13 = 1,97 A 
5. Lange Pause bzw. Ausdruck U, = 0,35 V 
6 


. Lange Pause bzw. Ausdruck P, = 85,67 W 


Rechenzeit rd. 40 Sekunden 


Die Parallelschaltung von zwei Spannungsquellen mit den Leerlaufspannungen Ujo, U2o und 
den Innenwiderständen R;,, R;2 speise auf einen Belastungswiderstand Ra. Gesucht sind die 
beiden Einspeiseströme und der Spannungsabfall am Belastungswiderstand. 


Uo = U20o = 12V 


Rı=1ı9 
Ro=20 
Ag =10n 


Bild 2.1.4 2-Maschen-Netzwerk 


Das allgemeine Netzwerk gemäß Bild 2.1.1 läßt sich auf das zu untersuchende Netzwerk ge- 
mäß Bild 2.1.4 durch folgende Vereinbarungen reduzieren: 

. Die Sternspannungen sind nicht vorhanden, d.h. U,= U, = U.=0 

. Die Sternwiderstände Rs und R, sind nicht vorhanden, d.h. Rs,=Rg,=0 

. Für die äußeren Spannungsquellen gilt die Zuordnung U,=0, Us = Uso, U3=- Uzo 


. Für die äußeren Widerstände gilt die Zuordnung R, =» (1° 10+'9 gesetzt), RR=R;,, 
R3=Ri2, Ra = Ra 


Son — 
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c) 
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Bei den vom Rechner ausgegebenen Ergebnissen gilt folgende Zuordnung: 

Der Strom |; ist in der gegebenen Schaltung nicht vorhanden. 

Der Strom I, entspricht dem Strom I; in der gegebenen Schaltung. 

Der Strom I; entspricht dem negativen Strom I, in der gegebenen Schaltung. 

Die Spannung U, entspricht der Spannung Ua in der gegebenen Schaltung. 

Die Leistung P, ist nicht vorhanden. 

R, kann theoretisch gleich unendlich oder auch gleich Null gesetzt werden. Programm- 
technisch sind jedoch die Extremwerte unzulässig, da diese zu einem unendlich großen oder 
unendlich kleinen Wert der Determinante führen. Daher ist R} entweder gleich 10! als 
Ersatz für unendlich oder gleich 10°! 2 als Ersatz für null zu setzen. 


Eingaben: EEX 10 STO1 


1 sto2 

2 sSTO3 
10 sTo4 
0 sTto5 

0 sSTO6 

0 sSTOA 
12 sToB 
-12 sSTOC 


Start [E] 


Ergebnisse: 1. Kurze Pause: Determinante D = 3,2 10"! 

2. Lange Pause bzw. Ausdruck: I,=0,00 A entfällt 
3. Lange Pause bzw. Ausdruck: 12=0,755A =1 

4. Lange Pause bzw. Ausdruck: I3=-0,38A=-1, 
5. Lange Pause bzw. Ausdruck: U, = 11,25 V = Ug 


6. Lange Pause bzw. Ausdruck: P, = 0,00 W entfällt 
Die Rechenzeit beträgt rd. 40 Sekunden 
Für EEX CHS 10 in STO 1 ist der Wert der Determinante D = 3,2:10°°. 


Für R,= 10°Q ergibt sich für die Determinante D = 3,2108 und nachfolgend die vorge- 
nannten Lösungswerte. Daraus ist ersichtlich, daß die gesetzten Extremwerte für R, kein 
„Umkippen” der Rechnung bewirken. 


Zu folgender Matrix A ist die Determinante D und die inverse Matrix A”! zu bestimmen: 


6 2 -1 
A=|5 4-2 
-3 1-4 
Eingaben: 6 STOI 
5 STO2 1. Spalte 
-3 STO3 
2 STO4 
4 STOS5 2. Spalte 
1 STO6 
-1 ST07 
-2 STO8 3. Spalte 
-4 STO9 


d) 


Start 


Ergebnis: D=-49 (Determinante der Matrix A) 
Rechenzeit rd. 5 Sekunden 


In einer kurzen Pause wird der Wert der Determinante mit D = -49 angezeigt. Nach Stop 
mit O in der Anzeige ist die inverse Matrix berechnet. Die Elemente sind spaltenweise in den 
Speichern STO 1 bis STO 9 abgespeichert. 


Start 


0,29 -0,14 0 sTOo1I STO4 STO7 
Ergebnis: A='=|-0,53 0,55 -0,14| = |STO2 STO5 STO8 
-0,35 0,24 -0,29 sto3 STO6 STO9Y 


Rechenzeit rd. 20 Sekunden 
Die Ausgangsmatrix A kann durch nochmaligen Aufruf der Inversion mit Start B wieder 


erzeugt werden. 


Die quadratische Matrix A soll mit der Spaltenmatrix B multipliziert werden: 
5 32 2 
A=|-1 4 0 B= 5 
6 -7 -3 1 
Gesucht ist die Spaltenmatrix C=A-B. 
Eingaben: 5 STO1I 


-1 STO2 1. Spalte 
6 STO3 
3 STOo4 
4 STOS 2. Spalte 
-7 STO6 


-3 STO9 


2 STO7? 
0 STO8 3. Spalte 
2 STOA | 


5 STOB Spaltenmatrix 
-1 STOC 
Start 
23 (abgespeichert in STO D) 
Ergebnisse: C'=| 18 (abgespeichert in STO E) 
20 (abgespeichert in STO |) 


Rechenzeit rd. 20 Sekunden 
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2.1.5 Netzwerke mit komplexen Widerständen 


Bei der Berechnung von Netzwerken mit komplexen Widerständen tritt neben den vier Grund- 
rechenarten in Polar- und kartesischer Koordinationsdarstellung als weitere Grundoperation die 
Berechnung des Parallel-Ersatzwiderstandes auf: 

Zı 

Zıt 2 


Z= 


In dem nachfolgend angegebenen Programm „Komplexe Netzwerkberechnung‘, siehe Tabelle 2.1.2, 
mit 107 Anweisungen sind die fünf Grundoperationen Addition, Subtraktion, Multiplikation, Divi- 
sion und Parallelschaltung zur Behandlung komplexer Netzwerke in Polarkoordinaten über die Start- 
adressen Label A bis Label E sowie in kartesische Koordinaten über Label a bis Label e abrufbar. 
Die Unterprogramme Label O bzw. Label 1 dienen der Koordinatenumwandlung für zwei komplexe 
Größen, die im Stackregister vorliegen. Dabei bleibt die Reihenfolge der beiden komplexen Größen 
im vierstufigen Stackregister unverändert. Mit dem Unterprogramm Label O wird von Polarkoordi- 


Tabelle 2.1.2 Anweisungsliste „Komplexe Netzwerkberechnung” 


Polarkoordinaten Kartische Koordinaten Unterprogramme 
0873 «#LBLO 
082 05Ba 874 >R 
062  65Ba 875 R+ 
0664 6551 876 R+ 
885 KIN 677 »R 
878 Rt 
eer 6SB8 73 Rt 
088 656 888 RTN 
889 6521 881 sLBLi 
818 RTN 882 >P 
083 Rd 
612 ker 864 R+ 
213 Rt 885 3P 
did x 850 xy 886 Rt 
815 Rd 851 087 Rt 
ale + 852 886 RTN 
87 Rr 852 889 sLBL2 
05 RN 636 5T03 
AU:B8 818 »LBLO @55 831 R+ 
626 E37 856 Ö 892  5$TO2 
621 Ri 857 6 893 Ri 
a2E * 858 894 ST 
023 Ry 835 R+ 
024 Kr 866 83e  STO8 
825 - 861 897 R+ 
826 Rt 062 626 RTN 
Ber KIN 863 899 xLBL3 
ANB 025 ALBLE 108  STO4 
822 6SB2 772 161 Kay 
636  65BR 182  STUS 
831 65863 1853  RCLB 
832 6SBC 1894 RELI 
0833  RÜLS 1865  RCELZ 
8354 RÜL4 186° RCL3 
035  65BO 187 RTN 
a36 RTh 872 RTN 


8] 
oO 


naten in kartesische Koordinaten umgewandelt, mit dem Unterprogramm Label 1 wird die umge- 

kehrte Umwandlung getätigt. 

Das Programm „Komplexe Netzwerkberechnung” ist kein problemorientiertes Lösungsprogramm, 
sondern es stellt lediglich die für die komplexe Rechnung erforderliche komplexe Funktion durch 
spezifizierte Unterprogrammaufrufe zur Verfügung. Bei einer konkreten Aufgabenstellung können 
diese Unterprogramme als Programmbausteine eingebaut werden [7, 8]. 


2.1.6 Anwendungsbeispiel einer komplexen Netzwerkberechnung 


Ein Netzwerk aus komplexen Widerständen (Bild 2.1.5) ist an eine Wechselspannung von U =220 V 
angeschlossen. Mit Hilfe des Programms ‚„‚Komplexe Netzwerkberechnung” sollen der Gesamtwider- 
stand sowie alle Teilströme und Teilspannungen bestimmt werden. 


Bild 2.1.5 Netzwerk mit komplexen Widerständen 


Eingabe Ergebnisse / Bemerkungen 


, [100 
=|oa 

-100 
Zu: | 200 


Parallel: Zul Zs [le] Za = (50 +j 100) 0 
ENTER 


) 
= |o Bee 
Addition: Z4 + Zz [a] Z3 = (100 +j 100) Q 


. | 100 ENTER 
2: | 50 
Parallel: Zu 112, [t][e] Z5 = (36 +52) 0 
2... | 50 ENTER 
1" | 100 


Eingangsimpedanz: 
Addition: Z3 + Z, [r][e] Z. = (136 + 2) 02 (Eingangsimpedanz) 
Ze = 136,01 80.8 g) 
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U: 


Ic 


22 


Eingabe 


0 
220 


Division: 


INJıc 


| 50 
100 


Multiplikation: I11'Z 


| 0 
220 


Subtraktion: U-U, 


, | 100 
50 
U, 
Division: —— 
22 
| -0,84 
1,62 


Subtraktion: Iı - I» 
| 0 
"(50 

Multiplikation: I2°Za 
ö | 54,46 


102,30 


Subtraktion: U2 - Uz 


„j 100 

200 
Division: — 
Za 


Bemerkungen / Ergebnisse 


Zweimaliges Rollen, damit die Impedanz Z, 
im X-, Y-Register vorliegt 


Eingangsstrom: A 
14 = 1,62 19.84 A (Eingangsstrom) 


Polarkoordinatenumwandlung 


U} = 180,84 e 1274" y 


Zweimaliges Rollen, damit der Minuend U, 
im X-, Y-Register vorliegt. 


U, = 102,30 d) 54.46° 


Polarkoordinatenumwandlung 


12=0,91 389” q 


Zweimaliges Rollen, damit der Minuend 1, 
im X-, Y-Register vorliegt. 


13=0,73 094° a 


U, = 36,29 e) 942° y 


Zweimaliges Rollen, damit der Minuend U; 
im X-, Y-Register vorliegt. 


U, = 80,84 e) 72.9 y 


Potarkoordinatenumwandlung 


14= 0,36 09-5 a 


Eingabe Tasten Ergebnisse / Bemerkungen 


-100 
2: | 200 


Polarkoordinatenumwandlung 


Multiplikation: 1424 Us = 80,84 e) 72:9 y 


100 
25: | 0 


Polarkoordinatenumwandlung 


15=0,81 e117.02° a 


2.2 Übergangsvorgänge 1. Ordnung 


Die Sprungantwort von Systemen mit einem Energiespeicher kann bei bekanntem Anfangswert 
und Endwert durch eine Übergangsfunktion 1. Ordnung beschrieben werden: 


t t 


xit)=x0e T+x, (1-6) (2.2.1) 


Die Übergangsfunktion enthält eine mit der Systemzeitkonstanten T abklingende Anfangsgröße xo 
und die sich mit der Systemzeitkonstante aufbauende Lösung für den Endzustand x. 
Sind die Vorzeichen von Anfangswert xo und Endwert x. verschieden, so schneidet die Übergangs- 
funktion x(t) zur Zeit to die Nullinie. Für die Zeit ty gilt 
X %o 

Ko 


%=TIn (2.2.2) 


2.2.1 Programmbeschreibung „Übergangsvorgang 1. Ordnung” 


Das in Tabelle 2.2.1 angegebene Programm für die zyklische Berechnung der Funktionswerte aus 
GI. (2.2.1) benötigt außer dem Stackregister und den Speichern STO A bis STO D für die Eingabe- 
werte X9, X, T und Ar keine Arbeitsspeicher. Das Zyklusprogramm wird über die Startadresse 
Label A aufgerufen. Der im X-Register vorliegende Wert wird als Startwert für die normierte Zeit 
T= u genommen. Mit den drei ENTER-Anweisungen wird in allen Stackregistern die normierte 


Zeit eingeschoben. 

Mit der Print-Anweisung in Zeile 008 wird die jeweilige Zeit t ausgedruckt. Das Ergebnis der 

GI. (2.2.1) wird in Zeile 022 ausgegeben. In den Zeilen 025, 026 wird die normierte Zeit r um 

die in STO D vorgegebene Schrittweite Ar erhöht und mit dem Sprung an den Programmanfang 
der nächste Berechnungszyklus eingeleitet. 

Die Berechnung der Zeit to für den Nulldurchgang der Funktion x(t) gemäß GI. (2.2.2) wird in 

den folgenden Programmteil beginnend mit Label a in Zeile 028 bis zur R/S-Anweisung in Zeile 038 
durchgeführt. In Zeile 035 wird in einer kurzen Pause die normierte Zeit 7. angezeigt. 
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Tabelle 2.2.1 Anweisungsliste „Übergangsvorgang 1. Ordnung” 


Zeitfunktiin_ @6l »LöLA 020 x 
802  ENTt 021 + 
aa3  ENTt 822 PRIX Dxtt) 
ea4  ENTT 023 SPC 
665 Keil 824 xer 
BE x 825  RELD 
867 PSE<FD], e26 + 
@u8 FPIX > 827__6TOR 
80° Ru 
Blö  RÜLA 829 RCLB 
ei er 03a RELR 
[36 ChS 831 - 
812 e 032  RCLB 
818 x 833 ® 
el5 LET 834 IN, 
ale ChS 835 PSE DI 
air 1 836  RELE 
als + 037 x ; 
615 RÜB 836 RS Dh 


2.2.2 Anwendungsbeispiel 


In der in Bild 2.2.1 angegebenen Schaltung wird 


zur Zeit t=0 der Schalter S geschlossen. Gesucht a = R 
ist der zeitliche Verlauf der Spannung uc {t). Ug N S 4 " 
T = R-C= 10° 5000  10°CF =5s 
Uco = -10V Bild 2.2.1 Einschaltvorgang 
U. 12 v eines RC-Kreises 
Uc () = -10V ar Er +12 V ( .. ) Zahlenwerte: 
Uo = 12V 
R=-1KN 

Eingaben: Anfangswert Uco:- 10 STOA C = 5000 »F 

Endwert Uco: 12 STOB uc(t=0)=-10V 

Zeitkonstante T: 5 STOC 

normierte Zykluszeit Ar: ı STOD 
Start: 0, 
Ergebnisse: tı = 0.00 s tg = 25.005 

uc(tı)= -10.00 V uc (te) = 11.85 V 

ta = 5.005 tz = 30.00 s 

uc (t2) = 3.91 V uc{t7)= 11.95 V 

ta = 10.005 tg = 35.00 

uc (ta) = 9.02 V Uc (te) = 11.98 V 

t4 bu 15.00 s tg = 40.008 

uc(ts)= 10.90 V uUc (te) = 11.99 V 

t5 = 20.005 

Uclts)= 11.60 V 


Die Rechenzeit für einen Zyklus beträgt rd. 3 Sekunden. 
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Die Zeit für den Nulldurchgang wird nach dem Aufruf des Programms in Label a ausgegeben. 


Start: [2] 


Ergebnis: To = 0,61 (Anzeige mit kurzer Pause) 
to - 3.03 Ss 


Der zeitliche Verlauf der Spannung uc {t} 
ist in Bild 2.2.2 aufgetragen. 


Bild 2.2.2 
Übergangsvorgang 1. Ordnung 


0 45 10 15 20s 25 


—t 


2.3 Freie Schwingung in elektrischen Stromkreisen 
2.3.1 Einschaltvorgang eines RLC-Kreises 


Der zeitliche Verlauf der Ströme und Spannungen in einem elektrischen Stromkreis nach einem 
Schaltvorgang läßt sich mit dem Überlagerungsgesetz aus der Summe eines stationären Anteils und 
eines Ausgleichanteiles bestimmen. Der Ausgleichanteil stellt einen freien Schwingungsvorgang in 
dem Stromkreis dar, durch den der vorgegebene Anfangszustand auf den von außen eingeprägten 
stationären Endzustand übergeführt wird [9]. 


Für den Gesamtstrom i in Bild 2.3.1 mit 


iFigti (2.3.1) 
gilt die Maschengleichung 
i dia ! fiad un +2 [ia (2.3.2) 
= + + ES 3 
usigR+tL dt +7 iadt ti Fra sdt 


it 


1 so IL: Bild 2.3.1 


Einschaltvorgang eines RLC-Kreises 


Setzt man die Summe der Terme mit dem freien Strom is; gleich Null, so verbleibt die für den 
stationären Endzustand gültige Gleichung übrig. Es gilt also unabhängig von der äußeren Spannung 
die Gleichung 

de 17 

— + | iydt (2.3.3) 


=4Rt+ 
O=%R LG c. 
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Beim Einschaltvorgang stellt die Gl. (2.3.3) nur einen Teil der Lösung von GI. (2.3.2) dar, während 
sie beim Ausschaltvorgang (Kurzschluß der Spannungsquelle) die gesamte Lösung liefert. 

Die Differentation der Gl. (2.3.2) liefert eine homogene Differentialgleichung zweiter Ordnung, 
ersten Grades der Form: 


u Rd 1 f. 
ren tajüro (2.3.4) 


Mit der allgemeinen Lösung: 


4ece Tre?‘ (2.3.5) 


R 1 R® L_ 
ea tr vie) rat (2.3.6) 
REN 200 vey-L- e 
bee Az: (2.3.7) 


L 
LT. 

FOREN 2 fuer (2.3.8) 
1[Cı ne, ©, 

ut lBerr Ze) wa 


Da die magnetische Energie der Induktivität proportional dem Quadrat des Stromes und die elek- 
trische Energie der Kapazität proportional dem Quadrat der Spannung ist, bleiben Strom und alle 
Spannungen im Schaltaugenblick unverändert, d.h. als Anfangsbedingung gelten diejenigen Werte, 
mit denen der vorhergehende Dauerzustand beendet wurde: 


i,(+0) = igl-0) = Io (2.3.10) 
Uct (+ 0) =Ug (- 0) — Uco (2.3.11) 


2.3.2 Anfangsbedingungen 


Charakteristische Übergangsvorgänge ergeben sich für zwei verschiedene Anfangsbedingungen: 


a) Die freie Schwingung wird eingeleitet, wenn die Aufladung des Kondensators auf + Uco 
beendet ist, d.h. der Strom ig, = 0 wird. Der Ausgleichsvorgang bewirkt zunächst einen 
Abbau der im Kondensator gespeicherten Energie. 

b) Die freie Schwingung wird eingeleitet, wenn der Stromfluß durch die Induktivität sein 
Maximum iso = Io erreicht, d.h. die Spannung u> = O wird. Der Ausgleichsvorgang bewirkt 
zunächst einen Abbau der gespeicherten magnetischen Energie. 


2.3.3 Aufgeladener elektrischer Energiespeicher 


Der Fall des aufgeladenen elektrischen Energiespeichers entspricht den Anfangsbedingungen 
uc+H (0) = Uco und is(0) =0. 
Aus Gin. (2.3.3) und (2.3.7) folgt für t=0: 


rn "fa 
&1=-02=C Us d0=m, (2.3.12) 


-n 
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Damit ergibt sich für den freien Strom is aus Gi. (2.3.3) 


U dr _ „-6t 
ig = 1; are —— (2.3.13) 
oder in hyperbolischer Schreibweise: 
U 
ig = Tr e*t sinh öt (2.3.132) 
Aus Gl. (2.3.7) ergibt sich für die freie Einschwingspannung uc+ am Kondensator: 
Uco „= ( et et et+ dt ) 
ua=ye (a 2 +6 2 (2.3.14) 
oder in hyberbolischer Schreibweise: 
Uco at f 
Uc= ze (a sinhöt +8 coshöt) (2.3.14a) 


Vor Inangriffnahme der numerischen Lösung der Gin. (2.3.13 und 2.3.14) sind drei Fälle zu unter- 
scheiden, die zu den drei charakteristischen Ausgleichsvorgängen: aperiodischer Verlauf, aperiodi- 
scher Grenzfall und gedämpfte periodische Schwingung führen: 


1. Aperiodischer Verlauf, d.h. 8 = reell 


Es gilt: 


R L 
27 Vs (2.3.15) 


Für die Programmierung geeignete Formen der Lösungsgleichungen sind: 


. Uco -7 st _ dt 
ke ee . (2.3.16) 
Uco a ed! ee ) 
“etc AR-AB Rr2l RA 
mit 

2L 

T=7 (2.3.18) 
1 R\ L 

=} (3) E (2.3.19) 


2. Aperiodischer Grenzfall, d.h. &=0 
Es gilt: 
R L 
ae y: (2.3.20) 


Da die Lösungsgleichungen (2.3.16 und 2.3.17) mathematisch unbestimmte Formen ergeben, 
müssen diese nach der L’Hospitalschen Regel bearbeitet werden: 
Wegen: 


=2t (2.3.21) 
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ergeben sich für die Programmierung geeignete Lösungsformeln zu: 


2U ade 
u Ze "= Tete T (2.3.22) 
_t 
Ge = Ur, (1+5)e T (2.3.23) 


3. Gedämpft periodische Schwingungen, d.h. 5 = jw 
Es gilt: 
R.,/t 
2<Vc 
Aus Gin. (2.3.16 und 2.3.17) ergeben sich mit der Eulerschen Formel @*!=coswt +j sinwt für 
die Programmierung geeignete Lösungsgleichungen zu: 


Uco — 
= ——e sinwt (2.3.24) 
wL 
Uco /L -5 ; 
U V ge in (wt +arctanwT) (2.3.25) 


2.3.4 Erregter magnetischer Energiespeicher 


Der Fall des erregten magnetischen Energiespeichers ist besonders wichtig bei den Untersuchungen 
von Überspannungen an elektrischen Betriebsmitteln als Folge von Kurzschlußabschaltungen. Dieser 
Fall entspricht den Anfangsbedingungen iz (0) = I, und uc+(0) =0. 

Aus Gin. {2.3.5 und 2.3.9) folgt für t=0: 


Fr 
lo (2.3.26) 
f2 
= lo, =r, (2.3.27) 
Damit ergibt sich für den freien Strom iz aus Gl. (2.3.5): 
St —6t st dt 
. „0 ar ( Case N a Ba ) 
= get (a 2 ö 2 (2.3.28) 
Aus Gl. (2.3.9) ergibt sich für die freie Einschwingspannung uc+ am Kondensator: 
I st _ „dt 
uc=- = get > (2.3.29) 


Entsprechend zum vorigen Abschnitt sind wieder drei Fälle für den Ausgleichsvorgang zu unter- 
scheiden: 


1. Aperiodischer Verlauf, d.h. 6 = reell 


I 
Ucı = 705° T(eöt_.endt) (2.3.30) 
._lo Fir &_[R 51 

u Be (&+5)e -(&-5) (2.3.31) 
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2. Aperiodischer Grenzfall,d.h. &=0 


BE oT 2333 

U4=77e „te (2.3.32) 
ot 

ig = lo (1-#)e F (2.3.33) 


3. Gedämpft periodische Schwingung, d.h. d = jw 


t 


a =& 
Ur = - zE e Tsinwt (2.3.34) 
lb [U -# 
4„=- — V se sin(wt-arctan@wT) (2.3.35) 
wLYycC 


2.3.5 Programmbeschreibung ‚Freie Schwingung C” 


Das in Tabelle 2.3.1 angegebene Programm umfaßt 223 Programmzeilen. Vor dem Programmstart 
über die Startadresse Label A sind folgende Eingaben erforderlich: 


Anfangsspannung des Kondensators Uc. in V>STOA 
Widerstand des Reihenschwingkreises R in 2=>STOB 
Kapazität des Reihenschwingkreises C in F»>STOC 
Induktivität des Reihenschwingkreises L in H > STOD 


Nach einigen vorbereitenden Rechnungen und Speicherbelegungen in den Zeilen 003 bis 025 wird 
als erstes Ergebnis die Zeitkonstante T gemäß Gl. (2.3.18) in Zeile 026 ausgegeben. Mit den Ab- 
fragen in den Zeilen 028 und 033 werden die Fallunterscheidungen für die Varianten der möglichen 
Lösungen: aperiodischer Verlauf, aperiodischer Grenzfall und gedämpft periodische Schwingung 
vorgenommen. Die Gleichungen für eine gedämpft periodische Schwingung Gin. (2.3.24 und 2.3.25) 
werden ab Adresse Label a in Zeile 124 ausgewertet. Beim aperiodischen Grenzfall, d.h. der Wurzel- 
radikant aus Gl. (2.3.19) verschwindet, wird die weitere Berechnung bei der Adresse Label b in 
Zeile 093 fortgesetzt. Der aperiodische Verlauf wird ab der Programmzeile 035 bis zu der Zeile 092 
berechnet. Für die beiden aperiodischen Berechnungsvarianten wird in den Zeilen 030, 031 die 
Schrittweite für die zyklische Berechnung auf 1/5tel der Zeitkonstanten T festgelegt und in STO 8 
abgespeichert. 

Für den Fall der gedämpft periodischen Schwingung werden Lösungspunkte an den Maximastellen 
der Sinusfunktion der Gin. (2.3.24 und 2.3.25) berechnet. Hierzu wird eine Umsteuerung mit Hilfe 
des Flags F2 zwischen dem Wert von arctan«wT + k- 180° und dem Wert 90° + k : 180° vorge- 
nommen. Der Wert arctanwT ist als Winkel im Gradmaß in STO 3 abgespeichert. 

Die Schleife für die zyklische Berechnung beginnt für den gedämpften Schwingungsfall mit Label c 
in Zeile 168, für den aperiodischen Grenzfall mit Label d in Zeile 102 und für den aperiodischen 
Fall mit Label e in Zeile 052. Für alle Varianten wird die aktuelle Zykluszeit t (bzw. der Winkel wt 
für den Fall der gedämpften Schwingung) in einer kurzen Pause angezeigt. Hierzu ist das Unterpro- 
gramm Label 3 von Zeile 218 bis 223 vorgesehen. Dadurch besteht die Möglichkeit, während der 
kurzen Pause durch Betätigung der R/S-Taste die Berechnung anzuhalten und eine spezielle Zeit 
oder auch eine spezielle Schrittweite At in STO 8 einzugeben. Der im X-Register vorliegende Wert 
wird bei wiederholter Betätigung der R/S-Taste als Zeit t oder als Winkel wt ausgedruckt und in 
den Speicher für die laufende Zeit t (bzw. für den Winkel wt) übernommen. 
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Tabelle 2.3.1 Anweisungsliste ‚Freie Schwingung C’ 


Starf ucg 
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del _*LBiR 

“a2 CF2 

a63 1 

664 0 

885  STOI 

@uE  RÜLD 

807 STOB 

eos STU2 

089 STOZ 

e18 RÜLC 

el ST:8 

eı2 * 

813  RCLB 

814  STO6 

815  $T07 

el6 ST=2 

017 2 

618  ST0S 

815 STx2 

820  5STx3 

al * 

022 x2 

023 - 

024 REL2 

825  STOß 

26 PRIS=DT 

827 R+ 

6285 %8? 

823  6T0a 885 
058 5 EIS 
831 STB 887 
832 R 888 
0855  X=8? 085 


056 * 034 
0833  STOI 895 
84a STx3 e36 
841 1 837 
842 6588 [E15 
843 RCLA 099 
0844 RCL3 188 
645 * 181 
846  STO4 182 
647 RCLE 103 
848 x 104 
843 STx5 165 
058 6 166 
0851  5T0i 167 
052 «LBLe 108 
0535 6583 162 
054 RELB 116 
055 RCL3 1il 
e56 = 112 


834  ETOb 098 
835 CHS es 
0836 X 882 
637 RCLD 833 


113 
114 
115 
116 
117 
118 
113 
128 
121 

122 
123 
124 
125 
126 
127 
128 
123 
138 
131 

132 
133 
134 
135 
156 
137 
138 
139 
148 
141 

142 
143 
144 
145 
146 
147 
148 
149 
158 
151 

152 
153 
154 
155 
156 
157 
158 
159 
168 
161 

162 
163 
164 
15 
166 
167 
168 


1837 


STE 

PRTX=D fo 
17% 

PRTX=D To 
2 


STO8 
RCL2 
RELI 
x 
TAN“ 
5TO3 


K) 
6588 
y 

e 

CHs 
5T0i 
#LBLi 
RCL3 
ST-i 


up 


Up 


UP 


Das Unterprogramm Label O ab Zeile 200 dient zur Erzeugung eines Merkmales zur Kennzeichnung 
der jeweils vorliegenden Schwingkreisvariante. Für den aperiodischen Fall wird eine Zeile 111... 
erzeugt, für den aperiodischen Grenzfall eine Zeile 222... und für den Fall der gedämpften Schwin- 
gung eine Zeile 333... 

Im Fall der gedämpft periodischen Schwingung wird zusätzlich zu der Zeitkonstanten T noch die 
Frequenz fo und die Periodendauer T, dieser Schwingung ausgegeben. Die zyklische Ergebnisaus- 
gabe erfolgt in Dreiergruppen als Zeit t bzw. Winkel wt, Strom is der freien Schwingung und 
Spannung am Kondensator uc+ der freien Schwingung innerhalb einer Endlosschleife, die durch 
Betätigung der R/S-Taste abgebrochen werden kann. 


2.3.6 Test- und Anwendungsbeispiele 


Es soll die freie Schwingung eines Reihenschwingkreises mit folgenden Daten ermittelt werden: 


Anfangsspannung des Kondensators Uc. = 1.000V, Kapazität C = 0,3uF, 
Induktivität L = 0,0844 H 


Widerstand a) R = 500 bR=2 y: c) R = 2.0002 


a) Eingaben: 1.000 STOA (Uco) 
50 sSTOB {R) 
0,3EEX- 6 STOC (C) 
0,0844 STOD (L) 


Start [A]: 
Ergebnisse: T=3.376 -10°°s 
fo = 999,09 Hz 
To = 1,0009 - 10”°s 
3333333333 Merkmal für gedämpfte Schwingung 


y= 2,70° „= 182,70° 

ir =0,09 A 1, =-0,08 A 
uc+ = 998,89 V Ucı = - 861,27 V 

„= 90,00° »= 270,00° 

ig = 1,75 A = -1,51 A 
uc+ = 43,81 V ucı =- 37,78 V 


Die Ergebnisse sind in Bild 2.3.2 dargestellt. 


DIL7 
r\r 


Bild 2.3.2 

Freie Schwingung beim aufgeladenen 
elektrischen Energiespeicher für den Fall 
einer periodisch gedämpften Schwingung 
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1.000 
03EEX- 6 
0,0844 


RCLD,RCLC, +, YX, 2, x 


Start [a]: 


Ergebnisse: T = 1,59...10” *s 
2222222222 


b) Eingaben: 


t=0s 
=0 

uc+ = 1.000 V 
t=3,18...10 °s 
ie=0,31 A 

ucr = 982,48 V 


Nach Anzeige t = 6,36 10°s Betätigung der 
R/S-Taste zur Änderung der Zykluszeit auf 1 ms: 


.IEEX —-3 sTO8 
.2EEX -3 R/S 


t=20-10 *s 
+ =0,67 A 
Uc+ = 642,17 V 
t=3,0:10 *s 
r=0,54 A 
Uc+ = 437,93 V 


Die Ergebnisse sind in Bild 2.3.3 mit den 
Kurven a) für den Sprung und.) für den 
Strom aufgezeichnet. 


ng 


c) Eingaben: 1.000 STOA (Uco} 
2.000 STOB {R) 
03EEX-6 STOC (C) 
0,0844 STOD {L) 


Merkmale für den aperiodischen Grenzfall 


1,0 
A 


| 0,6 


0,4 


0,2 


0 0 Odms 10 


0 022 04 086 


—t 


Bild 2.3.3 Freie Schwingung beim auf- 
geladenen elektrischen Energiespeicher 
für den aperiodischen Fall 


a) Spannung beim aperiodischen Grenzfall 


L 
a2yE 
b) Spannung beim allgemeinen aperio- 


dischen Fall R = 2000 2 
c) Strom beim aperiodischen Grenzfall 


„/L 
a-ayl 


d) Strom beim allgemeinen aperiodischen 


Fall R= 20002 


Start [A]: 


Ergebnisse: T =8,44 10" °s 
1131111111 Merkmal für den aperiodischen Fall 


t=00s 
4=0,0A 2. 
uc+ = 1.000 V h = 3,376 : 10” °s 
N -5 4=027A 
® = 1,688 10° °s ucı = 982,53 V 
4=0,16 A 
ucr = 995,06 V 


Die Ergebnisse sind in Bild 2.3.3 mit den Kurven b) für die Spannung und d) für den Strom 
aufgezeichnet. 
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Tabelle 2.3.2 Anweisungsliste „Freie Schwingung L” 


Start ig 881 
802 
063 
284 
805 
Bas 
607 
808 
869 
818 
el 
B1z 
@13 
6id4 
815 
eis 
#17 
e1& 
&19 
82a 
821 
622 
823 
024 
625 
8265 
827 
628 
623 
838 
831 
832 
833 
&34 
835 
83E 
837 
838 
039 
846 
041 
842 
043 
344 
845 
B46 
647 
048 
849 
858 
851 
852 
053 
854 
855 
856 


*LBLH 113 
CF2 114 
1 115 
[) 116 
5Tal 117 
RÜLD 118 
5T08 119 
5TÜR 128 
RCLE 121 
5TÜ3 122 
sT+8 123 
* 124 
RÜLE 125 
STüe 126 
star 127 
ST#z 126 
& 125 
STX2 138 
STx3 131 
® 132 
5Tüs 133 
a8 134 
- 135 
5T0s 136 
RÜLE 137 
sToR 135 
PRTA=DT 133 
Ry 146 
ar: 141 
&Tla 142 
Ei 143 
ST+& 144 
R+ 145 
a=07 146 
ETüb ——e 891 HLBLb 147 
CHS 832 £ 148 
Y& 8983 65Bö 149 
RCLE 894 RELR 158 
* 895  5Tx4 15 
sTül 836 8 152 
STx3 8987 5STüs 153 
1 098 xLBLd 154 
65B& 8992 6583 155 
RÜLA 188 RCLZ 156 
RCLS 161 * 157 
+ 182  STOS 158 
5T04 103 STÜS 159 
RÜLC 184 RCLi 168 
x 185  6SB2 161 
$T05 186  STAS 162 
8 187 ST-5 163 
ST0i 188  RCL4 164 
*LBLe 189 x 165 
6563 118  RCL3 166 
RCLI 111 x 167 
RCLZ 112 CHS 168 


PRTSDT, 
€ 


sTüs 


cHs 
PRTXDucf 
RCL4 
RCLO 


Up 


UP 
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2.3.7 Programmbeschreibung „Freie Schwingung L” 


Das Programm dient zur Beschreibung der freien Schwingung in einem RLC-Reihenschwingkreis 
mit einem endlichen Strom als Anfangsbedingung (erregter magnetischer Energiespeicher) und 
energielosem Anfangszustand des Kondensators. Das in Tabelle 2.3.2 angegebene Programm umfaßt 
224 Programmzeilen. Vor dem Programmstart über die Startadresse Label A sind folgende Eingaben 
erforderlich: 

> STOA 
Widerstand des Reihenschwingkreises R in >» STOB 
Kapazität des Reihenschwingkreises CinF >STOC 
Induktivität des Reihenschwingkreises L inH > STOD 


Anfangsstrom der Induktivität Io in A 


Der prinzipielle Aufbau des Programms entspricht demjenigen nach Tabelle 2.3.1 für den aufge- 
ladenen kapazitiven Energiespeicher als Anfangsbedingung. Die zyklische Berechnung der Konden- 
satorspannung uc+ und des Kreisstromes i; wird nach den Gin. (2.3.30 bis 2.3.35) durchgeführt. 
Ein zusätzliches Unterprogramm Label 4 dient zur Ausgabe des Stromes für die beiden aperiodi- 
schen Lösungsvarianten. 


2.3.8 Anwendungsbeispiel 


Berechnung der freien Schwingung nach Kurzschlußabschaltung einer Hochspannungsleitung. 

Eine 110-kV-Leitung werde in 50 km Entfernung von einer Kurzschlußstelle durch einen Leistungs- 
schalter abgeschaltet. Der Lichtbogen erlischt bei einem Reststrom von 10 A. Es soll der zeitliche 
Verlauf der Spannung an der Abschaltstelle anhand einer r-Ersatzschaltung in Abhängigkeit von 
der Zeit berechnet werden. 


Leitungsdaten: R’= 0,2 2/km Xc = 300 KR: km 
X, =0,4.Q2/km I=50 km 


Der Kurzschlußwiderstand an der Fehlerstelle sei 902. 
Der Ersatzschaltplan ist in Bild 2.3.4 dargestellt. 


Io t:0 ir R= 


i iu a Bild 2.3.4 
| ah Ersatzschaitplan der Leitung 
-z ox: = 0,265 uF mit Kurzschluß am Leitungsende 
Eingaben: 10 STOA (Anfangsstrom iso) 
100 STOB (Widerstand R) 
0,267 EEX — 6 STOC (Kapazität C) 
63,7 EEX — 3 STOD (Induktivität L) 


Start [A]: 


Ergebnisse: Zeitkonstante T = 1,274 ms 
Frequenz der freien Schwingung fo = 1218,59 Hz 
Periodendauer der freien Schwingung To = 0,82 ms 
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Merkmal für die Lösungsvariante einer gedämpft 


periodischen Schwingung: 3333333333 


wt= 90,00° 
ucı =-4195,48 V 
4 =-0,87 A 
wt= 174,15° 
uc+ = - 368,06 V 
= -7,36 A 
wt= 270,00° 
Uct = 3040,26 V 
= 0863A 
wt= 354,15° 
uUct = 266,71 V 
= 5,33 A 


Yr Uef 


0 0 
-5 -2500 
-10 -5000 t 
0 360 70 —- wt 1440° 
0 0,82 1,64 —t 3,28ms 


Bild 2.3.5 Freie Schwingung beim erregten magne- 
tischen Energiespeicher für den Fall einer periodisch 
gedämpften Schwingung 


In Bild 2.3.5 sind die freien Schwingungen der Kondensatorspannung uc+ und des Kreisstromes ir 


in Abhängigkeit von der Zeit aufgetragen. 


Das Isoliervermögen der beiden Schalterpole muß der Überlagerungsspannung aus der anliegenden 
50-Hz-Wechselspannung und der Kondensatorspannung u.; mit einer Frequenz von ca. 1,2 kHz 


standhalten. Die Amplitude der Kondensator- 


spannung ist entscheidend abhängig von der ausge- = H 
blasenen Stromamplitude während der Kurzschluß- 
stromabschaltung und Lichtbogenlöschung im uf 100 i£ 5 
Schalter Kuren = imax V:). Der ohmsche , N 
Widerstand an der Kurzschlußstelle beeinflußt über 
die Abklingzeitkonstante der hochfrequenten 
Überlagerungsspannung nur indirekt deren SUR FR 
Amplitude und damit die Spannungsbean- 
spruchung des Schalters. Bild 2.3.6 zeigt den -2000 -10 
Verlauf für den aperiodischen Grenzfall bei 

Bild 2.3.6 


n 
R=2/ C= 980,570. 


2.4 Netzwerk-Umwandlung 


m 


Freie Schwingung beim erregten magnetischen 
Energiespeicher für den aperiodischen Grenz- 


fall mi = — 
allmit R 2y& 


L 


Eine grundlegende Aufgabe bei Netzwerk-Umwandlungen ist die Umwandlung von Dreieck-Schaltun- 
gen in äquivalente Stern-Schaltungen und die Umwandlung von Stern-Schaltungen in äquivalente 
Dreieck-Schaltungen (Bild 2.4.1). Die Lösung der Umwandlungsgleichungen erfordert bei komplexen 
Widerständen einen erheblichen rechentechnischen Aufwand, so daß ein programmierter Umwand- 
lungsalgorithmus eine wesentliche Hilfe erbringt [10]. 


Bild 2.4.1 


Bezeichnungen zur Netzwerkumwandlung 


a) Dreieck-Schaltung b) Stern-Schaltung 


3 
NE 


Für die komplexen Widerstände gilt in kartesischer Koordinatenschreibweise: 
Dreieck-Schaltung: 


Zıa=RıatiXı2; Zas=RastiXas; Zi RsıtiXa (2.4.1a,b, c) 


Stern-Schaltung: 
Zı=Rı+iX4; Za=RoatiXa; Za=RatiXa (2.4.2a,b, c) 


2.4.1 Stern-Dreieck-Umwandlung (\ = A-Umwandlung) 


Die Impedanzen einer Dreieckschaltung ergeben sich aus den Stern-Impedanzen zu: 


_EıZ2atZaZatZ3Zı a BudElpHar"X4 Ko) + RalHi Rt Ba] 
e12 7 2 er 2 2 2 
au Rat X (2.4.32) 
Pe FRE Ra . 
J 1 2 R2+xX2 a 
a+ ‘3 
_ZıZ2tZoZ23t232ı _ Rı (RaRa = X2Xa) + Xı (R2Xa + Ra%2) 
Fi 175 21 19 2 Zn ne Tee at 
> En (2.4.3b) 
Rı (R2X3 + RaX2) = Xı (RaRa — X2X3) > 
R7+X 
_Zı&atZ2Z3tZ23Z2 Rz (RaRı =X3X4) + X2(R3Xı + RıX3) 
> u es) A 7 ei; 
- rn (2.4.3c) 
sr |. % Mel FF Rah Ra (BR Re u 
J 3 ı 2 2 : 
R,+X 


2.4.2 Dreieck-Stern-Umwandlung (A > A -Umwandlung) 


Die Impedanzen einer Sternschaltung ergeben sich aus den Dreieck-Impedanzen zu: 


_  2ı2'2aı _ARzıRı2 = X31%ı2) R+ XzıRı2 + Rzı%ı2) X 
Z, = zZ +Z +Z = 2 2 f 
2127 223" £3ı R’+X 
(2.4.48) 
er (X31 Rı2 + Rz1Xı2) R= (Raı Rı2 = Xa1 Xı2) X 
R?+x2 
nn 223 Zıa _ [RıaR2a = X12%22) R + (X12R23 + Rı2X23) X E 
2 Zu2t+ZastZaı R?+x2 
(2.4.4b) 
u (X12Ra3 + R12X23) R- (Rı2aR23 — X12X23) X 
J R? + x2 
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Za1'Z2a 2 _ (RasRaı = X23Xa1} R + (X2aRaı + R2aXa1) X & 


este: 
I ZuatZastZaı R? +x? 
(2.4.4c) 
PR: (X23Ra1 + R2a3X31) R - (RaaRaı — X23%X31) X 
| R?+x2 
R= Rıa + Ras B Raı und X = X1a + X23 r Xaı (2.4.5a, b) 


2.4.3 Programmbeschreibung „Netzwerkumwandlung’’ 


a) A -, A-Umwandlung 


Das in Tabelle 2.4.1 angegebene Programm beginnt mit der Startadresse Label A für die Stern- 
Dreieck-Umwandlung in Programmzeile 001. Mit der folgenden Anweisung GSB 2 wird ein Unter- 
programm zur Berechnung des Produktes Z; ' Z, aufgerufen. Die eigentliche Produktbildung wird 
in einem weiteren Unterprogramm mit der Einsprungadresse Label c vorgenommen. Da dieses 
Unterprogramm auch zur Addition komplexer Größen mitbenutzt wird, ist eine Flag 2-Abfrage 
eingebaut, durch die bei Einstieg über Label c die in Zeile 124 vorgenommene Addition mit den 
Anweisungen in den Zeilen 127 bis 129 wieder rückgängig gemacht wird. Nach Ausführung der 
Anweisung LAST X in Zeile 129 entspricht die Konfiguration im Stack-Register wieder dem Stand, 
nach dem die Anweisung in Zeile 123 ausgeführt war. Das Flag 2 darf daher nicht vom Tastenfeld 
aus gesetzt werden. 

Das in Polarkoordinaten übergebene Ergebnis wird in kartesische Koordinaten umgewandelt und 
der Realteil in STO 7, der Imaginärteil in STO 8 abgespeichert. Mit den folgenden Anweisungen 
007 bis 016 werden auf entsprechende Weise die Produkte Z> ' Zy3 und Z3' Z, gebildet und je- 
weils den Inhalten von STO 7,STO 8 aufaddiert. In diesen Speichern steht damit der Zähler aus 
den Gin. (2.4.3a bis c) zur Verfügung. In den Zeilen 017, 018 wird die indirekte Adressierung für 
die Auswertungen der Gin. (2.4.3a bis c) im Unterprogramm Label 4 vorbereitet. Beim ersten Auf- 
ruf dieses Unterprogramms wird durch die Anweisung in den Zeilen 107 und 109 die Impedanz Z, 
als Imaginärteil in das Y-Register und als Realteil in das X-Register geladen. Nach Umwandlung in 
Polarkoordinaten wird das Unterprogramm Label d aufgerufen, um die Division gemäß GI. (2.4.3a) 
durchzuführen. 

Das Ergebnis, also die Dreieckimpedanz Z}>, wird mit Real- und Imaginärteil in den Sekundär- 
speicher STO 1’ und STO 2’ abgespeichert. 

Ab der Programmzeile 026 wird in entsprechender Weise für die Berechnung der Dreieckimpedan- 
zen Z., und Z,, verfahren. Das logische Programmende ist mit der R/S-Anweisung in Zeile 041 
erreicht. Die gesuchten Dreieckimpedanzen stehen dann in dem aktiven Speicherbereich von STO 1 
bis STO 6 jeweils als Real- und Imaginärteil bereit. Die Eingabeimpedanzen stehen nun unverändert 
im Sekundärbereich zur Verfügung. 


b) A-, A-Umwandlung 


Das Dreieck-Stern-Umwandlungsprogramm beginnt mit der Startadresse Label B in Programm- 
zeile 042. Mit den Anweisungen in den Zeilen 043 bis 054 wird durch zweimaligen Aufruf des 
Unterprogramms Label a die Summe aller Dreieckimpedanzen für die Nenner der Gin. (2.4.4a bis c) 
gebildet und nach Umwandlung in Polarkoordinaten als Betrag in STO 7 und Phasenwinkel in 
STO 8 abgespeichert. Nach dem Einsprung in das Unterprogramm bei Label a wird in Zeile 120 
das Flag 2 gesetzt. Hierdurch ist es möglich, einen Teil des Unterprogramms c mit zu benutzen. 
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Tabelle 2.4.1 Anweisungsliste „Netzwerkumwandlung” 


A+d& 


d>ı 


sul 
882 
063 
884 
885 
886 
687 
686 
883 
eia 
ei 
Blz 
813 
014 
@15 
dıe 
e17 
eı1e 
819 
028 


822 
623 
824 
825 
e26 
827 
828 
023 
838 
831 
832 
033 
834 
835 
636 
837 
838 


842 
845 
044 
845 
846 
847 
84E 
849 
058 
851 
e52 
053 
054 
055 
856 


FES 
sTol 
xEY 
sT0Z 
F2S 


Nana IE 


ALELE 


5Tu7 

Ser 
sTo8 
65Bil 
6$B5 


857 P2S 
858 | STOI 
859 er 
ese | STO2 
061 P2S 
662  65B2 
863 65885 
864 85 
865 | ST03 
866 Key 
867 | 5T04 
868 *5 
869 6553 
878  65B5 
871 P2S 
872 | STO5 
0873 Kay 
874 | STO6 
UP 876 A#LBLI 
877 RCLE 
ereE  RCLS 
879 pP 
088  RCLZ 
esi RELI 
[1:3 >F 
685  6SBe 
184 __RIN 
UP 685 sLBLZ 
866  RLLE 
e87  RCLI 
888 PP 
0869 RiL4 
e9e RELS 
891 pP 
____092  6SBe 
up 893 RTN 
0934 #LBL3 
835  RCL4 
096 RÜLS 
897 pP 
898  RCLE 
08939  RCLS 
186 >P 
181 3:7 
182 RTN 
UP 183 xLBL4 
184 RELS 
185  REL? 
186 pP 
10?  RCLi 
ı188 DS2I 
189  RCLi 
118 P 
111  6$Bd 
_—_ Ui RE 


up 


118 _RIN 


U+8 


A«xB 


AU:8 


Ausgabe 


- 


113 sLBLS 
114 RCLB 
115  RCL? 
1l6 G6SRd 
117 +R 
119 «LBia 
126 SFZ 
121 »LBic 
122 KEY 
123 Kr 
124 + 
"125 F2? 
126 6708 
127 LSTK 
128 z 
129  LETK 
138 
131 xLBLö 
132 Ri 
133 + 
134 Rt 
135 RTN 
136 LBbLa 
137 EraY 
138 R+ 
139 * 
148 Ri 
141 ev 
142 = 
143 kt 
144 RTN 
145 _ sLBLC 
146 [d 77 
147 l 
148 STül 
149 #LBL6 
158 RCLI 
151 PSE=DJI 
152 RCLi 
153 PRTX=DI3; 
154 ISZI 
155 Fo? 
156  6T08& 
157 SF8 
158 670% 
159 &LBLE& 
168 sPpc 
161 cFe 
162 #LBL$ 
163  RCLI 
164 7A 
165  x#Y? 
166  6T06 


Nach der Flag-Abfrage in Zeile 125 wird dies automatisch zurückgesetzt. Mit den Unterprogramm- 
Aufrufen GSB 1 und GSB 5 in den Programmzeilen 055, 056; 062, 063; 069, 070 werden die 
Sternimpedanzen gemäß den Gin. (2.4.4a bis c) berechnet. 

Die Ergebnisse werden nach Austausch der Primär-Sekundär-Speicher wieder sinngemäß in STO 1 
bis STO 6 abgespeichert. Das logische Programmende ist mit der R/S-Anweisung in Zeile 075 er- 
reicht. Da die beiden Umwandlungsprogramme die Ergebnisse jeweils in dem gegenüber der Start- 
konstellation sekundären Speicherbereich ablegen und mit diesem aktuellen Zustand enden, kann 
eine Rückverwandlung des Systems unmittelbar durch Betätigung der entsprechenden Starttaste 
angeschlossen werden. Damit ist gleichzeitig ein einfacher Test der Umwandlungsgenauigkeit mög- 
lich. 

Das Umwandlungsprogramm benötigt insgesamt 144 Programmzeilen. Zur Anzeige oder Ausgabe 
der Ergebnisse kann ein zusätzliches Hilfsprogramm über die Startadresse Label C aufgerufen 
werden. Mit diesem Programm von Programmzeile 145 bis 166 werden die sechs Speicher STO 1 
bis STO 6 paarweise, durch eine Leerzeile getrennt, ausgedruckt bzw. jeweils im Anschluß an die 
kurz angezeigte Speichernummer angezeigt. 

Falls das Hilfsprogramm generell zur Ausgabe der Ergebnisse herangezogen werden soll, ist es 
zweckmäßig, die beiden R/S-Anweisungen in den Zeilen 041 und 075 durch eine Sprunganweisung 
GTOC zu ersetzen. 


2.4.4 Test- und Anwendungsbeispiele 


a) Eine Sternschältung mit den komplexen Widerständen Z,, Z, und Z, soll in eine äquivalente 
Dreieckschaltung umgewandelt werden. 


Zahlenwerte Z,=(10+j20) 02 Eingaben: 10 STO 1 Re {Z,} 
Z, = (40 - 550) Q 20 STO2 Im{Zı} 
23 = (20 + j10) 02 40 STO3 Re{Z} 


-50 STO4 Im {Z} 


20 STO5 Re{Z;} 
10 STO6 Im {Z3} 


Start [A]: 


Nach dem Programmstop können die errechneten Dreieckimpedanzen aus den nun aktivierten 
Sekundärspeichern STO 1’ bis STO 6’ einzeln oder mit Start automatisch abgerufen 
werden: 


RCL1:  112,002=Re{Zı>} 
RCL2: - 46,0002= 1m {Zı2} 
RCL3: 62,00 02 = Re (Za3} 
RCL4: - 104,00 02 = Im {Z.3} 


RCLS: 23,90 2 = Re (Z31} 
RCL 6: 34,88 N2= Im {Z31} 


Die Rechenzeit für die Umwandlung beträgt rd. 30 Sekunden. Für die Rückverwandlung der 
Dreieckschaltung in eine äquivalente Sternschaltung kann unmittelbar durch Betätigung der 
Starttaste B erfolgen. Danach stehen die errechneten Sternimpedanzen in den Speichern 
STO 1 bis STO 6 zur Verfügung. 
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b) Eine Dreieckschaltung mit den komplexen Widerständen Z}2, Za3, Zaı soll in eine äquiva- 
lente Sternschaltung umgewandelt werden: 


Zahlenwerte Z}2= 100% Eingaben: 100 STO 1 Re {Z42} 
Z33 5500 0 STO2 Im {Zı2} 
Zu = 500 0 STO3 Re{Zs} 


50 STO4 Im {Zus} 


0 STO5 RefZy} 
-50 STO6 Im{Z} 


Start [e}: 

Nach dem Programmstop können die errechneten Stern-Impedanzen aus den nun aktivierten 
Sekundärspeichern STO 1’ bis STO 6’ einzeln oder mit Start automatisch abgerufen 
werden: 


Start [e}: 
zZ, =(0-j50) 0 
Z, = (0 +j50) 0 
Z3 = (25+j0) 2 


Rechenzeit rd. 15 Sekunden 


2.5 Harmonische Analyse 
2.5.1 Numerische Berechnung der Fourier-Koeffizienten 


Das theoretische Rüstzeug zur numerischen Berechnung der numerischen Fourier-Koeffizienten ist 
in Band 1 dieser Reihe dargestellt. Bei der hier angegebenen Programmvariante werden vier Symme- 
trieeigenschaften berücksichtigt und zusätzlich eine Effektivwert- und Klirrfaktorberechnung durch- 
geführt [11]. Bezüglich der Symmetrieeigenschaften werden folgende Unterscheidungen mit Hilfe 
einer in STO A abgespeicherten Kennzahl 0, 1, 2, 3 oder 4 getroffen: 


O keine Symmetrie, Stützpunktintervall 27, n< 16 
ungerade Symmetrie, d.h. f(x) =-f{-x), Stützpunktintervallm, n< 17 
gerade Symmetrie, d.h. f(x) = f{-x), Stützpunktintervallm, n< 17 


{STO A} = j a & 
alternierend, d.h. f(x) = - f(x - m), Stützpunktintervallm, n< 17 


Son — 


um 7 versetzte Symmetrie, d.h. f(x) = f(x - m), 
Stützpunktintervall m, n< 17 


Im Gegensatz zu dem in Band 1 angegebenen Programm werden hier für alternierende Funktionen 
nur die wirklich vorhandenen ungeraden Koeffizienten und für um r-versetzte Funktionen nur die 
wirklich vorhandenen geraden Koeffizienten berechnet und ausgegeben. Die Zahl der Stützwerte 
wird in STO B eingegeben. Es können je nach Symmetriebedingung maximal 16 oder 17 Stütz- 
werte verarbeitet werden. Damit ergibt sich bei dem Stützpunktintervall 27r (0 in STO A) mit maxi- 
mal 16 Stützwerten ein minimaler Stützpunktabstand von 22,5° und bei dem Intervall m (1 bis 4 

in STO A) mit maximal 17 Stützwerten ein minimaler Stützpunktabstand von 11,25°. 

Für das Intervall 2rr werden die numerischen Koeffizienten A, bis maximal Ag (cos-Glieder) und 
B, bis maximal B, (sin-Glieder) berechnet, während für das Intervall m die Koeffizienten A, bis 
maximal A,s und B} bis maximal B,s berechnet werden. 
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2.5.2 Effektivwert und Klirrfaktorberechnung 


Der Effektivwert einer periodischen, nicht-sinusförmigen Funktion u({t) ist gleich der positiven 
Wurzel aus der Summe der Quadrate aus dem Gleichanteil U, und aller Effektivwerte der sinus- 
förmigen Oberschwingungen User: 


Urmsfus ru: „Hu rt (2.5.1) 
Mit: 
Ük 
Ukert = — (2.5.2) 
v2 
Als Unterscheidungsmaß, inwieweit sich eine periodische, nicht-sinusförmige Funktion u(t) von 
der Sinusform unterscheidet, verwendet man den Klirrfaktor K. Unter dem Klirrfaktor K versteht 
man daher das Verhältnis aus dem Effektivwert der Oberschwingungen zu dem Effektivwert der 
Gesamtschwingung. 


2 2 2 
_ Un on + Uger + Ugen +. 
2 2 2 2 
Wut Vor + Uzen + Ugan + --- 


Bei der Fourier-Analyse fallen als Ergebnis die numerischen Fourier-Koeffizienten A, für die 
Amplituden der Cosinusglieder und B, für die Amplituden der Sinusglieder an. Die Amplitude 
der Gesamtschwingung C, mit der Kreisfrequenz kw ergibt sich dann als Wurzel aus der Quadrat- 
summe der Einzelamplituden: 


84 /AR + BL (2.5.4) 


Mit Gl. (2.5.4) läßt sich für den Effektivwert schreiben: 


K (2.5.3) 


em 
Um =) Brit (2.5.5) 
k=2 


K= ——— (2.5.6) 


In dem Programm „Harmonische Analyse 2’ sind für die drei Terme in den Gin. (2.5.5 und 2.5.6) 
folgende Speicher vorgesehen: 


n n 
STO 17: S'ck= D tAR+BR) (2.5.7) 
k=2 =2 
STO 18: C?=A?+B3 (2.5.8) 
STO 19: Co=Ay (2.5.9) 
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2.5.3 Programmbeschreibung „Harmonische Analyse 2” 


Der Eingabeteil des Programms, siehe Tabelle 2.5.1, wird durch Betätigung der Taste E aufgerufen 
und führt an den physikalischen Anfang des Programms. In Programmzeile 002 wird das Unter- 
programm Label d ab Zeile 211 aufgerufen. 

In diesem Unterprogramm wird in Zeile 215 das Flag 1 immer dann gesetzt, wenn der Inhalt von 
Speicher A ungleich Null ist. Vor dem Rücksprung in das rufende Programm wird bei Berücksichti- 
gung der Symmetriebedingung der Wert 1 in das X-Register übernommen. Ist keine Symmetriebe- 
dingung zu berücksichtigen, so wird der Wert 2 übergeben. Anschließend wird die Schrittweite Ax 
berechnet und in Zeile 009 nach STO C abgespeichert. Mit den Anweisungen 010, 011 wird das 
Indexregister null gesetzt. Die Eingabeschleife umfaßt die Zeilen 012 bis 024. 

Innerhalb der Eingabeschleife wird die jeweilige Koordinate des Stützwertes x; ausgedruckt und der 
Index i bei den Eingabeunterbrechungen in Zeile 018 angezeigt. 

Nach dem Programmstop in Zeile 018 ist der Index i im X-Register und der Winkel im Y-Register 
abgespeichert. Mit diesen Werten kann vom Tastenfeld aus gegebenenfalls eine Berechnung der 
Stützwerte erfolgen. Der Inhalt des X-Registers wird nach Betätigung der R/S-Taste als Stützwert 
übernommen. Die Eingabeschleife wird verlassen, wenn die in STO B vorgegebene Anzahl der 
Stützpunkte abgearbeitet wurde. Mit Label C in Zeile 025 beginnt das eigentliche Berechnungs- 
programm. Mit den Anweisungen 027 bis 032 werden die Sekundärspeicher STO 7’, STO 8’ und 
STO 9° für die Berechnung des Effektivwertes und des Klirrfaktors null gesetzt. Die Flag 1 Abfrage 
in Zeile 037 wird bei alternierende — oder um r-versetzte Funktionen übersprungen, da bei diesen 
sowohl A-Koeffizienten wie auch B-Koeffizienten vorhanden sein können. Falls das Flag 1 nicht 
gesetzt war, d.h. die Stützpunkte sind über dem Intervall O, 27 verteilt, wird das Flag 2 gesetzt und 
damit eine Weiche gestellt, um zunächst die A-Koeffizienten und anschließend die B-Koeffizienten 
zu berechnen. 

Ist dagegen das Flag 1 gesetzt, wird in der Zeile 045 durch eine Äquivalenzabfrage festgestellt, ob 
die vorliegende Funktion als ungerade Funktion oder als gerade Funktion behandelt werden soll. 
Ist in STO A eine 1 abgespeichert, d.h. die Funktion ist ungerade, erfolgt ein Sprung nach Label B 
in Zeile 052, so daß nur die B-Koeffizienten berechnet werden. Anderenfalls werden nur die 
A-Koeffizienten berechnet. In beiden Fällen wird das Programm in Zeile 173 mit der Berechnung 
des Effektivwertes und des Klirrfaktors fortgesetzt. 

Mit den Anweisungen 3, GSB 4 in den Zeilen 048 und 049 bzw. 4, GSB 4 in den Zeilen 053 und 
054 werden die in STO E abzuspeichernden Ordnungszahlen der Anfangskoeffizienten ermittelt: 


Bei {STO A} =0,1,2 und4 beginnen die A-Koeffizienten mit Ayo, bei STO A = 3 dagegen 
mit A. 
Die B-Koeffizienten beginnen bei {STO A} =0,1,2 und 3 mit B, und bei STOA=4 mit B.. 


Die Berechnung der Fourier-Koeffizienten beginnt mit dem Sprung nach Label a in Zeile 062. 

Hier werden zunächst das Indexregister und der als Summierspeicher vorgesehene Speicher D null 
gesetzt. Die Summierschleife erstreckt sich von Programmzeilte 067 mit Label b bis zu der Rück- 
sprunganweisung GTO b in Zeile 109. Innerhalb dieser Schleife ist mit den Anweisungen in den 
Zeilen 082 bis 101 die Korrektur des ersten und letzten Summanden für die Analyse symmetrischer 
Funktionen in dem Intervall [O, 7] eingebaut. Die Multiplikation mit einer Sinus- oder Cosinusfunk- 
tion des Abszissenwertes wird über die Flag O-Abfrage in Zeile 073 gesteuert. Jeder berechnete 
Koeffizient wird durch die Anweisungen in den Zeilen 110 bis 115 dahingehend überprüft, ob sein 
Betrag größer als die gesetzte Schranke e = 10”? ist. Falls dies nicht der Fall ist, werden die An- 
weisungen für die Ausgabe bis Label O in Zeile 154 übersprungen und die Ordnungszahl zur Be- 
rechnung des folgenden Koeffizienten gebildet. 
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Tabelle 2.5.1. Anweisungsliste „Harmonische Analyse 2‘ 


Start 1 861 _&LBLE 057 + 113 
082 65Bd 114 
8653 I 115 
884 8 116 
885 L 117 
886 x 116 
667° 6583 119 
888 * 128 
889 STOC 121 
818 8 122 
si STOI 123 
812 «LBic 124 
813  RCLC 125 
814  RCLI 126 
815 'x 127 
8i6 PRIX 128 
#1? LSTX 129 
e18 R-S 138 
819  STOs 131 
e2e 1821 132 
821  RELI 133 
022  RCLB 134 
823 3%» 135 
824 __6Tüc 136 

Start 2 &5 _ xLBLC 137 
826 SPC 138 
827 & 139 
828 FZS 148 
e2s  STür 141 
83a STÜe 142 
831  STÜs 143 
832 P&S 144 
0833 RÜLA 145 
834 3 146 
835 XzY? 147 
056 67Ü1 148 
827 Fi? 149 

—0838  6T08 158 
032 #LELI 151 
048 SF2 152 

—841 6Tün 153 
842 «LBLE 154 
043 RÜLH 155 
844 1 156 
045 x=Y? 157 
846  6TOB 158 

»647 #LBLA 159 
848 3 168 
8402 6554 161 
858 SF& 162 
851  ETO« 163 
0852 #LBLB 164 
853 4 165 
054 6584 166 
855 l 167 
056 RiLE 168 


3 169 x©Y? 
CHS 178  6T0a 
x)Y? 171 F2? 
6T08 172  6TOB 
2 173 25 
3 174 | RCL? 
$Tol 175 | RCLE 
65B3 176 + 
ST#i 177 2 
ESBd 178 & 
* 179 | RCL9 
RCLE 188 + 
X=Y? 181 Tk 
ETO1 182 | PRTISDCgH 
x=0? 183 | RÜL? 
6701 184 | ST+& 
2 185 | RÜLE 
ST 186- = 
sLBL] 187 5 
DSPa 188 PN 
RCLE 189 PRIX=DK 
FRTX=Di 198 RS 
RELB P 191 #LBL3 
DSP3 182  RCLB 
PRTZ=DA,B; 193 Fi? 
STxi 124 6Ti8 
1 [: 95 RTN 
$ 136 aLBid 
sToI 197 1 
RÜLE 138 = 
x=8° 199 RTN 
8701 UP Zz08 «LBLs 
DSZI ze1  RCLA 
1 202  x=Y? 
x=y? 203 8T01 
ETO1 284 8 
D521 —285 6Tüe 
sLöll 286 &LBLI 
RCLD 207 l 
ST+i 208 «LBLO 
SPÜ 289  STOE 
*LBLö 218 RTN 
2 Up 211 #LBLd 
4 212 eFi 
sTul 213  RCLA 
ISZi 214 x#89? 
2 215 SFI 
RÜLR 2l6 R+ 
Kr? 217 1 
1523 218 Fi? 
RCLs 219 RTN 
1 228 Rs 
- 22 2 
6SB3 222 RTN 


65Ba 
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Falls der berechnete Koeffizient größer als die gesetzte Schranke von e = 10° ° ist, wird in der 
Zeile 134 die Ordnungszahl des Koeffizienten ausgegeben. Mit den Anweisungen in den Zeilen 117 


bis 131 wird der Koeffizient = berechnet und die Koeffizienten A;,, B; mit dem Faktor 2 wert- 
richtig gebildet. Die berechneten Fourier-Koeffizienten werden mit drei Kommastellen in Zeile 137 
ausgedruckt und gegeneinander durch eine Leerzeile abgesetzt. Mit den folgenden Anweisungen 
von 138 bis 152 werden die Quadrate der berechneten Fourier-Koeffizienten durch indirekte Adres- 
sierung den Sekundärspeichern STO 7’ bis STO 9’ nach Maßgabe der in STO E abgelegten Ord- 
nungszahl zugeordnet. 

In den Zeilen 155 bis 157 wird die indirekte Adressierung für den Ordnungszahlenspeicher STO E 
vorbereitet. Falls eine alternierende — oder um rr-versetzte Funktion vorliegt, wird die Ordnungs- 
zahl um einen Doppelschritt heraufgesetzt. In Zeile 163 wird die neue Ordnungszahl in das X- 
Register gerufen. Mit Hilfe nachfolgender Unterprogrammaufrufe GSB 3 und GSB d wird die 
obere Grenze für die Ordnungszahl festgelegt und in Zeile 169 mit dem um 1 verminderten Inhalt 
von STO E verglichen. Ist die Grenze noch nicht erreicht, so wird der nächste Berechnungszyklus 
bei Label a in Zeile 062 begonnen. Anderenfalls werden je nach dem Zustand von Flag 2 die 
B-Koeffizienten berechnet oder die abschließende Berechnung wird eingeleitet. Hierzu stehen in 
STO 7’ die Summe aller Quadrate der Oberschwingungsamplituden gemäß GI. (2.5.7), in STO 8° 
das Quadrat der Grundschwingungsamplitude gemäß Gi. (2.5.8) und in STO 9’ das Quadrat des 
Gleichanteils gemäß Gl. (2.5.9) zur Verfügung. Der Effektivwert nach Gl. (2.5.5) wird in Zeile 182 
ausgegeben. Anschließend wird der Klirrfaktor gemäß Gl. (2.5.6) berechnet und in Zeile 189 aus- 
gegeben. Das logische Programmende ist mit Zeile 190 erreicht. 

Über die zweite Startadresse Label C in Zeile 025 kann das Programm ohne neue Dateneingabe 
wiederholt gestartet werden. Dies ist von Vorteil, wenn nur einzelne Stützwerte in den Speichern 
STOO bis STO 16 verändert werden. Vor Ausgabe der Sinus-Koeffizienten B; wird durch zwei 
Leerzeilen mit den Anweisungen in den Zeilen 060, 061 eine Trennmarkierung auf dem Druck- 
streifen erzeugt. 


2.5.4 Test- und Anwendungsbeispiele 


Am Beispiel einer Sägezahnfunktion mit der Amplitude y = 1 soll im folgenden eine Fourier- 
Analyse für die fünf möglichen Funktionsarten durchgeführt werden. 


1. Sägezahn über dem Intervall 27 gemäß Bild 2.5.1: 
Eingaben: 0 STOA 


16 STOB 
Start [e} 


Nach Anzeige O Eingabe y, 


05, IR/S 
1/16, |R/S 
2/16, |R/ 

3/16, |R/S 


Nach Anzeige 1 Eingabe y, 


Nach Anzeige 2 Eingabe Y2 


Nach Anzeige 3 Eingabe y, 


ee) 2||n 
(2) a 


Nach Anzeige 14 Eingabe y,, = 14/16, |R/ 


Bild 2.5.1 Im Intervall 27 periodische 
Nach Anzeige 15 Eingabe y,, = 15/16, [R/S Funktion {STO A}=0 r 
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Ergebnisse: exakte Werte der Fourier-Koeffizienten:" 


Gleichanteil: Au = 0,500 a9 = 0,5 
B}; = -0,314 b} = - ı = - 0,318 
1 
B- = - som. 
5 0,151 b, = 0,159 
1 
B- = - ser es 
3 0,094 bz a 0,106 
sin-Glieder 1 
Bi Kot: B, = -0,063 ba == 7, = 0,080 
1 
B = - =-— =- 
5 0,042 bs Er 0,064 
a Een ge 
Bs = -0,026 be == 57 = 0,063 
1 
B- = - =- — =- 
7 0,012 b, ur 0,045 
Ver = 0,564 Effektivwert (exakter 1/./3 = 0,577) 


K 


0,527 Klirrfaktor 
Rechenzeit rd. 15 Minuten 


2. Sägezahn über dem Intervall m als ungerade Funktion gemäß Bild 2.5.2: 


Eingaben: 1 STOA y 
17 STOB 


-In 


e a IH 
1116, ER 
R/S 


Start [e} 


Nach Anzeige O0 Eingabe y, 


Nach Anzeige 1 Eingabe y, 


Nach Anzeige 2 Eingabe y, = 2/16, Bild 2.5.2 Ungerade Funktion mit 
Nach Anzeige 3 Eingabe y„ = 3/16, |R/S {STO A}=1 
h h lot) =-f(-wt) 


Nach Anzeige 15 Eingabe y,,„ = 15/16, |R/S 
Nach Anzeige 16 Eingabe y,g = 1, R/S 


1) Anmerkung: Die numerischen Fourier-Koeffizienten liefern in der Synthese die Stützpunkte exakt wieder. 
Außerhalb der Stützpunkte treten jedoch Abweichungen zu dem vorgegebenen Kurvenverlauf ein, die um so 
größer sind, je stärker der Anteil der Oberschwingungen höherer Ordnungszahlen ist {verursacht durch Ecken 
und Spitzen). Da bei der numerischen Fourieranalyse unterstellt wird, das nur ein begrenztes Spektrum 
existiert (abhängig von der Zahl der Stützpunkte), ist die Abweichung bei den höheren Harmonischen be- 
sonders groß. 
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Ergebnisse: exakte Werte der Fourier-Koeffizienten: 


B, = 0,835 b, = 2 = 0,637 
B> = -0,314 b> Bi 5 = - 0,318 
2 
Bz2 = 0,206 b; = * 0,212 
2 
B4 = 0,151 ba = Ar = - 0,159 
2 
Bs = 0,117 bs = Er Zi 0,127 
2 
Bs = -0,094 bs = 67 = - 0,106 
2 
B;, = 0,076 b, = 77 = 0,091 
sin-Glieder 2 
r 8 ” —- 0,063 bg ee er —-0,080 
B. sinkwt: Br 
B; = 0,051 bg = 2 - 0,071 
2 
Bio = -0,042 bjo = 107 - 0,064 
2 
B}1ı = 0,033 b}1 = Fr 0,058 
2 
Bja = - 0,026 bia —ze 127 —- 0,053 
2 
Bıa = 0,019 b4a = Ian 0,049 
2 
Bja = -0,012 bja =- rm - 0,045 
2 
Bis = 0,006 bis = u 0,042 
Yerr = 0,550 Effektivwert (exakter Wert 1/\/3 = 0,577) 
K=- 0,579 Klirrfaktor 


Rechenzeit rd. 15 Minuten 


Sägezahn über dem Intervall als gerade 
Funktion gemäß Bild 2.5.3: 


Eingaben: 2 STOA 
17 STOB 


Bild 2.5.3 


Gerade Funktion mit {STO A} = 2 
flwt) = f(- wt) 


Da die gleichen Stützwerte wie beim vorhergehenden Fall 2. gelten, kann das Eingabepro- 
gramm umgangen werden. Die Berechnung der Fourier-Koeffizienten wird dann über Label C 
gestartet. 


Start [c}: 


Ergebnisse: exakte Werte der Fourier-Koeffizienten: 
Gleichantel: Au = 0,5 . = 05 
2 2 
Ar = 0,407 4 =” (2) = —-0,405 
2 2 
Aa = 0,046 a, = -(2) = 0,045 
2 2 
A; = 0,018 5 - -(2) = -0,016 
2 2 
cos-Glieder Ag = -DDID an =” (7) = 0,008 
; ; 2 
A'ensktiti | A, = 9,002 a = -(2) = -0,005 
97 
2 2 
Aıı= 0,005 a = -() = -0,003 
2 2 
Aıs3= 0,004 a3 = -() = -0,002 
2 2 
As” - 0,004 dı5 7 (=) - - 0,002 
Yerr = 0,578 Effektivwert (exakter Wert 1/3 = 0,577) 
K = 0,126 Klirrfaktor 


Rechenzeit rd. 15 Minuten 


Sägezahn über dem Intervall als alternierende Funktion gemäß Bild 2.5.4: 


Eingaben: 3 STOA 
17 STOB 


Gleiche Stützwerte wie im Fall 2. 


Bild 2.5.4 


Alternierende Funktion mit (STO A} = 3 
f(wt) =-flwt-m) 
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Start [ec]: 


Ergebnisse: exakte Werte der Fourier-Koeffizienten: 
2 
Aı = 0,407 a4 = —- (2) = -0,405 
2 2 
Az = -0,046 a3 = -(2) = —-0,045 
2 2 
As = -0,018 4% 7 -(2) = -0,016 
2 2 
i R 2 
Auspskabi]| 1, = -0007 29 = -(2) = -0,005 
9r 
2 2 
Ayı = - 0,005 441 7 (3) = - 0,003 
212 
Ars= -0,004 a1 = -(Z.) = -0,002 
2 \2 
Ass — 0,004 dı5 7 -(2,) = -0,002 
2 
Bı = 0,635 b, = 7 = 0,637 
= -— 
B3 = 0,206 bz = = - 0,212 
B;s = 0,117 b; = 2% 0,127 
5 ‚ 5 57 ; 
B, = 0,07 b, = 2 = 0,09 
sin-Glieder = ‚076 un m 1 
Be sinkt | 5, = Hd Se en, 
9n 
Bıı= 0,033 b4] = eo; 0,058 
11 D 11 117 ‚ 
Bı3= 0,019 b13 = ER 0,049 
13 ‚ 13 137 4 
2 
Bıs = 0,006 b1s = ı* 0,042 
Yerr = 0,564 Effektivwert (exakter Wert 1/\/3 = 0,577) 
K = 0,329 Klirrfaktor 


Rechenzeit rd. 15 Minuten 


5. Sägezahn über dem Intervall m als um r-versetzt symmetrische Funktion gemäß Bild 2.5.5: 


Eingaben: 4 STOA 
17 STOB 
y 
Gleiche Stützwerte wie im Fall 2. 
10 
Bild 2.5.5 058, * —e ut 
0 


Um n-versetzt symmetrische Funktion 


-1, 
{STO A}=4, flot) = flot-m) SE. 


Start [c} 


Ergebnisse: exakte Werte der Fourier-Koeffizienten 
Gleichanteill: An = 0,500 %. = 05 
Bz = 0,314 b, = - ı = -0,318 
1 
Ba = 0,151 b, = a 0,159 
1 
B; = -0,094 b&, = ee 0,106 
sin-Glieder 1 
Bas = -0,063 bs = - > = -0,080 
B. sin kwt: ar 
1 
Bjo = -0,042 bio = nn 0,064 
1 
B}2= -0,026 b42 = nn 0,053 
1 
B}4= -0,012 b44 = er 0,045 
Verr = 0,564 Effektivwert (exakter Wert 1/\/3 = 0,577) 
K = 1,000 Klirrfaktor (keine Grundschwingung) 


Rechenzeit rd. 15 Minuten 


6. Thyristorsteuergeräte mit Phasenanschnitt und Phasenabschnitt der sinusförmigen Stromkurve 
haben den Vorteil, daß sie keine Grundschwingungsblindleistung erzeugen [12]. Jedoch können 
hier die Oberschwingungen höherer Ordnung besonders große Werte gegenüber der Grund- 
schwingungsamplitude annehmen. An einem Beispiel gemäß Bild 2.5.6 soll ein Strom i=1 sin wt 
mit 1=1A durch periodischen Phasenabschnitt bei 1/4 + km und Phasenanschnitt bei 
3/4 + km gesteuert werden. Für 17 Stützwerte soll eine harmonische Analyse durchgeführt 
werden. 


Eingaben: 1 STOA {ungerade Funktion) 
17 STOB (Anzahl der Stützwerte) 
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Start [E] 


Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 
Nach Anzeige 


Nach Anzeige 12 mit Ausdruck 135,00°: 
Nach Anzeige 13 mit Ausdruck 146,25°: 
Nach Anzeige 14 mit Ausdruck 157,50°: 
Nach Anzeige 15 mit Ausdruck 168,75°: 
Nach Anzeige 16 mit Ausdruck 180,00°: 


0 mit Ausdruck 0°: 


1 mit Ausdruck 
2 mit Ausdruck 
3 mit Ausdruck 
4 mit Ausdruck 
5 mit Ausdruck 
6 mit Ausdruck 
7 mit Ausdruck 


8 mit Ausdruck 


9 mit Ausdruck 101,25°: 
10 mit Ausdruck 112,50°: 
Nach Anzeige 11 mit Ausdruck 123,75°: 


11,25°: 
22,50°: 
33,75°: 
45,00°: 
56,25°: 
67,50°: 
78,75°: 
90,00°: 


R/S 


|| 
oa|llo 


Pe] 
m 
42] 


x x x 
< I |< | |I< 
[>] 


(halber Sinuswert) 


x 
N 
< 


21|2/|9 
a 
olıon| Io 


e) 
— 
[72] 


BISIGIGIENGIE 


x\Ix||x | [x 
nn] In] IN 
<||< | |< |I< 
EIEIEIER 


2 
S 


je] 


E] (halber Sinuswert) 


3 
Be 
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22/117 
— I 
on 


2) 
> 
[2] 


Der zweimalige Zeilenvorschub zu Beginn der Rechnung zeigt an, daß nur Sinus-Koeffizienten 


Ergebnisse: 

existieren. 
B} = 0,186 
Bz = 0,14 
B;s = 0,094 
B, = -0,094 
Ba = -0,042 
B}1 = 0,042 
Bı3 = 0,012 
Bis = -0,012 
ler = 0,278 
K= 0,881 


Rechenzeit rd. 8 Minuten 


Effektivwert 
Klirrfaktor 


Grundschwingungsamplitude sin wt 


Oberschwingungsamplituden sin kwt 


In Bild 2.5.6 sind die Schwingungsanteile bis zur Ordnungszahl 7 eingetragen. 

Aufgrund der begrenzten Anzahl Stützwerte, die zu einer Schrittweite von minimal 11,25° führt, 
ergeben sich insbesondere in der Umgebung von Unstetigkeitsstellen Ungenauigkeiten der Approxi- 
mation. In dem dargestellten Beispiel wurde z.B. der Funktionswert an den Sprungstellen bei 45° 
und 135° gleich dem halben oberen Funktionswert gesetzt. Eine stetige Verbindungslinie der Stütz- 
punkte kappt einen Teil der Spitze der gegebenen Funktion ab. Dadurch ergibt sich beim Effektiv- 
wert ein um rd. 8 % zu kleiner Wert gegenüber dem exakten Wert: 


„1 V 2, 
-iz 1-7=1:0,3014 


RER EBEN ERZUR TEN 
0 4 Q8_.y 85 


Bild 2.5.6 Harmonische Analyse einer sinusförmigen Stromkurve mit Phasenanschnitt- und 
Phasenabschnittsteuerung 
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3 Nachrichtentechnik 


3.1 Schwingkreisberechnung 


Ausgehend von dem allgemeinen Fall der Parallelschaltung einer Induktivität mit ohmschen Wider- 
stand in Reihe und einer Kapazität mit ohmschen Widerstand in Reihe, Bild 3.1.1, wird ein Pro- 
gramm zur Berechnung des komplexen Gesamtwiderstandes, der Stromverteilung und der Resonanz- 
frequenz angegeben [13]. 


Bild 3.1.1 Parallel-Schwingkreis 


Für den Gesamtstrom I gilt: 


I=UY=-U(G+jB) (3.1.1) 
mit dem Wirkleitwert G und dem Blindleitwert B: 
R R 
G=4J+7 (3.1.2) 
zii zZ 
Sur 
L 
B-Ö- (3.1.3) 
2, z, 


und den Beträgen der Zweigimpedanzen: 


Zı= YR? + (wL)? (3.1.4) 
2,= IR (4) (3.1.5) 
2, ) 2 wC . 


Der Resonanzfall liegt vor, wenn der Schwingkreis nur Wirkstrom aus dem Netz bezieht, d.h. wenn 
der Blindleitwert B=0 wird. Damit folgt für die Resonanzfrequenz f,: 


1 /L-cR? 
#, (3.1.6) 


2n YLC ) L-CR3 
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3.1.1 Programmstruktur „Schwingkreisberechnung” 


Im ersten Teil des Programms mit der Startadresse Label A werden die Gin. (3.1.1 bis 3.1.5) ausge- 
wertet. Für die Berechnung müssen folgende Eingabewerte bereitgestellt sein: 


Indirekte Adresse {1} 


Widerstand R, in Q > STOA 20 
Induktivität LinH > STOB 21 
Widerstand R, in 0 > STOC 22 
Kapazität CinF >STOD 23 
Frequenz finHz> STOE 24 
Spannung Uin V>STOO 0 


Die Schrittweite der über STO | gewählten Variablen wird in STO 9 abgespeichert. 
Bis zum Programmende in Zeile 077 werden folgende Werte in Polarkoordinaten ausgegeben: 


Impedanz Z inS®2 (Betrag und Winkel) 


Gesamtstrom | inA (Betrag und Winkel) 
Teilstrom I} inA (Betrag und Winkel) 
Teilstrom la inA (Betrag und Winkel) 


Falls einer der Eingabewerte variiert werden soll, z.B. Z, 1 =g(f), so kann dies durch die einmalige 
Eingabe einer Schrittweite in das X-Register nach dem Programmstop und anschließender Betätigung 
der R/S-Taste bewirkt werden. Hierzu ist vorher das Indexregister entsprechend der absoluten 
Speicheradresse der Variablen (z.B. 24 STO I für variable Frequenz f) zu setzen. Über die Start- 
adresse Label B wird die Resonanzfrequenz nach GI. (3.1.6) bestimmt. 


3.1.2 Programmbeschreibung „Schwingkreisberechnung” 


Das Programm, siehe Tabelle 3.1.1, umfaßt 114 Anweisungszeilen. Der Programmteil mit Start- 
adresse Label A endet mit einer R/S-Anweisung in Zeile 077. Die folgenden Anweisungen bis zu 
Label B in Zeile 090 dienen der zyklischen Berechnung von Ortskurvenpunkten für die mit STO | 
indirekt adressierte Variable. Bei {STO I} = 20 wird z.B. der Widerstand R, als Variable mit der 
in STO 9 vorgegebenen Schrittweite verändert. Die Auswertung der Gin. (3.1.2 bis 3.1.5) erfolgt 
mit den Anweisungen 002 bis 044. Mit den Anweisungen in den Zeilen 047 bis 050 wird die Impe- 
danz Z mit Betrag und Phasenwinkel ausgegeben. Die Impedanz Z; ist in den Speichern STO 5, 
STO 6 und die Impedanz Z; ist in den Speichern STO 7, STO 8 mit Betrag und Phasenwinkel 
abgespeichert. Die Anweisungen 056 bis 059 liefern den Gesamtstrom 1 nach Betrag und Phase. 
Anschließend werden noch die Teilströme |, und I» ausgegeben. Vor dem Programmende ist in 
Zeile 076 noch eine Clear X Anweisung eingefügt. Falls vor Betätigung der R/S-Taste ein Wert un- 
gleich Null eingegeben wurde, so wird dieser Wert als neue Schrittweite für die zyklische Berech- 
nung nach STO 9 abgespeichert. Mit der Pause-Anweisung in Zeile 085 wird die gewählte Adresse 
der zyklischen Variablen angezeigt. Der aktuelle Wert dieser Variablen wird in Zeile 087 ausgegeben 
und die Berechnung bei Label A am Programmanfang fortgesetzt. 

Mit dem Programmteil Label B von Zeile 090 bis Zeile 114 wird die Resonanzfrequenz gemäß 

G!. (3.1.6) bestimmt. 


53 


Tabelle 3.1.1 Anweisungsliste „Schwingkreisberechnung‘ 


ZT 881 «LBLA 058 CHS 
082 RCLB 859 PRISDYP, 
0883  STOI a68 SPC 
0884 RÜLE 861 RCLA 
885 2 0862  RCLS 
886 x 863 = 
887 Pi a4 FRTISDI 
088 x 865 RÜLE 
869 STK Be6  CHS 
818  RCLD 867 PRIKDYpn, 
LIRN x es SFC 
2 18 8693 RCL6 
813  STO2 878 RELF 
814 RCLC e71 * 
815 >P 72 PRINDIL 
als  STr 873  RCLE 
87 .Y 074 PRIXDYn 
eis STO8 e75 SPC 
es xy 075 LK 
826 x 877 __R#S 
21 STO4 = 678 #0” 
Bez RCLI 873  ETO& 
023  RÜLR [Lee RCL® 
624 Pr 881 #LELB 
ass $T0s 882  STO3 
BEE Ner 883 STri 
827 STÜE 864 RELI 
828 KEY 885 PSE > Adresse 
822 x2 @B6 RÜLi 
es STOS 88? _ FRTS=D[i]=Parameter 
e31 RüL2 888  SFÜ 
032  RCL4 a89 __ETÜH 
633 = fr 098 »LBLE 
634 RÜLI 651 RCLB 
835  RÜL3 892  RCLO 
836 > 893  RÜLA 
837 = 894 ve 
038 RCıa 035 x 
0839 RÜLS 896 = 
848 5 897  RCLB 
841 RCLE 08986  RCLD 
042 RCL4 0899 RELE 
043 * 188 se 
844 + 181 x 
845 F 162 - 
646 18 183  RCLB 
847 PRIS=DZ 184 R 
848 ey 105  RÜLD 
849 ins 186 & 
658 FRIX=DY, 107 + 
51  xer 188 VE 
es2 SPC 189 Pi 
853  RCLA 110 + 
054 Kar 111 2 
055 * 112 = 
0856  PRTX=DI 113 PRTX=Dfr 
057 xer 114 RS 
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3.1.3 Testbeispiel 


Eingabe: Widerstand R, = 102: 10 STOA 
Induktivitätt L= 0,1H: 0,1 STOB 
Widerstand Rz, =202: 20 STOC 
Kapazität C=-10uF: 10 EEX-6 STOD 
Spannung U=100 V: 100 STOO 


Start [8] 


Ergebnis: f, = 161,62 Hz 
Rechenzeit rd. 3 Sekunden 
Eingabe: f, = 161,62 Hz in STO E 


Start [A}: 

Ergebnisse: Z, = 340 0 gie! 
1,2029 Ad 

I1r= 0,98 A e1°438 

I2r= 1,00 A e178.52 


Rechenzeit rd. 15 Sekunden 


Eingabe: 10 STOE {Frequenz f in Hz) 
24 STOI {Adresse des Parameter-Speichers f) 
Start [A]: 


Ergebnise: Z = 11,862 ei 31,78’ e 

I = 843A 9179 

In = 847A e 192.14 

I2 = 0,06A ei®9.28 
Eingabe: 40 R/S (Dieser Wert wird als Schrittweite nach STO 9 übernommen} 
Ergebnisse: Anzeige 24 als Adresse des Parameter-Speichers (Frequenz f} mit folgendem 

Ausdruck: 

f=50Hz z 
236450 ce 
= 2,74 A e)62.97 Im{l} 


= 3,03A e"172. 
= 0,31 A e'96. 


In Bild 3.1.2 ist die Ortskurve für den 
Gesamtstrom | aufgetragen. 


Bild 3.12 Ortskurve für den Gesamtstrom I 


3.2 Wellen- und Betriebsparameter passiver Filterschaltungen 
3.2.1 Allgemeines Reaktanz-Halbglied 


Eine elektrische Filterschaltung (Siebschaltung) ist ein Vierpol, der die Aufgabe hat, von einem am 
Eingang ankommenden Frequenzgemisch nur gewisse Frequenzbänder möglichst ungeschwächt 

zum Ausgang durchzulassen, während die übrigen Frequenzen möglichst stark geschwächt werden 
sollen. Die Frequenzbereiche, in denen das Filter durchlässig sein soll, heißen Durchlaßbereiche (DB), 
die anderen Sperrbereiche (SB). An der Grenze zwischen Durchlaß- und Sperrbereichen liegen die 
Grenzfrequenzen w,. Die vier wichtigsten Filterarten sind: 


Filterart Durchlaßbereich Sperrbereich 


Tiefpaß 
Hochpaß 


Bandpaß 


Bandsperre 


Die Filterwirkung einer Siebschaltung hängt von den Beschaltungswiderständen am Ein- und Aus- 
gang ab. Deshalb sind zur Beurteilung der Filterwirkung die Wellenparameter für reflexionsfreie 
Beschaltung und die Betriebsparameter für nicht reflexionsfreie beliebige Beschaltung von Interesse 
[14, 15]. 

Als allgemeine Berechnungsgrundlage werden zur Beschreibung des verlustlosen Grundhalbgliedes 
gemäß Bild 3.2.1 folgende Normierungskennwerte eingeführt: 


Ro = VZa' Zu (3.2.1) 
Za 


FE 
NT z 


(3.2.2) 


Ro ist der Kennwiderstand und £2 ist die normierte Frequenz (Kennfrequenz). 


= 
A 


Bild 3.2.1 Halbglied mit T- und r-Seite 


a) Tiefpaß b) Hochpaß c) Bandpaß d) Bandsperre 


Bild 3.2.2 Filter-Halbglieder 
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Für die Grundfilterarten Tiefpaß, Hochpaß, Bandpaß und Bandsperre gemäß Bild 3.2.2a bis d gelten 


die Normierungskennwerte: 


&o für Tiefpaß 
wg i 
on für Hochpaß (Rrr EUGEN. KR ) 
Arp 
Q2= P 


für Bandpaß 


3.2.2 Berechnung der Filterelemente 


für Bandsperre (Res =- 3= ) 
BP 


(3.2.3) 


(3.2.4) 


Mit den Reaktanzelementen L,, C, und L,, C, lassen sich die vier Grund-Filterarten durch Kom- 


vination von jeweils zwei Halblgiedern als r- oder T-Schaltung aufbauen {Bild 3.2.2). 


Sind der Kennwiderstand R, und die Grenzfrequenzen gegeben, so sind die einzelnen Reaktanzen 
der Halbglieder durch folgende Gleichungen bestimmt: 


i Ro 
Tiefpaß: L,=-5; C5= 
3 
1 
Hochpaß: C,= Lu = 
ochpa u Fer b 
Ro 
BE T 
Bandpaß: 
i Be = 
BEN og b” 
a2 @gı _ 
L, [6] gı Wg2 C,= 
Bandsperre: 5 
0 
Lb Ga Ga = 


3.2.3 Berechnung der Wellenparameter 


13.2.5a, b) 


(3.2.6, b) 


{3.2.7a, b) 


(3.2.8a, b) 


(3.2.9a, b) 


(3.2.10a, b) 


Für den Wellenwiderstand der T-Schaltung mit den normierten Parametern 2 und Ry gilt allgemein: 


Zwr =Rovi-0°? 


(3.2.11) 
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Der reelle Lösungsbereich gilt für den Durchlaßbereich (DB) und der imaginäre Lösungsbereich gilt 
für den Sperrbereich (SB): 


| RRVN1-02? fürlaısı 


m: \irova?-1 frıaı>ı wen 
Für die n-Schaltung gilt entsprechend: 
| sr für |Q] 1 
Zune Y' gi (3.2.13) 
-j——— für >ı 
2? -1 


Das Verhältnis der Eingangsspannung U} einer Filterschaltung zu der Ausgangsspannung U, bei 
Abschluß mit dem Wellenwiderstand nennt man die Spannungs-Wellenübertragungsfunktion F,.: 


U : 

Eu® m. = en + Id (3.2.14) 
=2 

Free "W=1-202+j20 11-2 (3.2.15) 


Die Wellendämpfung in Neper ergibt sich aus dem natürlichen Logarithmus des Betraganteils der 
Spannungs-Wellenübertragungsfunktion. 


0 für |Q] < 1 
ay = (3.2.16) 


In|FyI=21n{YR?-1+|2) für IQ] >1 


Für das Wellenphasenmaß b,, gilt: 


22v1-02 
| arctan — für [2] s 1 
by =! 1-20 (3.2.17) 
7 für [2] >1 
180° 
etz ba (3.2.18) 


3.2.4 Berechnung für Betriebsparameter 


Als Betriebsparameter sind die Eingangswiderstände Z.r und Z.„ und die Betriebsübertragungs- 
funktion Fg von besonderer Bedeutung. Die Betriebsübertragungsfunktion Fg setzt die vom 
Generator maximal abgebbare Wirkleistung zur Wirkleistung am Abschlußwiderstand zueinander 
ins Verhältnis: 


Prnax Uo / 22 = e’B*ibs 
5. 


Fs*® P, "zu, Z, (3.2.19) 


Die Betriebsübertragungsfunktion stimmt für den Fall der beidseitigen Anpassung Z)= Z, = Z,, 
mit der Spannungs-Wellenübertragungsfunktion überein. 
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Bild 3.2.3 Allgemeines Vierpol-Übertragungsglied 


q 
(Bild 3.2.3), so folgt für die Eingangswiderstände Zer, Zer 


Fan 1-20 +j2Nal1-0°) af 
0 aı-2)+j2n == 


R 
Setzt man beim T-Grundglied Z}=Za =R = — und beim r-Grundglied Z}=Z3=R=sRo 


s(1-20°)+j2Q 
Zen = Ro —— 3.2.21 
-1L-202+j20511-N2 


Die Eingangswiderstände von T- und r-Gliedern haben reziproke Frequenzgänge, Bild 3.2.4. Für die 
Betriebsübertragungsfunktion der T-Schaltung gilt: 


For ze ®T=1-202+j0 au-ar +2] (3.2.22) 
Mit der Betriebsdämpfung in Neper: 

Br = tin [1-20 + 02[an-a+2] (3.2.23) 
und dem Betriebsphasenmaß: 


atı-2?)+z 
bar = arctan | ——————— 


1202 (3.2.24) 

a (3.2.25) 
Entsprechend gilt für die Betriebsübertragungsfunktion der 7-Schaltung: 

For=e" "1-202Hjalsıı-nd+4 | (3.2.26) 


Die Betriebsübertragungsfunktion für das r-Grundglied stimmt mit derjenigen für das T-Grundglied 


Ro 
überein, wenn anstelle der Fehlanpassungsvariablen q = Bi die Fehlanpassungsvariable s = 
gesetzt wird. 


Ro 


1 
a RR RR Bestg 1a, 
dem 1 
"0 = 213; = 


a) T-Schaltung b) n-Schaltung 


Bild 3.2.4 T- und r-Grundglieder zur Berechnung der Betriebsparameter 
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3.2.5 Programmstruktur ‚Filterberechnung” 


Das Programm enthält vier Startadressen für die Berechnung der Filterarten: Tiefpaß (Label A), 
Hochpaß (Label B), Bandpaß (Label C) und Bandsperre (Label D) sowie zwei Startadressen für die 
Auswahl der Schaltungsart in Form einer T-Schaltung (Label E) oder r-Schaltung (Label e). 


Vor dem Programmstart sind folgende Eingabedaten bereitzustellen: 


Kennwiderstand Ro in 2 > STOA 
(untere) Grenzfrequenz fg} inHz > STOB 
(obere) Grenzfreauenz f,, inHz > STOC 


Startfrequenz f in Hz > STOD 
Anpassungsfaktor s bzw. q > STOE 
Schrittweite uer Kennfrequenz AR 
für die zyklische Berechnung > STOI 
Falls der Inhalt von STO E gleich null ist, wird START, 


die Berechnung der Eingangsimpedanz und der 


Betriebsparameter übersprungen. Der Rechnungs- / Berechnung und Ausgabe 


Filterelemente 


ablauf ist in Bild 3.2.5 anhand eines Flußdia- 
gramms aufgezeigt. In der äußeren Zyklusschleife 


wird die Kennfrequenz als Laufvariable nach MaßR- 

; : ’ 2 Kennfrequenz 
gabe der in STO E abgelegten Schrittweite verän- Unterscheidung 
dert. Bei der Schrittweite null in STO E wird nur DB — sB 
ein Lösungsvektor bestimmt. 


Betriebsparameter 
23, Ye, dB, Y8 


Bild 3.2.5 Flußdiagramm „Filterberechnung”’ 
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Wellenparameter 
Zu, Ps Aw. Pu 


3.2.6 Programmbeschreibung „Filterberechnung” 


Das Programm umfaßt 223 Programmazeilen, siehe Tabelle 3.2.1. Es wird nur der primäre Speicher- 
bereich von STO O bis STO 7 benutzt. Die Tiefpaßberechnung beginnt in Zeile 001 mit Label A. 
Als Unterscheidungsmerkmal gegenüber dem Hochpaß wird das Flag O gesetzt. Mit dem Unterpro- 
grammaufruf GSB O in Zeile 008 wird die Kreisfrequenz w, gebildet und in STO O abgespeichert. 
Mit den folgenden Anweisungen werden die Gin. (3.2.5 und 3.2.6) ausgewertet und in Zeile 018 
über den Unterprogrammaufruf GSB 6 die Ausgabe der Werte für die Induktivität und Kapazität 
eines Halbgliedes veranlaßt. Ab Label 2 in Zeile 054 beginnt der Berechnungszyklus mit der Kenn- 
frequenz f2 als Parameter. Falls der Wellenwiderstand imaginär ist (Sperrbereich), wird vor dem 
Ausdruck des Wertes eine Zeile mit dem Merkmal 111... erzeugt. In Programmzeile 126 werden 
die Wellen-Übertragungsparameter a,,, y%, ausgegeben. Die Auswertung der Gin. (3.2.20 und 3.2.21) 
endet mit der Ausgabe des Eingangswiderstandes Z, in Zeile 156 nach Betrag und Phasenwinkel. 
Anschließend werden die Betriebs-Übertragungsparameter ag, pg berechnet und in Zeile 168 aus- 
gegeben. Das Unterprogramm Label 3 in Zeile 182 bis 211 wird zur Berechnung der Bandpaß-, 
Bandsperre-Elemente mit den Anweisungen in den Zeilen 025 und 033 aufgerufen. 


3.2.7 Test- und Anwendungsbeispiele 


Mit den Filterkenndaten Kennwiderstand R, = 600.2, Grenzfrequenz f, = I1kHz für Hoch- oder 
Tiefpaßfilter und Ro = 600.9, fg} = 1kHz, fa = 1,2 kHz für Bandpaß- und Bandsperrefilter sollen 
die Wellen- und Betriebsparameter in Abhängigkeit von der Frequenz ermittelt werden. 

Um einen umfassenden Überblick über die charakteristischen Unterschiede der einzelnen Filter- 
parameter zu gewinnen, werden die numerischen Ergebnisse in Diagrammen dargestellt. 


a) Tiefpaß 


Eingaben: 600 STOA (Kennwiderstand Ro in Ohm) 
1000 STOB (Grenzfrequenz fg in Hz) 
0 STOD_ (Startfrequenz f in Hz) 
1 STOE (Anpassungsfaktor q bzw. s) 
2 STOI (Schrittweite AR für den zyklischen Rechnungsablauf) 


Start [e} Programmvorbereitung für T-Schaltung (Flag 1 wird gelöscht) 


Start [A]: 


Ergebnisse: L,= 0,095 H 


C, = 0,265 uF Halbglied 
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Tabelle 3.2.1 Anweisungsliste „Filterberechnung‘ 


Tefpen BGBl nLELA 57 
Gr SFü 858 
Be __ ETUI 853 

268 
605 Fo 861 
u6s "LELI 62 
80° RCLE 863 
gas ESEU 864 
6a STORE 865 
ale RÜLH 866 
e1ı x a67 
PIE 858 
B12  STur 869 
öl4  RüLa 078 
215 RiLe ar 
FIT # 72 
Br 570 fe er3 
gie GB | 9 a74 
eis  Rcın (Ca e75 
B26  RCLB 875 
e21 F a7? 
BEE__ETOE 78 
62 878 
924 8A 
025 6583 681 
Bes Rlıs B82 
827 Pius 083 
BER RÜLS 654 
G22  RÜL? @35 
a3e__6TOI 86 

o7 
BR 88 
623 6882 ass 
624 RÜLT 896 
@35  RÜLS 831 
a3E  RÜLE as2 
837° sLEL) [La 833 
638 PRSTD(Ca FO 
3 Füc ji | 885 
6 ce I | 886 
641 iv 897 
B4E LSTE a98 
643 RÜLE 893 
644 F 188 
645 x 181 
a4 RÜLE 182 
47 R 183 
a4 = 184 
042 »LELF ies 
0658 Fa? 186 
851 6702 167 
2 186 
03 6 189 
054 HLEL2 i16 
05 FSEDN 111 
056  STO 112 
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FRTS 113 
See 114 
cFzZ 1:5 

sTos 1ls 

ST+E jr 

z: 116 

STüsr 119 
cHs 128 

sTtor 121 

1 122 
* 123 

5To2z 124 

STas 125 

5T+7 126 

>80? 1er 
SPC 128 

x>07 129 

6701 138 
EEN 131 

1 132 

& 133 
Enit 134 
F 135 

* 136 
FRTE. =D 111... 137 
Ry 138 
CH> 138 
SsF2 148 
"LELI 141 
IN 142 

TuS 143 
Fi 144 
1A 145 

ECLA 146 

x 147 

RÜLZ 146 

X78? 14% 

eTcı 158 

Rr 151 
Fi: 152 
CHS 153 

Rt 154 

#LELI 155 

Ry 156 

FRIÄ=DZ, 157 
SPÜ 158 
Fe? 159 

6Tol 168 

RCL3 161 

RCL6 162 

x 163 
RCL?Y 164 
3 165 
Rt 166 

1 167 
6T‘4 168 


*LBL1 169 #LBL1 
SF? 178 RCıI 
1 171 x=8° 
ö 172 R+S 
e 173 ST 
RÜL3 174 RELI 
RCLI 175 6702 
ABS UP 176 «LBL@ 
+ 177 Pi 
*LBL4 178 x 
LN 173 2 
2 188 x 
x 181 RTN 
esBbe>|MliP 182 #LBL3 
RÜLE w 163  RELC 
x=8° 184 RCLB 
6701 185 - 
RCLE i6eE  6SB8 
STx5 167 STC2 
RCL? 1868  RÜLB 
RCLE 189 6588 
$Tx5 198 Rn 
x 181  RELE 
’F 192 65Be 
RCL5 193 F 
RÜL? 134 5703 
»P 185  RELR 
Äey 196 = 
CHS 187  $STOS 
Ri 135  RÜLR 
= 19% RÜLZ 
R+ 208 E 
+ 201  5TO4 
Rt 202 RÜLZ 
Fif 263 RÜLH 
&T05 204 x 
1-X 285 1-8 
x7 20€ Tor 
CHS 287 RCL3 
er 208  RÜLR 
*LBL5 28% N 
RELR 216 STE 
x UP z1 RIN 
esbE > 212 #LBL6 
RcLE 213 PRIX 
RCLE 214 xer 
x 215 TRTX 
RCLE 216 SPC 
17A 217 RTN 
t 
Rei 219 er 
x 228 R7S 
RCL> 
FR P2 [43 w 
N 223 RS 
GEBE > 


1. Zyklus 2. Zyklus 3. Zyklus 4. Zyklus 
0.000 2.000 4.000 6.000 


1111111111. 1111111111. 1111111111. 
1039.230 2323.790 3549.648 


2.634 4.127 4.956 
180.000° 180.000° 180.000° 


Yıy A 

zZ, ng 600.000 1033.887 2322.732 3549.315 
Pe 0.000 -270.512 -270.014 -270.002° 
ag in Np 0.000 2.087 4.159 5.375 
v8 0.000° -150.255 -118.968° -109.190° 


Rechenzeit rd. 20 Sekunden je Zyklus 


n Anmerkung: Der nach 1111... folgende Wert ist imaginär. 


Startet man das Programm mit den gleichen Eingabedaten über den Vorbereitungsschritt [tel 
so werden die Tiefpaß-Filterparameter der n-Schaltung berechnet. 


Start [et Programmvorbereitung für r-Schaltung (Flag 1 wird gelöscht) 


Start [A]: 


Ergebnisse: L,=0,095H | 
C,= 0,265 ur | 


1. Zyklus 2. Zyklus 3. Zyklus 4. Zyklus 
0000 | 2.000 | 4.000 6.000 


1111111111. 1111111111. 1111111111. 
Zwin 600.000 -346.410 -154.919 -101.419 
ay in Np 0.000 2.634 4.127 4.956 
mw 0.000° 180.000° |. 180.000° 180.000° 
Zz.inQ 600.000 348.200 154.990 101.428 
Ye 0.000° 270.512° 270.014° 270.002° 
ag inNp 0.000 2.087 4.159 5.375 
v8 0.000° -150.255° -118.968° -109.190° 


Rechenzeit rd. 20 Sekunden je Zyklus 


Halbglied 


Mit den berechneten Kenndaten für ein Halbglied gemäß Bild 3.2.2a lassen sich folgende Filter- 
schaltungen realisieren: 


Tiefpaß (T) Tiefpaß (m) 
lo lo la 


+ 2. a) T-Schaltung + Ey ' en b) #-Schaltung 


Bild 3.2.6 Tiefpaß in T- und #-Schaltung für Ro = 6000, f,= 1kHz 


In den folgenden Bildern 3.2.7 bis 3.2.12 sind die Kennlinien von Wellenwiderstand, Dämpfung, 
Phasengang und Eingangswiderstand für die beiden Tiefpaß-Grundschaltungen in Abhängigkeit von 


der Frequenz dargestellt. 


1000 
Q 


0 0,5 1,0 15Kkhz 20 


Bild 3.2.7 Wellenwiderstand eines Tiefpasses 
in T-Schaltung (Ro = 600.02, f, = 1 kHz) 


2,5 
Np 
2,0 
a 
m 
10 —— 
IN: — 
D 0 1,5 kHz 2,0 
— 
Pm—h 
0 0,5 1,0 1,5 20 
— N 


Bild 3.2.9 Wellen- und Betriebsdämpfung 
eines Tiefpasses (Ro =600. 0, f, = 1 kHz, 
R,=R,=Ro) 
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0 05 1,0 1,5 kHz 2,0 
——l 
Bild 3.2.8 Wellenwiderstand eines Tiefpasses 
in m-Schaltung (Ro = 600.22, f, = 1 kHz) 


0 05 1,0 15 kHz 2,0 

— nn 

0 0,5 1,0 1,5 2,0 
—N 


Bild 3.2.10 Phasengang der Wellen- und 
Betriebsübertragungsfunktion eines Tief- 
passes (Ro =600. 02, f, = 1 kHz, 
R,=R,=Ro) 


200 


mm fi) 
400 NR 70 


— Re Ze 


0 200 


0 400 600 1200 N 1600 
—e Relz,} 
Bild 3.2.12 Eingangswiderstand eines Tief- 
passes in n-Schaltung (Ro = 600.0, 
fg =1KHz, RR=R,=Ro) 


Bild 3.2.11 Eingangswiderstand eines Tief- 
passes in T-Schaltung (Ro = 600.22, 
f,=1kHz, R=R,=Ro) 


b} Hochpaß 
Eingaben: 600 STOA (Kennwiderstand Ro in Ohm) 
1000 STOB (Grenzfrequenz f, in Hz) 
500 STOD (Startfrequenz f in Hz) 
1 STOE (Anpassungsfaktor q bzw. s) 
05 STOI (Schrittweite AS für den zyklischen Rechnungsablauf) 


Start [e]: Programmvorbereitung für T-Schaltung (Flag 1 wird gelöscht) 


Start [B]: 


Ergebnisse: L,=0,095 H 
C,= 0,265 uF 


1. Zyklus 2. Zyklus 3. Zyklus 4. Zyklus 


1111111111. 1111111111. 1111111111. 

1039.230 670.820 0.000 519.615 

äy in Np 2.634 1.925 0.000 | 
Yw 180.000° 180.000° 180.000° R 
Zz.inQ 1033.887 667.656 268.328 483.735 
Pe 270.512° 272.424° 296.565° 7.125° 
ag in Np 2.087 1.256 0.347 0.008 
B 150.255° 173.884° -135.000° | -60.255° 


Rechenzeit rd. 20 Sekunden je Zyklus 


Halbglied 


Für die Berechnung der r-Schaltung wird die Startfrequenz in STO D auf 500,001 Hz verändert. 
Dadurch wird eine Error-Unterbrechung bei 2 =- 1 entsprechend f = 1kHz wegen Zu, = ® 
vermieden. 


Start Flle} Programmvorbereitung für m-Schaltung (Flag 1 wird gesetzt) 
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Start [8]: 


Ergebnisse: L,= 0,095 H 
C,= 0,265 uF 


1. Zyklus 2. Zyklus 


1111111111. 1111111111. 
-346.411 -536.659 
2.634 1.925 0.000 0.000 
180.000° 180.000° -179.676° | -59.999° 
Zz.inQ 348.201 539.202 1341.654 744.204 
Pe -270.512° -272.424° 63.434° -7.125° 
ag inNp 2.087 1.258 0.347 0.008 
YB 150.255°° 173.885° -134.999° | -60.255° 


" Von dem ausgegebenen Wert (Hauptwert) sind 360° abzuziehen. 


Halbglied 


3. Zyklus 4. Zyklus 


Rechenzeit rd. 20 Sekunden je Zyklus 


Mit den berechneten Kenndaten für ein Halbglied gemäß Bild 3.2.2b lassen sich folgende Filter- 
schaltungen realisieren (Bild 3.2.13): 


Co Co 
-—|| > 
L 
F a) T-Schaltung 


Bild 3.2.13 Hochpaß in T- und r-Schaltung für Ro = 600%, fa= 1KHz 


b) n-Schaltung 


In den Bildern 3.2.14 bis 3.2.19 sind die berechneten Größen für den Hochpaß in Abhängigkeit 
von der Frequenz dargestellt. 


1,0 


kN 
zZ 
Wo 
Rı 
0,50 
0,25 
0 , 
0 0,5 10 15 kHz 2,0 0 05 10 1,5 kHz 2,0 
—f —f 
Bild 3.2.14 Wellenwiderstand eines Hoch- Bild 3.2.15 Wellenwiderstand eines Hoch- 
passes in T-Schaltung (Ro = 600.0, passes in n-Schaltung (Ro = 600 0, 
f,= 1 kHz) fg = 1 kHz) 
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1,0 1,5 kHz 20 


Bild 3.2.16 Wellen- und Betriebsdämpfung 
eines Hochpasses (Ro = 600 02, f, = 1 kHz) 


-750 
0 


250 500 
—e Re l2e} 


750 2 1000 


Bild 3.2.18 Eingangswiderstand eines 
Hochpasses in T-Schaltung (Ro = 600 Q, 
a = 1 kHz) 


c) Bandpaß 


—N 


Bild 3.2.17 Phasengang der Wellen- und 
Betriebsübertragungsfunktion eines Hoch- 
passes (Ro =600. 02, f„= 1 kHz) 


1600 
2 


1200 
Im{Ze} 


-400 


0 400 00 
—— Ref2e} 
Bild 3.2.19 Eingangswiderstand eines 
Hochpasses in n-Schaltung (Ro = 600 0, 
f, = 1 kHz) 


1200 2 1600 


Eingaben: 


Start BG: 


600 
1000 
1200 

800 

1 
2 


STOA 
STOB 
sSTOC 
sSTOD 
STOE 
STOI 


(Kennwiderstand Ro in Ohm) 

(Untere Grenzfrequenz f,} in Hz) 

(Obere Grenzfrequenz fy> in Hz) 

(Startfrequenz f in Hz) 

(Anpassungsfaktor q bzw. s) 

(Schrittweite AQ für den zyklischen Rechnungsablauf) 


Programmvorbereitung für T-Schaltung (Flag 1 wird gelöscht) 
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Start [c} 


Ergebnisse: L, = 0,477H 
C, = 44,209 nF ö 
L,= 0016H Halbglied 
C, = 1,326 uf 


0] 1. Zyklus 2. Zyklus 3. Zyklus 4. Zyklus 


1111111111. 1111111171. 1111111111. 

Zu [9] 2012.461 670.820 1374.773 
aw [Np] 3.850 1.925 0.000 3.134 
Pw 180.000° 180.000° 60.000° 180.000° 
2. [9] 20 10.936 667.656 483.735 1371.229 
Pe 270.826° 272.424° -7.125° -270.159° 
as [Np] 3.759 1.258 0.008 2.751 

5 123.227° 173.884° 60.255° -137.265° 


Rechenzeit rd. 20 Sekunden je Zyklus 


Start [el Programmvorbereitung für r-Schaltung (Flag 1 wird gesetzt) 


Start [c} 

Ergebnisse: L, = 0,477H 
C, = 44,209 nF : 
L,= 0,016H Halbglied 
C.= 1,326 ur 


11111111311. 1111111111. 
-178.885 -536.656 692.820 


1.925 0.000 
180.000° 180.000° 60.000° 
179.021 539.200 744.208 
-270.028° -272.424° 7.125° 
3.759 1.258 0.008 2.751 
123.227° 173.884° 60.255° -137.265° 


Rechenzeit rd. 20 Sekunden je Zyklus 


Mit den berechneten Kenndaten für ein Halbglied gemäß Bild 3.2.2c lassen sich die Filterschaltungen 
in Bild 3.2.20 realisieren: 
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lo e CL L 


| 
2C, 


Ip 
am 
o— — 


a) T-Schaltung 


b) m-Schaltung 


Bild 3.2.20 Bandpaß in T- und z-Schaltung für Ro = 600.0, fyı = 1kKHz, fga= 1,2 kHz 


In den folgenden Bildern 3.2.21 bis 3.2.26 sind die berechneten Größen für den Bandpaß in Ab- 


hängigkeit von der Frequenz dargestellt. 


2,0 
kn 


06 08 10 1,2 1,4KH21,6 
—f 
Bild 3.2.21 Wellenwiderstand eines Band- 
passes in T-Schaltung (Ro = 600 0, 
fg = 1 kHz, fg2 = 1,2 kHz) 


10 
Np 


0 
1,0 1,5 kHz 2,0 
— f 
— nd 
-® -95 Ada 35 7 
—— N 


Bild 3.2.23 Wellen- und Betriebsdämpfung 
eines Bandpasses (Ro, = 600.02, f,, = 1 kHz, 
fga = 1,2 kHz) 


"06 08 10 12 16KHZ1,6 
— 
Bild 3.2.22 Wellenwiderstand eines Band- 
passes in r-Schaltung (Ro = 600 92, 
fgı = 1 kHz, fg2 = 1,2 kHz) 


0 0,5 1,0 1,5 kHz 2,0 


—f 


H— 1 

-o© -95 -104 35 7 
Bild 3.2.24 Phasengang der Wellen- und 
Betriebsübertragungsfunktion eines Band- 
passes (Ro = 600 22, f,} = 1 kHz, 
fg2 = 1,2 kHz) 
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- 500 
0 250 500 70 N 100 
— Re l2,} 
Bild 3.2.25 Eingangswiderstand eines Band- D 50 1000 1500. 2000 
passes in T-Schaltung {Ro = 600 02, —- Reiz.) 


fgı = 1 kHz, fga = 1,2 kHz) 
Bild 3.2.26 Eingangswiderstand eines Band- 
passes in n-Schaltung (Ro = 600 8, 
fgı =1 kHz, fga = 1,2 kHz) 


d) Bandsperre 


Eingaben: 600 STOA (Kennwiderstand R, in Ohm) 
1000 STOB (Untere Grenzfrequenz f,, in Hz) 
1200 STOC (Obere Grenzfrequenz f,> in Hz) 
900 STOD_ (Startfrequenz f in Hz) 
1 STOE (Anpassungsfaktor q bzw. s} 
-1I STOI  (Schrittweite AR für den zyklischen Rechnungsablauf) 


Start [Le] Programmvorbereitung für T-Schaltung (Flag 1 wird gelöscht) 


Start [D]: 


Ergebnisse: L, 0,016 H 
1,326 uF | 

0,477H 
44,209 nF 


($) 
» 
l) 


Halbglied 


r 
o 
l 


[$} 
(22 
n 


2. Zyklus 3. Zyklus 4. Zyklus 


-0.538 k 
1111111111. 
532.272 505.590 701.478 


aw [NP] 0.000 0.000 1.992 
vw 54.973° -65.158 180.000 


z. [9] 511.050 452.860 697.445 

Ye -6.439° 17.473° 272.106° 

as [Np] 0.005 0.012 1.329 2,797 
Ya 55.165° -65.475° 171.404° 136.501° 


Rechenzeit rd. 20 Sekunden je Zyklus 
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Start Flle} Programmvorbereitung für m-Schaltung (Flag 1 wird gesetzt} 


Start [Do]: 


Ergebnisse: L, = 0,016H | 
C,= 1,326 uF 
[= Dam Halbglied 


C, = 44,209 nF 
0.462 -0.536 -1.538 -2.538 
1311111111. 1111111111. 
Zwin® 676.346 712.039 -513.202 -257.158 
ay inNp 0.000 1.992 3.167 
54.973 -65.158 180.000 

i 704.432 794.947 516.170 257.789 

-270.147° 


6.439° -7.473° -272.106° 
136.50 1° 


0.005 
55.165° 


171.404° 


Rechenzeit rd. 20 Sekunden je Zyklus 


Mit den berechneten Kenndaten für ein Halbglied gemäß Bild 3.2.2d lassen sich die Filterschaltungen 
in Bild 3.2.27 realisieren: 


La la 
Ca Co 
Ip 
2 
Zi T 
a) T-Schaltung b) r-Schaltung 


Bild 3.2.27 Bandsperre in T- und n-Schaltung für Ro = 600%, fyı = 1kHz, fa = 1,2 kHz 


In den folgenden Bildern 3.2.28 bis 3.2.33 sind die berechneten Größen für die Bandsperre in 
Abhängigkeit von der Frequenz dargestellt. 
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0 0,5 1,0 1,5kHz 2,0 
— f 
Bild 3.2.28 Wellenwiderstand einer 
Bandsperre in T-Schaltung 
{Ro =600. 02, fyı = IkKHz, fg2 = 1,2 kHz) 


10 
Np 


—-—o 
un S 


OD 05 910 A1SKHz 20 


—f 


0092 01 02 1,-1 -029 .-014 


m -m ——() 
Bild 3.2.30 Wellen- und Betriebsdämpfung 
einer Bandsperre (Ro = 600 92, f,, = I1kHz, 
fga = 1,2 kHz) 


250 500 750 2 1000 
—- reelle} 


Bild 3.2.32 Eingangswiderstand einer 
Bandsperre in T-Schaltung (Ro = 600 0, 
fgı = 1 kHz, fg2 = 1,2 kHz) 
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9 0,5 1,0 1,5KHz 2,0 


BE 
Bild 3.2.29 Wellenwiderstand einer 
Bandsperre in n-Schaltung 

(Ro = 600.9, fgı = IKHz, fga = 1,2 KHz) 


-180 
0 0,5 10 j 1,5 kHz 2,0 


— 


002 09n 1 |-1 -029 -014 


0 —— N 
Bild 3.2.31 Phasengang der Wellen- und 
Betriebsübertragungsfunktion einer 
Bandsperre (Ro =600.02, f.ı = 1 kHz, 
fg2 = 1,2 kHz) 


0 500 1000 1500.02 2000 
—e Re [2.} 
Bild 3.2.33 Eingangswiderstand einer 
Bandsperre in r-Schaltung (Ro = 600 0, 
fgı = 1 kHz, fg2 = 1,2 kHz) 


3.3 Aktive Filter 


Im Niederfrequenzbereich erfordert der Aufbau von Filtern mit passiven Bauelementen in Form 
von LC-Reaktanzvierpolen große Induktivitäten, die nur schwierig zu realisieren sind. Die Verwen- 
dung von Induktivitäten läßt sich jedoch umgehen, wenn man RC-Schaltungen mit einem gegen- 
gekoppelten Operationsverstärker als aktives Bauelement kombiniert. Derartige Schaltungen werden 
als aktive Filter bezeichnet. Im folgenden wird für die drei Grundtypen aktiver Filter: „‚Tiefpaß, 
Hochpaß und Bandpaß’’ das Übertragungsverhalten beschrieben sowie ein Programm zur Berech- 
nung der Sprungantwort und des Frequenzganges entworfen [16, 17, 18]. 

Der Operationsverstärker als aktives Glied der Filterschaltung wird als Elektrometerverstärker, d.h. 
als Verstärker mit extrem hohem Innenwiderstand und definierter Grundverstärkung vo eingesetzt 
(Bild 3.3.1). 


Bild 3.3.1 


Operationsverstärker als Elektrometerverstärker 
mit einstellbarer Grundverstärkung vo 


(3.3.1) 


Für k=O0 erreicht die Grundverstärkung die Leerlaufverstärkung des Operationsverstärkers. Der 
für die Filteranwendung wichtige Einstellbereich der Grundverstärkung liegt im Bereich 3 sk<1 
entsprechend 1 < vg S 3. Durch den hohen Innenwiderstand des Operationsverstärkers wird die 
zu verstärkende Signalquelle praktisch nicht beeinflußt. 


3.3.1 Aktiver Tiefpaß 


Aus der Vielzahl möglicher Realisierungen aktiver Tiefpaßfilter wird als exemplarisches Beispiel für 
die Programmierung ein Filter 2. Ordnung mit Einfachmitkopplung gemäß Bild 3.3.2 ausgewählt 
[19, 20]. Der besseren Übersichtlichkeit wegen wird die interne Gegenkopplung des Operationsver- 
stärkers zur Realisierung der einstellbaren Grundverstärkung vo gemäß Bild 3.3.1 in den Filter- 
schaltungen nicht mehr besonders eingezeichnet. 


Bild 3.3.2 


iR 
| Sa ) 
Aktives Tiefpaßfilter mit Einfachmitkopplung 
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Bei vorgegebener Eingangsspannung u; lassen sich für die beiden unbekannten Ströme i,, ia und 
für die unbekannte Ausgangsspannung u, folgende drei Maschengleichungen formulieren: 


ueuR+l fi2dt+u2 (3.3.2) 


B 8 u2 
un Rttiiiz) Rt (3.3.3) 


Te u2 
fa, -i,) dt= Vo (3.3.4) 


Durch Differentiation der GIn. (3.3.2 und 3.3.4) und Umordnung der Gl. (3.3.3) erhält man: 


u Hactusu, (3.3.5) 

in DR-i Rt =u (3.3.6) 
1 2 

ne -iaa-n=0 (3.3.7) 


Aus diesem Differentialgleichungssystem erster Ordnung erhält man nach Elimination der Ströme i 
i, und iz eine Differentialgleichung zweiter Ordnung für die Ausgangsspannung u>: 


" ‚1 
u2 = | oo (vo-3) -U2+tvo un | (3.3.8) 
mit 
1 @o 
wone: fo = = (3.3.9, b) 


Zur Lösung dieser Differentialgleichung müssen die Anfangsbedingungen u, (t =+0) und 
u; (t= +0) bestimmt werden. Für die sprungförmige Eingangsspannung u; gilt: 
| U, fürt20O 


u, (t) = | Ofürt<o (3.3.10) 


Unter der Voraussetzung, daß die beiden Kondensatoren zur Zeit t=0 ungeladen sind, gilt: 


u(t=+0)=vou., (0) =0 (3.3.11) 
Aus der GI. (3.3.7) folgt für die Ableitung der Ausgangsspannung: 
‚ Vo. _, 
Uz = c (ir = ig) (3.3.12) 


Da zur Zeit t=0 wegen uz (0) =0 der Ausgang wie kurzgeschlossen wirksam ist, muß i} = ia sein, 
so daß auch u, zur Zeit t=+0 gleich Null sein muß (der zweite Eingangswiderstand liegt dem 
ausgangsseitigen Kurzschluß parallel}: 


u, (t=+0)=0 (3.3.13) 


Mit den formulierten Anfangsbedingungen kann das in Band 1 dieser Reihe angegebene Programm 
zur Lösung von Differentialgleichungen 2. Ordnung y’' = f(x, y, y’) herangezogen werden. 

Der Frequenzgang ergibt sich durch Transformation der Differentialgleichung GI. (3.3.8) in den 
Bildbereich mit p = jw: 


1 
p?u, = ws Ip u (vo -3) ua + voun| (3.3.14) 
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Damit ergibt sich für den Frequenzgang der Verstärkung v: 


ce 
nD 


vo 
v- rg Tr For pe (3.3.15) 
1- (ao +15, 3 "vol 


Zur Normierung wird die Verstärkung v auf die Verstärkung v{w=0) = vo bezogen: 


Var - 1 (3.3.16) 
w . WW 
I +0,83 -W) 
Die normierte Verstärkung V wird in der logaritnmischen Maßeinheit Dezibel angegeben: 
V*=20db-IgV (3.3.17) 


Der Filtertyp wird durch die Grundverstärkung v9 bestimmt. Für vg = 3 ergibt sich eine unge- 
dämpfte Schwingung mit der Kennfrequenz fg = In Abhängigkeit von der Grundverstärkung 
lassen sich folgende Filtertypen realisieren: 


1 
an AC' 


vo = 1,268 Bessel-Filter 

vo= 1,586 Butterworth-Filter 

vo = 2,234 3 db-Tschebyscheff-Filter 
vo=3 ungedämpfte Schwingung 


3.3.2 Aktiver Hochpaß 


Das im vorhergehenden Abschnitt beschriebene 
Tiefpaßfilter läßt sich in ein Hochpaßfilter um- 
wandeln, indem man die Widerstände und Kon- 
densatoren gegenseitig austauscht (Bild 3.3.3). 


Bild 3.3.3 Aktives Hochpaßfilter mit 
Einfachmitkopplung 


Analog zu den Gin. (3.3.2 bis 3.3.4) ergibt sich für das aktive Hochpaßfilter folgendes Gleichungs- 
system: 


ut lndetiatu (3.3.18) 
-4[: ser t fin -inar+ 32 (3.3.19) 

u ey c. I, ia Vo 3. 

U = Vo R (ir 7 ia) (3.3.20) 

Durch Differentiation der Gin. (3.3.18 und 3.3.19) und Umordnung erhält man: 
na+lR+Ug= u (3.3.20) 
1 ‚1 e 
hc aotru rn (3.3.22) 
HR-iaR-u2..=0 (3.3.23) 
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Aus diesem Differentialgleichungssystem erhält man nach Elimination der Ströme iz, iz und ia 
eine Differentialgleichung 2. Ordnung für die Ausgangsspannung u>: 
ei 2[. 1 vo 7] 
= 22 -3)-u+7Zu 3.3.24 
u2=@9 [% Far vo-3)-W Pr: 1 { ) 
mit 
1 


= ‚2.25 
ac (3 } 


Wo 


Die GI. (3.3.24) unterscheidet sich von der entsprechenden Gl. (3.3.8) lediglich in dem letzten 
Term durch = gegenüber u, in Gl. (3.3.8). Dadurch kann zur Berechnung der Sprungantwort des 
w 


Hochpaßfilters ein strukturell gleiches Programm wie zur Berechnung der Sprungantwort des Tief- 
paßfilters benutzt werden. Wesentliche Unterschiede ergeben sich jedoch bei den Anfangsbestimmun- 
gen uz(t=+0) und uz(t=+0): 

Unter der Voraussetzung ungeladener Kondensatoren zur Zeit t=0 wird der Eingangssprung un- 
mittelbar am parallelen Eingangswiderstand des Operationsverstärkers wirksam, so daß am Ausgang 
der mit vo verstärkte Spannungssprung erscheint: 


u,(t=+0)=voU, (3.3.26) 
Aus Gl, (3.3.22) folgt für die zeitliche Ableitung der Ausgangsspannung: 
v 2v 
Up=% tin u a, (3.3.27) 


Wegen der sprungförmigen Eingangsspannung gemäß GI. (3.3.10) ist die zeitliche Ableitung der 
Eingangsspannung für t= +0 gleich null. 
Da die beiden Kondensatoren wie kurzgeschlossen wirksam sind, gilt für die Ströme i, und iz: 


U 

in = = +ia (3.3.28) 

: U,=u U, vol, U, 

ig = Pe war (1-vo) (3.3.29) 
U 

in = = (2-vo) (3.3.30) 


Gin. (3.3.29 und 3.3.30) in Gl. (3.3.27) eingesetzt ergibt: 


; 1 U, 
u = ac Voll vol RC 2 Vol2 = vo) (3.3.31) 
un (t=+0)= U, woVo lvo -3) (3.3.32) 


Mit den Gin. (3.3.26 und 3.3.32) sind die beiden Anfangsbedingungen zur numerischen Lösung der 
Differentialgleichung GI. (3.3.24) mit dem vorhin schon angesprochenen Lösungsprogramm erklärt. 


Der Frequenzgang ergibt sich wieder durch Transformation der Differentialgleichung Gl. (3.3.24) 
in den Bildbereich mit p = jw: 
p?u - «3 |pu In -3)-up tv u] (3.3.33) 
2 0 2 0, Yo 2*Vo w2 1 Kr 


76 


Für den Frequenzgang der Verstärkung ergibt sich: 
u v 
BE re ee, (3.3.34) 
wg2 .wo 
1-(5) tig Wo 3) 
Hier wird zur Normierung als Bezugsgröße die Verstärkung bei = ® mit v(w=«) = vo gewählt. 


v 1 
Ve — — (3.3.35) 
1-(Z) +jg wo 3) 


Für die Umrechnung der normierten Verstärkung in Dezibel gilt wieder Gl. (3.3.17). 


3.3.3 Aktiver Bandpaß 


Ein Bandpaß als selektives Filter entsteht durch die Kombination eines Tiefpaß- mit einem Hoch- 
paßfilter (Bild 3.3.4). 


Bild 3.3.4 
Aktives Bandpaßfilter mit Einfachmitkopplung 


Das Bandpaßfilter läßt sich durch vier Maschengleichungen beschreiben: 


. 17. 
u, Su Rt 2% dt 3.3.36) 
Br I u2 
u7=i, R+tc in-iz<ig)dt+, (3.3.37) 
ERENSTASE: U2 
(41 -ia-ig) 2R = vo (3.3.38) 
uR+tluoR+Ww=u, (3.3.39) 


Durch Differentiation der Integralgleichungen und Umordnung ergibt sich: 


in R tig c =u; (3.3.40) 
“ u2 u 
R + =u; 

iq IACv % u] (3.3.41) 
. } R u2 

IE 2R -ia 2R -i32R - vo =0 (3.3.42) 
a RtiaR +Uu2 =u (3.3.43) 

Nach Elimination der Ströme erhält man: 
”_,2[,- _1 vo, 
U2=@, [u Fr: wo-3)-u+ = ui] (3.3.44) 
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Auch die GI. (3.3.44) stimmt im strukturellen Aufbau mit den entsprechenden Gleichungen für den 
aktiven Tiefpaß Gl. (3.3.8) und den aktiven Hochpaß GI. (3.3.24) überein. Hier erscheint lediglich 
die Eingangsspannung u; in der Form der ersten Ableitung nach der Zeit dividiert durch wo. 

Zur Lösung der Differentialgleichung Gl. (3.3.44) müssen wieder die Anfangsbedingungen u, (t=+0) 
und u; (t= +0) bestimmt werden. Bei energielosen Kondensatoren ist der Eingang des Operations- 
verstärkers als kurzgeschlossen wirksam, so daß die Ausgangsspannung zur Zeit t=+0 auch gleich 
null sein muß: 


u,(t=+0)=0 (3.3.45) 
Aus Gl. (3.3.41) folgt für die zeitliche Ableitung der Ausgangsspannung: 

U SU Vo Up Inc -4Rvo (3.3.46) 

u re i7 R Vo (3.3.47) 
Aus Gl. (3.3.40) folgt für die zeitliche Ableitung des Stromes iz: 

“di _t,,. 1 

un Tc Tor (3.3.48) 
GI. (3.3.48) in Gl. (3.3.47) eingesetzt ergibt: 

u; (t=+0)= U, wgVo (3.3.49) 


Mit den Gin. (3.3.45 und 3.3.49) sind die beiden Anfangsbedingungen zur numerischen Lösung der 
Differentialgleichung Gl. (3.3.44) erklärt. 

Der Frequenzgang ergibt sich wieder durch Transformation der Differentialgleichung Gl. (3.3.44) 
in den Bildbereich mit p = jw: 


2... 00452 ı 3), - „BL 
P’ua= ws E U2 Wo-3)-U*+tvg 2 un] (3.3.50) 
Für den Frequenzgang der Verstärkung ergibt sich: 
u v 
v- = rn en (3.3.51) 
' BEREALE Pia. 
3-votj a m» 


Hier wird zur Normierung als Bezugsgröße die Verstärkung bei der Bandmittenfrequenz w = wg 


.. _YO 5 
mit 3-0 gewählt. 
V 1 
V= vo (3-vo) = Tr: (3.3.52) 
wo w 
1+ji 
J 3 Vo 


Für die Umrechnung der normierten Verstärkung in Dezibel gilt wieder Gl. (3.3.17). 
Die Bandgrenzen sind definiert als diejenigen Frequenzen, bei denen die Verstärkung auf das 


Zzfache (entsprechend 3 db) abgesunken ist. Setzt man hierzu den Imaginärteil im Nenner von 


G!. (3.3.52) gleich +1, so ergibt sich aus der Differenz zweier quadratischer Gleichungen als Band- 
breite B für den aktiven Bandpaß: 


B=fo(3- vo) (3.3.53) 


Im ungedämpften Fall mit v, = 3 geht die Bandbreite des selektiven Filters gegen null. Dieser 
Sachverhalt ist anhand der Ergebnisse in Bild 3.3.12 noch besonders zu beobachten. 
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3.3.4 Programmbeschreibung „Sprungantwort aktiver Filter” 


Die Anweisungsliste zur Berechnung der Sprungantwort aktiver Tiefpaß-, Hochpaß- und Bandpaß- 
filter ist in Tabelle 3.3.1 angegeben. Als Startadressen sind die Label A, B und C vorgesehen. Zu- 
nächst werden die jeweiligen Anfangsbedingungen u; {t =+0) und uz (t = +0) berechnet. Diese 
Berechnungen sind mit dem Sprung nach Label O in Zeile 033 abgeschlossen. Hierzu sind in STO A 
die Grundverstärkung vo, in STO B die Kennfrequenz f, und in STO C die Höhe des Eingangs- 
sprungs U, abgespeichert. Die gewählte Schrittweite At wird in STO D vorgegeben. Das Programm 
besteht von Label b in Zeile 037 bis zur Zeile 185 aus dem allgemeinen Programm zur Lösung von 
Differentialgleichungen 2. Ordnung aus Band 1 dieser Reihe gemäß der dortigen Tabelle 4.5.2. 
Lediglich zur Normierung der Ausgangsspannung auf vo sind die Anweisungen 146 und 147 einge- 
fügt. Die im Unterprogramm Label a von Zeile 150 bis 185 eingebaute automatische Schrittweiten- 
steuerung erweist sich hier als besonders vorteilhaft. Die Differentialgleichung ua= fit, U2, Up} ist 
in dem Unterprogramm Label E von Zeile 186 bis Zeile 204 eingebaut. Der Funktionswert f (u,) 
wird über STO 9 zur Verfügung gestellt. Das Unterprogramm Label c dient zur Berechnung der 
Kreisfrequenz wo (in STO 8 abgespeichert) aus der in STO B abgelegten Kennfrequenz fo. 


3.3.5 Programmbeschreibung ‚„Amplituden- und Phasengang” 


Die Berechnung von Amplituden- und Phasengang der drei aktiven Filterarten Tief-, Hoch- und 
Bandpaß wird in einem Programm gemäß Tabelle 3.3.2 zusammengefaßt. Bei Start über Label A 
wird ein Tiefpaßfilter gemäß Gl. (3.3.16), bei Start über Label B wird ein Hochpaßfilter gemäß 
GI. (3.3.35) und bei Start über Label C wird ein Bandpaßfilter gemäß GI. (3.3.52) berechnet. Zur 
Kennzeichnung der jeweiligen Filterart werden unmittelbar nach dem Programmaufruf die beiden 
Flags O und 1 in folgende definierte Zustände gebracht: 


TP HP BP 


FlagO 1 16) 1 
Flag 1 1 1 0 


Der mit 1 bezeichnete Zustand bedeutet Flag x gesetzt und der mit O bezeichnete Zustand bedeutet 
Flag x nicht gesetzt. 

Das gemeinsame Unterprogramm beginnt mit Label 3 in Zeile 016 und endet mit der Rücksprung- 
anweisung in Zeile 078. Vor dem Programmaufruf muß der Verstärkungsfaktor v, in STO A und 
die Kennfrequenz fo in STO B eingegeben werden. Für zyklisch fortschreitende Berechnungen 
kann in STO D die Schrittweite Af eingegeben werden. Als Startfrequenz wird der im X-Register 
vorhandene Wert genommen. In den Zeilen 017, 018 wird der jeweils aktuelle Frequenzwert ange- 
zeigt bzw. ausgegeben. Die Verstärkung wird als Amplitudenverhältnis u„/u, in einer kurzen Pause 
in Zeile 066 angezeigt. Anschließend wird die normierte Verstärkung in db berechnet und in 

Zeile 070 ausgegeben. Der Phasenwinkel in Grad wird in Zeile 073 ausgegeben. Mit Hilfe der 
Speicherregisterarithmetik-Anweisung in Zeile 076 wird der um Af erhöhte Frequenzwert in 

STO O gebildet und der nächste Berechnungszyklus gestartet. 
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Tabelle 3.3.1 Anweisungsliste ‚„Sprungantwort aktiver Filter“ 


654 
655 
[e7 
057 
58 
e59 
ce 
a6l 
a6? 
&62 
664 
Bes 
[73 
ter 
08 
x [135 
RELh ara 

x Al 
srüor er 

[) &73 
5798 874 
ers 
676 
877 
ers 
er8 
e8a 
881 
as: 
a8? 
aad 
as8 
BREe 
es7 
ade 
869 
e36 
es1 
ea2 
693 
94 
[5 
096 
087 
096 
098 
106 
ı@1 
i62 
102 
184 
105 
10€ 


Tiefpaß au! ALELA 


EuEe 


“LBLE 
003  RÜLÄ 
aa STOE 
611 3 
iz = 

ul? 685Be 
014 
u18 


Hochpan v&s 


x 
RELC 
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ST 
[3-72 
Fa? 
STOb 
aLELE 
RÜLT 
€ 


E 
RCLL 


RÜLS 


UP 


187 
18€ 
108 
11&e 
111 
112 
113 
114 
115 
116 
117 
ile 
113 
126 
121 
122 
t23 
124 
125 
126 
ler 
12€ 
129 
138 
121 


_ 


nn 
a a Cu fa fa 


NEN at 


149 
158 
151 
152 
153 
154 
155 
156 
157 
158 
153 


FSE 
FRTS=>t 
RcL2 


FSE=$>W 


FÜLH 


PRTx=D „2 


ETük 
*LBLa 
RCL? 
RCLe 
=Y? 
6T0Z 
RÜLE 
RCLS 


up 


Up 


168 ABS 
161 . 
162 i 
163 FE) 
164 8Y? 
165  6TÜl 
166 . 
167 5 
168 $F8 
—163 6702 
178 «LBLi 
171 R+ 
172 . 
173 6 
174 5 
175 Kr? 
176 6T02 
177 2 
178 sr& 
179 xLBL2 
188 RCLD 
181 Fo? 
182 “ 
183 PSE=DAtL 
184 STOD 
185 RTN 
186 &LBLE 
187 Ry 
188 CHS 
189 RCLA 
198  RCL3 
191 x 
192 + 
193 Key 
1934 RELE 
195 * 
196  RCLA 
197 3 
198 = 
199 x 
288 + 
261  RCLB 
282 x® 
283 x 
204 RTN 
2085 sLBLc 
2086  RCLB 
287 2 
288 x 
209 Pi 
218 x 
2il RTN 


3.3.6 Test- und Anwendungsbeispiele 


Es sei die Sprungantwort, der Amplitudengang und der Phasengang aktiver Tiefpaß-, Hochpaß- und 
Bandpaßfilter mit der Kennfrequenz fg = 100 Hz für die Filtertypen Bessel mit vo = 1,268, Butter- 
worth mit vo = 1,568, Tschebyscheff mit vo = 2,234 und für den ungedämpften Fall mit vo =3 
zu berechnen. 

Eingaben für die Berechnung der Sprungantwort mit dem Programm nach Tabelle 3.3.1: 


1,268 

vo= sn in STOA 
3,000 

fo=100Hz in STOB 

U,=1V  insToc 


At =0,001s in STOD 


Start la]: Tiefpaß 


Ergebnisse: 


Bessel-Filter Tschebyscheff-Filter ungedämpftes Filter 
1,268 STO A 2,234 STOA 3STOA 


0,0003 s?) 0,0005 s?) 0,0010 s® 
0,0113 0,0452 0,1909 


0,0005 s 0,0010 s 0,0020 s 
0,0412 0,1639 0,6905 


0,0008 s 0,0015 s 0,0030 s 


0,0846 0,3299 1,3081 


0,0010 s 0,0020 s 0,0040 s 
0,1375 0,5186 1,8080 


0,0013 s 0,0025 s 0,0050 s 
0,1963 0,7089 1,9995) 
0,0015 s 0,0030 s 0,0060 s 
0,2584 0,8838 1,8095 


1) Durch die automatische Schrittweitensteuerung wurde die in STO D mit 0,001 vorgegebene Schrittweite 


durch zweimalige Halbierung auf 0,00025 reduziert. 
2) Durch die automatische Schrittweitensteuerung wurde die in STO D mit 0,001 vorgegebene Schrittweite 
durch einmalige Halbierung auf 0,0005 reduziert. 
Hier wurde die vorgegebene Schrittweite unverändert übernommen. 
Der Rechenfehler im ungedämpften Fall beträgt nach dem Sten Berechnungszyklus 0,05 %. Dieser ist durch 
Vergleichsrechnungen mit kleinerer Schrittweite nachweisbar. 


3 
4 


Die Ergebnisse sind in Bild 3.3.5 mit den Kennlinien für Bessel-Filter a) Butterworth b) 
Tschebyscheff c) und den ungedämpften Fall d} aufgetragen. 
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Tabelle 3.3.2 Anweisungsliste ‚Amplituden- und Phasengang aktiver Filter‘ 


Tiefpaß el »LELA 848 Fi? 
882 Ssf& 041 cCTöe 
883 Sri ber 670] 
084 6SE3 843 #LBL& 
885 EToh 844 x 
HochpaN BU6 wLELb 645 1 
897 CFö 0846  RÜl 
888 <Fi 647 x2 
089 65SB3 048 - 
Bi&a___6CTOE 842 aLBLI 
00» 
e12 sF& 851 178 
e13 CFi 852  ENTt 
014 6563 0533  ENTT 
815 __ETOC 0654  RÜLA 
up 615 #LBLS 855 x 
el? PSE 05€ FSE$DUL 
ie  PRTS=DF 57 R+ 
819 $TO8 858 Fi? 
6828 RCLB 859  ETOA 
621 ® 868 3 
822 Fo? 861 ENTt 
023 6T08 862 RCLA 
[,: 178 863 = 
0825 «LBLd 864 x 
826 $Tol 65 #LBLE 
627 Fi? 866 LOG 
826 £6TÜ8 867 74 
e23 ENT 868 B 
03 1X 863 x 
e31 - er8 PRIX DU" 
6832 xLEL& 871 KEY 
833 2 872 cHh® 
934 FÜLR 873 FRTA=DY 
635 = 874 SFC 
036 Fe? 075  RCLL 
037 gToe 876 ST+r& 
[es CHS 077 Rlcd 
633 xLBLE ars RTN 
878 R-& 


Kurve a) vo = 1,268 (Bessel-Filter) 

Kurve b) vo = 1,568 (Butterworth-Filter) 
Kurvec) vo= 2,234 (Tschebyscheff-Filter) 
Kurve d) vo = 3,000 (ungedämpfte Filter) 


Bild 3.3.5 Normierte Sprungantwort aktiver Tiefpaßfilter mit fo = 100 Hz 
und der Verstärkung vo als Parameter 
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Start [B]: Hochpaß 


Die Ergebnisse sind in Bild 3.3.6 für die verschiedenen Filtertypen aufgetragen. 


Bild 3.3.6 


paßfilter mit fo = 100 Hz und der Ver- 
stärkung vo als Parameter (Parameter 


-10 siehe Bild 3.3.5) 
0 2 & 6 ) 10 12ms 14 


Start [c}: Bandpaß 


Die Ergebnisse sind in Bild 3.3.7 für die verschiedenen Filtertypen aufgetragen. 


Bild 3.3.7 


paßfilter mit fy = 100 Hz und der Ver- 
stärkung vo als Parameter (Parameter 
siehe Bild 3.3.5) 


2% 8 8 % Wmık 


Die Berechnung von Amplitudengang und Phasengang erfolgt mit dem Programm gemäß 
Tabelle 3.3.2. Hierzu sind folgende Eingaben erforderlich: 


1,268 
Grundverstärkung Vo = Re | in STOA 
3,000 ! 
Kennfrequenz fo = 100 Hz in STOB 
Schrittweite Af= 25Hz in STOD 
Anfangswert der Frequenz = 25Hz im X-Register 


Normierte Sprungantwort aktiver Hoch- 


Normierte Sprungantwort aktiver Band- 
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Start : Tiefpaß 


Ergebnisse: 


1,268 STO A 2,234 STOA 
25 Hz 
-0,28 db 0,38 db 
-24,79° -11,54° 
50 Hz 50 Hz 
-1,18 db 1,49 db 
-49,11° -27,05° 
75 Hz 75 Hz 
vr -2,74 db 2,83 db 
fr) -71,39° -52,71° 
f 100 Hz 100 Hz 
vr -4,77 db 2,32 db 
fr) -90,00° 
f 125 Hz 125 Hz 
v* -6,99 db 
104,56° 


25 Hz 
0,56 db 
0° 

50 Hz 
2,50 db 
0° 

75 Hz 
7,18 db 
0° 


100 Hz 


Error wegen v>=® 
Nach CL x Eingabe 125 


Start [A]: 


125 Hz 
5,00 db 
180° 


Die Ergebnisse der Frequenzgangberechnungen für die verschiedenen Typen aktiver Tiefpaßfilter 
sind als Diagramme in Bild 3.3.8 für den Amplitudengang und in Bild 3.3.9 für den Phasengang 


dargestellt. 


Die entsprechenden Ergebnisse für den aktiven Hochpaß und den aktiven Bandpaß werden durch 


Aufruf der Programmadressen Label B und Label C gewonnen. 


Die Ergebnisse sind für den Hochpaß in den Bildern 3.3.10, 3.3.11 und für den Bandpaß in den 


Bildern 3.3.12, 3.3.13 dargestellt. 


20 
db 


ug 
150 Hz 200 


0 50 100 


—f 
Bild 3.3.8 Normierte Verstärkung V* 
eines aktiven Tiefpaßfilters mit 
fo = 100 Hz und der Verstärkung 
vo als Parameter (Parameter siehe Bild 3.3.5) 
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0 50 100 


—f 


150 Hz 20 


Bild 3.3.9 Phasengang der Verstär- 
kung eines aktiven Tiefpaßfilters 
mit fo = 100 Hz und der Verstär- 
kung vo als Parameter (Parameter 
siehe Bild 3.3.5) 


0 50 100 150 Hz 200 


— f 


Bild 3.3.10 Normierte Verstärkung 
eines aktiven Hochpaßfilters mit 

fo = 100 Hz und der Verstärkung 
vo als Parameter (Parameter siehe 
Bild 3.3.5) 


10 1S0Hz 200 


—— 


Bild 3.3.12 Normierte Verstärkung 
eines aktiven Bandpaßfilters mit 

fo = 100 Hz und der Verstärkung vo 
als Parameter (Parameter siehe 

Bild 3.3.5) 


3.4 Übertragungsverhalten elektrischer Netzwerke 


3.4.1 Übertragungsfunktion 


0 50 100 150 Hz 200 
—— f 
Bild 3.3.11 Phasengang der Ver- 
stärkung eines aktiven Hochpaßfilters 
mit fo = 100 Hz und Verstärkung vo 
als Parameter (Parameter siehe 
Bild 3.3.5) 


0 50 100 150 Hz 200 
—-f 
Bild 3.3.13 Phasengang eines aktiven 
Bandpaßfilters mit fo = 100 Hz und 
der Verstärkung v9, als Parameter 
(Parameter siehe Bild 3.3.5) 


Die Übertragungsfunktion F({p) ist die Laplace Transformierte der Gewichtsfunktion g(t) eines 
linearen zeitinvarianten Übertragungsgliedes [21, 22, 23]. Das Faltungsintegral der Gewichtsfunk- 
tion mit der Eingangsgröße u(t) ergibt die Ausgangsgröße v{t): 


ı 


v(t) = [st-n ulr)dr=git) * ut) 


(3.4.1) 
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Durch Transformation in den Bildbereich geht das Faltungsprodukt in eine Multiplikation der 
Laplace-Transformierten über [24]: 


v({p) =F(p) -u(p) (3.4.2) 


Da das Übertragungsverhalten linearer zeitinvarianter Übertragungsglieder durch die Übertragungs- 
funktion vollständig beschrieben wird, kommt der Berechnung dieser Funktion besondere Bedeu- 

tung zu. Sie läßt sich mit Gl. (3.4.2) als Quotient der Laplace-Transformierten der Ausgangsfunk- 

tion v(p) und der Eingangsfunktion u(p) bestimmen: 


Flip) =- 4 
(p=- (3.4.3) 


Die Übertragungsfunktion F(p) auf der imaginären Achse p = jw der komplexen p-Ebene heißt 
Frequenzgang F(jw) [25]. Der Frequenzgang ist besonders bei elektrischen Netzwerken von Be- 
deutung, da er das Übertragungsverhalten des Netzwerkes durch Amplitudengang und Phasengang 
beschreibt. Beide Kennlinien können mit Hilfe von Frequenzanalysatoren aus einem erregenden 
kontinuierlichen Frequenzspektrum aufgenommen werden. 


Bild 3.4.1 
RLC-Netzwerk 


Die programmierte Berechnung einer Übertragungsfunktion soll an einem RLC-Netzwerk 
(Bild 3.4.1) am Beispiel der Spannungs-Übertragungsfunktion F,{p) aufgezeigt werden [26]: 
Für die vier Maschen des gegebenen Systems lassen sich folgende Kirchhoffschen Gleichungen 
angeben: 


t 


-uo+ Ron + = |iy-ig)dr =0 
Cı. 
0 

t t 
4 K -i fürs I ig) + Raflia -i fe -ig)dr= 0 
a da la TEN a Nas 

{) 0) 

t 
1 a r 
na [ia-i2Jar+ Rı ig =0 (3.3.1) 
2. 

o 

r I a Ce 
Rzigt Rolie-iz)+L lie mia) =0 
R; ig = u4{t) 


Zur Ermittlung der Übertragungsfunktion des Systems wird das Integrodifferential-Gleichungs- 
system Gl. (3.3.1) aus dem Zeitbereich in den Bildbereich mit der komplexen Variablen p = jw 
transformiert. Unter der Voraussetzung, daß das Netzwerk zur Zeit t=0 keine Energie gespeichert 
hatte, d.h. alle Anfangswerte verschwinden, gilt im Bildbereich: 
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o° 
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[g) 
- 
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1 
{ +( +Lp+R,+ I 1 I2-(R,+Lp) I 0 
c ı*\c 2 pP)? c 3 2 a” 
I e % (3.3.2) 
= + {R + )' =0 
GP > 1 p/”? 
-(Ra+Lp} 1a +(Rz+Ra+tLpll; = 0 
Für den Strom Iz gilt im Bildbereich: 
ee, +tLp+R a -{Ra + Lp) 
Cıp Cıp ? Cp . 
1 
Us {p) 0 ee 0 
op cp 
0 -(R3+Lp) R,+R3+ Lp 
Pl = m (3.3.3) 
1 1 
Rot <— m een 0 0 
° CP Cıp 
— irre ie 
cp cı er’ cp 27 -Pp 
1 1 
0 ee Rı+t —— 0 
Czp ' &p 
0 - (Rz + Lp) 0 Rz+ Rz+Lp 
Für die Übertragungsfunktion F,{p) gilt: 
U, {p) l3 {p} Dz 13 
F EEE LELEL =R 3.3.4) 
De m 
R2 L 
1+ = P% 
Fl a ———— (3.3.5) 
Cı1C2Pp 
RR N 
2 L 
F,(p) re le... ıy 
‘pl = Ro+R2 R Ro‘ R 
o*R2 Ra 0 j LRRRFIE: KIEL © 
i Rı lite )erltlonst ats) * Rolcı*Ca 
(3.3.6) 


=. Li. 2| [2e Ro 
+ + RzCı (+ sl R262}+P IERSEZEZ TE 


C2|*Ro R2C, C2}+p®LRoc,C, 
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3.4.2 Eingangsimpedanz 


Neben der Übertragungsfunktion ist der Frequenzgang der Eingangsimpedanz eines elektrischen 
Netzwerkes eine wichtige Beschreibungsgröße für das Betriebsverhalten des Netzwerkes. Daher 
wird auch eine allgemeine Beziehung für die Berechnung und Programmierung der Eingangsimpe- 
danz Z, hergeleitet. Ersetzt man in dem Gleichungssystem für den Bildbereich GI. (3.3.2) die 
linke Spalte durch die rechte Seite des Systems, so ergibt sich eine Bestimmungsgleichung für den 
Eingangsstrom I; (p). 


1 1 1 
<—— +Lp+Ra+ — = -(RztLp) 
Pr Cap an 
1 1 
= Ritt =— 0 
Czp "7 &p 
- (Rz +Lp) 0 Rz+Ra+Lp 
I (pl) = mc aaı Were ng: (3.3.7) 
N 
Für die Eingangsimpedanz Z, gilt: 
Uolp)_ Du 
zZ = = — 3.3.8 
Ze(p) (Pl Dzu 
Für die Zählerdeterminante des Eingangsstromes I; (p) gilt: 
Rı Ra | (= ) 2 1 L 
Dzu= ——— I |< +1) C,p+ = C,p? + + —+ p 
2 &,Cap? |\Ra z 2" RıR3 Rı RıRz (3.3.9) 


Rz 


c R L 
+LCıCzp?+L2 pP +R2CıCap?+ rlrn)an+ &chp? 
1 1 3 3 


Damit folgt für die Eingangsimpedanz: 


Ten AS. 1 (3.3.10) 
+p ıl2xc. (er = )|+ R2Cı62 | P®ıcic, 
Der Zähler von Gl. (3.3.10) ist mit dem Nenner aus GI. (3.3.6) identisch. Dies bringt den Vorteil, 
daß die Berechnung dieses Terms sowohl für die Übertragungsfunktion als auch für die Eingangs- 
impedanz in einem gemeinsamen Unterprogramm durchgeführt werden kann. 


3.4.3 Programmbeschreibung „Übertragungsfunktion” 


Das Programm nach Tabelle 3.4.1 gestattet sowohl die zyklische Berechnung der Übertragungs- 
funktion wie auch die zyklische Berechnung der Eingangsimpedanz in Abhängigkeit von der Fre- 
quenz. Zur Berechnung der Übertragungsfunktion wird das Programm über Label A gestartet, 
während die Berechnung der Eingangsimpedanz den Programmstart über Label B erfordert. Das 
Programm umfaßt vier Unterprogramme mit den Einstiegadressen Label a, Label b, Label c und 
Label D. 
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Tabelle 3.4.1 Anweisungsliste ‚„Übertragungsfunktion” 


857 RUL4 


ee2 KCLU 858 x UP 
803 TuE 859  5T07 
004 65Be 868 x 
665 &LELO 861  ST-9 
II 173 d62 Rüin 
687  RELA 663 RÜLi 
605 RCL3 &64 3 
LE * 865  RÜLE 
eie Plid 856 e 
öl x 867 RCLS 
el2 RÜLE @66  RELE 
613 RCLS a9 RÜLÜ 
814 * ara + 
ei 1 el x 
ele + 872 + 
817 *F 673  RELE 
el8 RÜLE a4 RCLI 
als RÜELS 675 + 
a2e EZ Br + 
v21  686« e77 Ri“ 
a22__6TÜ 078 CLE 
Ze 623 _#LELE ar8 x 
0243 RÜLE B86  RüLc 
025  STüs as) x 
26 656% 082 Ri 
827 83 x 
@28 054 - 
a2: 685 Fils 
636 ase ® 
“31 Ber rl 
032 655 65ER 
833 063 Riıe 
034 ese RCS 
B3E EIN °F 
G3E e82 Reli 
B37 6983  RÜLS 
B38 aa4 >F 
838 @95  6SEc 
846 es6__6Tüi 
631 8937 sLöLe 
u42 e38  6SBO 
43 x a9 PRTE=D|F 
034 RÜLE 186 x e 
845 Rlis 1ei PRRÜDY 
846 * 182 SFC 
847 + 183 RELE 
a4 RÜLH 184 ST+8 
049 UP 165 «lLbla 
eSe RÜLz 186  RÜLS 
ast RCLE 187 PRIX =Df 
852 x 188 2 
053  STüe 189 x 
854  RPLLE 118 Pi 
055 x 111 x 
856 + 112  STÜs 


YO NO N N Ta ia am ME Ca a OT TITTEN 


nn a m 
[= RD die a Eal Bu » Sur BE De ar De 


u La hal lu Cl a A 
FM EM DPI te TATEN EI TI RO GI TFA AN ru BD 


am too nm 


TC A Fi 


u Par Pay Pas Pas Pas har Pas Pin Das Dis Dis a me ma m mn nn a m nn ne dh m ne u u be he ne ne en nd hun 
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Für beide Berechnungsvarianten wird der gleiche Eingabedatensatz benötigt. Die Eingabedaten 
sind folgenden Speichern zugeordnet: 


STOA: L inH STOO: Ro, in® 
STOB: C, in F STO1: R, in® 
STOC: GC; in F STO2: RR, in 
STOD: f inHz STO3: Rz in® 


STOE: Af in Hz 


Als Arbeitsspeicher werden nur die Primärspeicher STO 4 bis STO 9 benutzt. Der gesamte 
sekundäre Speicherbereich wird nicht aktiviert. 

Nach dem Programmstart über Label A in Zeile 001 wird die in STO D abgelegte Startfrequenz 
nach STO 8 übernommen. In dem Unterprogramm Label a wird die Kreisfrequenz « gebildet und 
nach STO 4 abgespeichert. Die RTN-Anweisung in Zeile 113 bildet eine gemeinsame Rücksprung- 
Anweisung für die Unterprogramme Label c und Label a. Zwischen den Anweisungszeilen 005 mit 
Label O und 022 mit GTO O wird eine Endlosprogrammschleife gebildet. Innerhalb dieser Schleife 
wird mit dem Aufruf des Unterprogramms Label b in Zeile 114 der Nenner aus Gl. (3.3.6) be- 
rechnet. Die folgenden Anweisungen von Zeile 007 bis 017 dienen zur Berechnung des Zählers. 

Mit dem Aufruf des Unterprogramms Label c wird die komplexe Division durchgeführt und der 
Betrag der Übertragungsfunktion sowie der Phasenwinkel ausgegeben. Nach Erhöhung der Frequenz 
um Af in STO 8 wird in STO 4 das neue « abgespeichert und der nächste Berechnungszyklus 
abgearbeitet. 

Die Berechnung der Eingangsimpedanz Z, wird über Label B in Zeile 023 gestartet. In der Pro- 
grammschleife zwischen den Anweisungen Label 1 in Zeile 027 und GTO 1 in Zeile 096 wird die 
GI. (3.3.10) ausgewertet. Zur Ausgabe von Betrag und Phasenwinkel wird das Unterprogramm 
Label c in Zeile 095 aufgerufen. 


3.4.4 Test- und Anwendungsbeispiele 


Vor dem Programmstart für die Berechnung der Übertragungsfunktion mit Label A oder der Ein- 
gangsimpedanz mit Label B sind folgende Eingabedaten erforderlich: 


Ohmscher Widerstand Ro in > STOO 
Ohmscher Widerstand R;, in =» STO1 
Ohmscher Widerstand Rz, inf >STO2 
Ohmscher Widerstand Rz in‘ => STO3 
Induktivität L inH >STOA 
Kapazität C,ı nF »>STOB 
Kapazität C; nF >STOC 
Startfrequenz f inHz >STOD 


Frequenz-Schrittweite Af in Hz > STOI 


Als Ergebnisse werden die Frequenz sowie bei Start über Label A die Übertragungsfunktion F,{f) 
nach Betrag und Phase oder bei Start über Label B die Eingangsimpedanz Z.(f} nach Betrag und 
Phase in einer Endlosschleife mit dem Frequenzintervall Af ausgegeben. 


Ro=Rı =R3=R3=600 2 
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600 
600 
600 
600 


02 
0,265 108 
0,265 10 
0 
200 


Eingaben: 


sSTOoo 
STO 1 
sto2 
sSTO3 


STOA 
STOB 
STOC 
sSTOD 
STOI 


Widerstände in Ohm 


Induktivität L in Henry 
Kapazität C, in Farad 
Kapazität C, in Farad 
Frequenz f in Hz 
Frequenzintervall Af in Hz 


Start [A]: Übertragungsfunktion F,, 


Ergebnisse: 

f= 0.00Hz f = 400.00 Hz ft = 800.00 Hz 
F,= 040 F,= 0.37 F,= 0.30 

v = 0.00° y = -31.03° y = -56.53° 
f = 200.00 Hz f = 600.00 Hz f = 1000.00 Hz 
F,= 039 Fh,= 0.33 F,= 0.26 

v = -15.97° v = -44.61° y = -66.89° 


Rechenzeit für einen Zyklus rd. 12 Sekunden 


In Bild 3.4.2 sind die Ortskurven 
für verschiedene Widerstands- 
varianten der Übertragungsfunk- 
tion angegeben. 


Bild 3.4.2 

Ortskurve der Übertragungsfunk- 
tion eines RLC-Netzwerkes gemäß 
Bild 3.4.1 mit Ro = Rı = 600.0, 
L=0,2H,C, =C3 = 0,265 uF 

a) Rz2=00, Rz=»QN 
b)R2=00, R3=6000 

c) Rz =6002,R3 =600% 

d) RR, =6002,R3,=-=Q 


Start E% Eingangsimpedanz Z. 


Ergebnisse: 

fs 0.00 Hz f = 400.00 Hz f = 800.00 Hz 
Z.= 1500.00 0 Z.= 1256.43 0 Z.e= 994.719 
y= 0.00° v = -19.05° y= -23.51° 

ft = 200.00 Hz f = 600.00 Hz f = 1000.00 Hz 
Z.= 1419.62 Z.= 1106.45 0 z.= 915.050 
v = -11,55° vo = -22.49° vo = -23.33° 


Rechenzeit für einen Zyklus rd. 15 Sekunden 
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ın Bild 3.4.3 sind die Ortskurven der Eingangsimpedanz für verschiedene Widerstandsvarianten 
angegeben. 


Bild 3.4.3 

Ortskurven des Eingangswiderstan- 
des eines RLC-Netzwerkes gemäß 
Bild 3.4.1 mit Ro = Rı =600.0, 
L=0,2H, C, = Ca = 0,265 uF 

a) RR=00, Rz3=-»RQ 
b)R3,=082, R3=6000 

c) Rz =6002,R3 =600% 

d) R,=6002,R;=-»0N 


| 
400 800 1200 1600 2000 
—m Re{Ze} 


3.5 Leitungsgleichungen 

Am Beispiel der Leitungsgleichungen läßt sich besonders anschaulich zeigen, daß bei Programmie- 
rung der exakten Vierpolgleichungen [27] alle möglichen Wertebereiche ohne rechentechnischen 
Mehraufwand abgedeckt sind. Dadurch kann ein erstelltes Programm sowohl für nachrichtentechni- 
sche Übertragungswege wie auch für energietechnische Leitungen angewendet werden. 


3.5.1 Berechnungsgrundlagen 


Die Primärform der Leitungsgleichungen lauten: 


U, = U2 cosh Y I+1, Zw sinh yl (3.5.1) 
14 = la cosh yI+ 32 sinh Y (3.5.2) 
ER BA 

mit dem Wellenwiderstand Zu: 

zw) en _ (3.5.3) 
G'+tjwcC' 
mit dem Übertragungsmaß Y: 

Y=VIR +Hj@L’) (G’+jwC”) (3.5.4) 
=(a+jßI=atjb (3.5.5) 


a ist die Übertragungsdämpfung und b der Übertragungswinkel der Leitung. 
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Für die Programmierung wird die Form: 


Zw 
U,=U, (cosh yı+ zZ sinh Y ) (3.5.6) 
= Z & 
22 
L,=l, (cosh yl+ Z= sinn y ı) (3.5.7) 
zu z 
benutzt. 


Für den Reflexionsfaktor gilt: 


22 "Zw _ Zw 


r ———- 


"Zutzu 22 


(3.5.8) 


Das Spannungsverhältnis a, wird in Dezibel angegeben: 
U} 
a,= 20 ig — db (3.5.9) 
U, 


Für die Leistung am Leitungsanfang gilt: 
Sı=Ulı =Pı+jQı (3.5.10) 


Für den Übertragungswirkungsgrad gilt: 


(3.5.11) 


3.5.2 Programmstruktur „Leitungsgleichungen’’ 


Mit Hilfe der nachfolgenden Eingabedaten und Eingabespeicher läßt sich ein Programm zur Aus- 
wertung der angegebenen Gleichungen konzipieren: 


Spannung U; in V am Leitungende => STOA 
Impedanz Z, in 2 am Leitungende => STOB 
Phasenwinkel y in ° der Impedanz Z, > STOC 
Frequenz f in Hz > STOD 
Länge I in km der Leitung > STOE 


Widerstandsbelag R’ in Q/km > STO 1 
Induktivitätsbelag L’inH/km > STO2 
Ableitungsbelag G’inS/km = STO3 
Kapazitätsbelag C’ in F/km > STO4 


Da die Speicher STO 1 bis STO 4 als Eingabespeicher vorgesehen sind, wird während der Berech- 
nung der sekundäre Speicherbereich aktiviert. Dadurch werden die Eingabespeicher nicht mit über- 
schrieben und stehen für wiederholte Programmstarts zur Verfügung. 

Die Ausgabe der Ergebnisse ist wie folgt gegliedert: 


1. Wellenwiderstand nach Betrag und Phase 

2. Reflexionsfaktor nach Betrag und Phase 

3. Übertragungsdämpfung und Übertragungswinkel 

4. Spannung am Leitungsanfang nach Betrag und Phase 


93 


. Logarithmisches Spannungsverhältnis 

. Strom am Leitungsanfang nach Betrag und Phase 
. Wirkleistung am Leitungsanfang 

. Blindleistung am Leitungsanfang 

. Übertragungswirkungsgrad 


oonoa{a 


Falls das Programm innerhalb der Berechnung gestoppt wird, ist darauf zu achten, daß der Winkel- 
modus-Radiant eingeschaltet ist und der sekundäre Speicherbereich aktiv ist. 
Folgende Unterprogramme können auch über einen GSB Aufruf eigenständig benutzt werden: 


Hyperbolischer Cosinus mit komplexem Argument Label a 
Hyperbolischer Sinus mit komplexem Argument Label b 


Komplexe Multiplikation in Polarkoordinaten Label c 
Komplexe Division in Polarkoordinaten Label d 
Komplexe Addition kartesischer Koordinaten 

mit Ergebnis in Polarkoordinaten Label e 


3.5.3 Programmbeschreibung „Leitungsgleichungen” 


Das in der Tabelle 3.5.1 angegebene Programm wird über die Startadresse Label A in Zeile 001 ge- 
startet. Zunächst wird der Winkelmodus auf Radiant umgestellt. Mit den Anweisungen 003 bis 008 
wird die Kreisfrequenz gebildet und nach STO 5 abgespeichert. Mit Hilfe der bezogenen Leitungs- 
daten in den Speichern STO 1 bis STO 4 werden die Längsimpedanz R + jwL und die Quer- 
admittanz G + jwC in Polarkoordinatenform berechnet und nach dem Speicherbereichswechsel 

in Zeile 018 den Speichern STO 0’ bis STO 3’ zugeordnet. In Zeile 027 wird das Unterprogramm 
Label d zur Berechnung der komplexen Division gemäß Gi. (3.5.3) aufgerufen. Der Wellenwider- 
stand wird in Polarkoordinatenform in dem Unterprogramm Label 2 ausgegeben. Anschließend 
wird der Reflexionsfaktor gemäß Gil. (3.5.8) berechnet und in Zeile 062 mit Betrag und Phase aus- 
gegeben. Mit den Anweisungen 063 bis 076 wird die Gl. (3.5.4) ausgewertet und in Zeile 077 die 
Übertragungsdämpfung a sowie der Übertragungswinkel b in Grad ausgegeben. Bei der Auswertung 
der Leitungsgleichung in der Form von Gi. (3.5.6) wird auf das in den Zeilen 042 und 044 nach 
STO 8 und STO 9 abgespeicherte Impedanzverhältnis Z3/Zy zurückgegriffen. In den Zeilen 080 
und 086 werden die Unterprogramme zur Berechnung der hyperbolischen Funktionen mit kom- 
plexem Argument cosh y ! und sinh y I aufgerufen. Die Spannung am Leitungsanfang wird in 
Zeile 098 nach Betrag und Phasenlage, in Bezug auf die Spannung U, als Bezugsachse in der kom- 
plexen Ebene, ausgegeben. Anschließend wird noch das logarithmische Spannungsverhältnis in 
Dezibel ausgegeben. Mit den Anweisungen 109 bis 120 wird die Leitungsgleichung für den Strom I; 
GI. (3.5.7) ausgewertet und das Ergebnis in Zeile 122 ausgegeben. Anschließend werden noch ge- 
mäß den Gin. (3.5.10 und 3.5.11) die Wirk- und Blindleistung am Leitungsanfang sowie der Über- 
tragungswirkungsgrad berechnet und ausgegeben. Das Hauptprogramm endet in Zeile 139. 
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Tabelle 3.5.1 Anweisungsliste „Leitungsgleichungen” 


Start 


dal 
@82 
68: 
004 
865 
Bac 
6a? 
aus 
aus 
ea 
öl 
al 
es 
BE) 
615 


xLELA 


657 
e56 
859 
806 
861 
862 
863 
864 
665 
Bö6 
667 
858 
869 
876 
671 
BTE 
673 
ara 
875 
are 
a7? 
er& 
07: 
asd 
851 
852 
883 
a84 
85° 
886 
ac? 
888 
863 


Rr 113 | STOS 
>F 114 ST«$ 
kt 115 x 
rt 1le | RCLC 
Eiba 5 117 D>R 
ESE2=D Ip, 118 | RiLe 
RCL3 118 | $T:5 
RCLZ 65Ba 
RcLi STxc h 
RCLE eSBE =D Ip 
6SBr 123 | ST-F I 
‘x 124 | RlL? 
Ä 125 | Rüis 
126 FE 
127 | FRTX=DPı 
ı2E | ST#8 
128 Aar 
se | PRTY=DU 
121 SFC 
132 | RiLS 
137 25 
134 
135 
12€ 
137 
x 138 
5 138 
R UP 148 
RÜ 141 
& 132 
& 143 
ä 143 
$ 14 
R 146 
147 
RC: 148 
149 ETn 
€ up 158 #LELE 
RÜi IE1  PRTS UP: 
152 ATS, A:B 
ste u 153 RD 
[SER r) 154 FRT& 
st; u 155 b>R 
5 156 SFÜ 
RCLA 157 RTN 
# UP: 155 «ukie 
LO cosh 159 SFZ 
E rise _C£Tüa UP: 
& UP. | 161 xLEib A+B 
x sinh FE CFE 
PRTX=D a, 163 »#LELG 
sec 164 ENT? 
RCL3 165 CHS 
RÜLE 166 er 
2:71 167 ENT 
RÜLR 166 R+ 


CHs 
xLBL& 
+ 


RT 
sin 


N 


Ft 
AeY 
+F 
RTN 
“LBLı 
zer 
Rs 
x 
Rs 
+ 
Ft 
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3.5.4 Test- und Anwendungsbeispiele 


1. Für eine Fernsprech-Freileitung mit den angegebenen Daten sind die Übertragungsparameter 
zu berechnen [9]: 


Leitungsdaten: R’= 2,16 22/km 
L’= 0,0018 H/km 
G’ = 0,5: 10°8 S/km 
C’ = 0,0068 aF/km 


I = 50 km 
Betriebsdaten: U, = 100V 
Z, = 2000 ed? Q 
f = 800 Hz 


Vor dem Programmstart sind folgende Eingaben erforderlich: 


100 STOA (Spannung U; in Volt) 
2000 STOB (lmpedanz Z, am Leitungsende in Ohm) 
20 STOC (Phasenwinkel der Impedanz Z; in Grad) 
800 STOD (Frequenz f in Hertz) 
50 STOE _({Leitungslänge | in km) 
2,16 STO1 (Widerstandsbelag in S2/km) 
0,0018 STO2 (Induktivitätsbelag in H/km) 
05EEX-6 STO3 (Ableitungsbelag in S/km) 
0,0068 EEX-6 STO4 (Kapazitätsbelag in F/km) 


Start [A]: 


Ergebnisse: Wellenwiderstand Zy, = 521.650 
Phasenwinkel p, = -6.29° 
Reflexionsfaktor r = 0.63 
Phasenwinkel p, = 13.93° 
Übertragungsdämpfung a = 0.11 
Übertragungswinkel b = 50.69° 
Spannung am Leitungsanfang U, = 78.78 V 
Phasenwinkel p,, = 19.25° 
Spannungsdämpfung a, = -2.07 db 
Strom am Anfang I} = 0.15 A 
Phasenwinkel yı, = 78.75° 
Einspeise-Wirkleistung P} = 5.83 W (Minuszeichen bedeutet 


Einspeise-Blindleistung Q; -9.89 var kapazitiv) 


0.81=81% 


Übertragungswirkgrad n 
Rechenzeit rd. 50 Sekunden 
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Bild 3.51 
Frequenzgang der Spannungsdämpfung 
einer Fernmeldeleitung 


In Bild 3.5.1 ist der Frequenzgang der Spannungsdämpfung im Intervall 100 Hz bis 20 000 Hz 
aufgetragen. Typisch für den nicht reflexionsfreien Belastungsfall ist der Wechsel von Span- 
nungsmaxima und Spannungsminima entlang einer Leitung. Dies wird hier durch die Oszilla- 
tion der Spannungsdämpfung zwischen +0,5 db und —9 db sichtbar. Der Bereich negativer 
Dämpfung besagt, daß die Spannung am Leitungsanfang kleiner ist als die Spannung am 
Leitungsende (—- 9,5 db entspricht einem Spannungsverhältnis U, : U, = 1:3). Die mit 
wachsender Frequenz scheinbar dichtere Aufeinanderfolge von Maxima und Minima ist 
durch den logarithmischen Frequenzmaßstab bedingt. 


Für eine 380 kV-Hochspannungsleitung von 400 km Länge sind für ohmsche Belastung der 
Wellenwiderstand und die Übertragungsparameter zu bestimmen [28]: 


Leitungsdaten: R’ = 0,03 2/km 
L’= 1,05 mH/km 
G’ = 0,6 - 1077 S/km 
C’ = 0,0107 F/km 


I = 400 km 
Betriebsdaten: U, = 220.000 V (Strangspannung) 
2, =313 0 (Belastungswiderstand) 
f = 50 Hz (Frequenz) 


Vor dem Programmstart sind folgende Eingaben erforderlich: 


220.000 STOA (Spannung U; in Volt) 
313 STOB (Impedanz Z, in Ohm) 
O0 STOC _(Phasenwinkel der Belastung in Grad) 
50 STOD (Frequenz in Hertz) 
400 STOE (Länge in km) 


0,03 STOI (Widerstandsbelag in S2/km) 

1,05EEX-3 STO2 (Induktivitätsbelag in H/km) 
0,6EEX-7 STO3 _(Ableitungsbelag in S/km) 
0,0107 EEX-6 STO4 (Kapazitätsbelag in F/km) 


9 
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Start [A]: 


Ergebnisse: Wellenwiderstand Zy, 
Phasenwinkel p, 


Reflexionsfaktor r 
Phasenwinkel yo, 


= 313,880 


Übertragungsdämpfung a = 


Übertragungswinkel b 


Spannungsanfang U; 
Phasenwinkel Y,,, 


Spannungsdämpfung a, 


Strom am Anfang I} 
Phasenwinkel y,, 


Einspeise-Wirkleistung P; 
Einspeise-Blindleistung Q} 


Übertragungswirkgrad n 


Rechenzeit rd. 50 Sekunden 


-2,09° 
0,02 
94,40° 
0,02 
24,15° 


= 228.118,69 V 


23,82° 
0,31 db 


= 709,40 A 


24,46° 


161.816.856,3 W = 161,8 MW (je Phase) 
-1.824.064,98 VAr = - 1,8 Mvar (je Phase) 


0,96 96 % 


4 Impulstechnik 


4.1 Zeitverhalten von Übertragungsgliedern 


In der Impulstechnik tritt häufig die Frage nach dem zeitlichen Verlauf der Ausgangsspannung von 
integrierenden oder differenzierenden Übertragungsgliedern bei pulsförmigen Eingangsspannungen 
auf. Zur Beantwortung dieser Frage ist die Übertragungsfunktion als Differentialgleichung zu for- 
mulieren, die dann unter Verwendung der in Abschnitt 1.3.3 angegebenen Programmbausteine für 
periodische Schaltfunktionen mit Hilfe des in Band 1 angegebenen Programms zur numerischen 
Lösung von Differentialgleichungen schrittweise gelöst werden kann. 


4.1.1 Integrierstufe 


Die Kombination des Differentialgleichungsprogramms für Differentialgleichungen erster Ordnung 
mit dem Programmbaustein für die Erzeugung einer Rechteckspannung soll am Beispiel eines Inte- 
griergliedes gemäß Bild 4.1.1 aufgezeigt werden [29]: 

iR 


—— 


oO I O 
ai ST 
Bild 4.1.1 Integrierglied 0 0 


Das Übertragungsverhalten der Schaltung in Bild 4.1.1 kann durch zwei Gleichungen beschrieben 
werden: 


u =iRtWw (4.1.1) 
du, W 

E eu a 

i=C dt "An (4.1.2) 


= u -w (1+ =)| 14.1.3) 


mit: 
u] = U, sign (sin wt) (4.1.4) 


für einen rechteckförmigen Verlauf der Eingangsspannung u.. 
Die Programmieraufgabe beschränkt sich auf den Einbau der Gin. (4.1.3 und 4.1.4) als Unterpro- 
gramm in das gegebene Differentialgleichungsprogramm. 


4.1.2 Programmbeschreibung „Integrierstufe’’ 


Das in Tabelle 4.1.1 angegebene Programm umfaßt bis zur Zeile 110 die Anweisungen zur numeri- 
schen Lösung von Differentialgleichungen 1. Ordnung nach dem Runge-Kutta-Verfahren gemäß 
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Tabelle 4.1.1 Anweisungsliste „Zeitfunktion eines Integriergliedes’’ 


847 RÜLF 033 SF 
682  RCLH @48 2 94 8702 
“83  STol 043 x 095 &LBL! 
ed  RÜLE e5e ST+@ 896 R+ 
ae5 STR 851  RELF 897 P 
Ge &LBLb 852  RCLZ 898 8 
ee? RÜLZ 053 + 899 5 
G66  RCLI 654  RCLI 108  xey? 
@6° cfö 055  RÜLE 161  6T02 
eia  65BE 056 + 162 2 
all RCLE 05?  5ST03 183 SFö 
[oe x e58  STOI 184 &LBL2 
613 STOS 853  6SBE 185 RCLC 
a4  STO8 @68 RELL 186 For 
615 & 861 x 187 x 
ale E @62  ST+8 188 FSE=DAt 
el? REL2 863 5 1893  STOC 
@16 + 064 ST:6 118 __RIN 
015  RCLC 865  RÜLE UP 111 #LBLE 
628 2 B6E  ST+Z 112 3 
e21 F 867 RÜLI 113 x 
022 RCLI 868 SPC 114 Fi 
623 + 869 PSE 115 x 
624 $ST03 878 PRTIS=Dt 116 5 
625 G6SBE 871  RELZ 117 & 
826  RÜLE 872 PSE 116 x 
ver x 873 PRTX=D u, 119 Kl 
023  STLE 874 ___6TOb 128 SIN 
023 2 UP 875 #iBla 121 ABS 
&38 x 876 RL? 122 LSTS 
e31  ST+6 87? RÜLE 123 xra" 
832 4 Fe x=Y? 124 % 
833 ® 08798 6T0E 125 RCi3 
634 RÜL2 eB6 - 126 x 
635 + @81 RÜLE 127 Ne 
636 RÜLZ 882  RCLS 126 RELD 
e37  GSBE 663 - 129 RECLS 
036  RCLÜ 684 3 138 : 
833 x 885 ABS 131 1 
048  STUr ade e 132 + 
e4 Fi? 887 1 123 X 
042 6Tid 886 5 134 - 
043 65Ba 069 Kart 135  RCLO 
044 Fe? e9e  6TOI 136 : 
845 ETÜk [® 3 13? RCLE 
846 «#LBL& 8932 & 133 + 

133 RIN 


Abschnitt 4,5.1 in Band 1. Mit den Unterprogrammaufrufen GSB E wird bei Label E in Zeile 111 
die aktuelle Zeit t im X-Register und die Spannung u, im Y-Register zur Verfügung gestellt. Mit 
den Anweisungen 112 bis 119 wird der Winkel wt in Grad gebildet. Die Anweisungen 120 bis 126 
dienen zur Erzeugung einer Rechteckspannung mit der in STO 9 abgespeicherten Amplitude. Der 
berechnete Funktionswert u, wird mit der Return-Anweisung in Zeile 139 an das rufende Pro- 
gramm übergeben. 
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4.1.3 Test- und Anwendungsbeispiel 


RC-Integrierglied: R= 10 kQ 
C=1yr 
U=100V 
f=50 Hz =T = 20 ms Periodendauer 


Vor dem Programmstart sind folgende Eingaben erforderlich: 


to =0s in STOA Anfangszeit 
U,=0V in STOB Anfangswert der Spannung u;(0) = uz(- 0) in Volt 
At=0,001s in STOC Schrittweite At 
R = 10?Q in STOD Ohmscher Widerstand 
C=10"°; in STOE Kapazitiver Widerstand 
Re = 100 in STO8 Belastungswiderstand 
U = 100 V in STO9 Amplitude der Eingangsspannung 
T=20:10”°s in STOI Periodendauer 


Start [A]: 


Als Ergebnisse werden die 
Wertpaare t, u» ausgegeben: 


t = 0,001s 

us = 7,977 V 

t = 0,002s 

u> = 16,615 V 

t = 0.0035 

u = 24,354 V 

t = 00108 Bild 4.1.2 Ausgangsspannung eines Integriergliedes 
un = 59,528 V bei rechteckförmiger Eingangsspannung 


Kurve a) Ra = 100 KR 
Kurve b) Rs = 1KkR 


In Bild 4.1.2 sind die Lösungswerte in Form eines Liniendiagramms aufgetragen. Bei der Kurve a 
mit Ra = 10K® ist auch die Verlagerung der Spannung u, Zu einem zeitachsen-symmetrischen 


RR 
Verlauf, bedingt durch das mit der Zeitkonstanten T=C Er == 10 ms abklingende Gleich- 
stromglied, festzustellen. 


4.2 Induktiver Lastkreis mit Freilaufdiode 


Bei der Realisierung von Gleichstromstellern mit Thyristoren oder bei getakteten Netzgeräten 
treten pulsförmige Ansteuerungen von induktiven Lastkreisen auf. Die Induktivität hat dabei die 
Aufgabe, die onmsche Last während der Pulspause aus dem gespeicherten magnetischen Energie- 
potential zu versorgen. Der Stromfluß wird während der Strompause über eine Freilaufdiode ge- 
mäß Bild 4.2.1 geführt. Für einen derartigen Anwendungsfall soll der zeitliche Verlauf des Stromes 
im Lastkreis bestimmt werden. 
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Für den Strom i im Lastkreis gilt die Maschengleichung: L 


d; j 
u=L+ir (4.2.1) I@ m (R 


dB (42.2) Bild 4.2.1 
dt L L Induktiver Lastkreis 
Für die pulsförmige Speisespannung u gilt mit der mit Freilaufdiode 
Periodendauer T: 
u 0s:ıs] 
u-j für. (4.2.3) 
0 5 <t<T 


Die pulsförmige Speisespannung mit der Pulsdauer 3 kann programmtechnisch durch eine Vor- 
zeichenabfrage des Sinus aus dem Argument in Grad 360 f t gewonnen werden: 


u=Ug (4.2.4) 


1 „.. sin360° ft>O 


= 2.5 
so sin3Ko°tr<o (4.2.5) 


Die Gl. (4.2.5) wird mit dem Übergabeparameter t im X-Register durch folgende Anweisungsfolge 
realisiert: 


-ooox-x[|< 
L----L[L[L[[|n 


£ 


LI, sinx 


4.2.1 Programmbeschreibung ‚‚Induktiver Lastkreis mit Freilaufdiode’’ 


Da es sich auch hier wieder um die Lösung einer Differentialgleichung 1. Ordnung handelt (G1.4.2.2), 
braucht lediglich ein Unterprogramm mit der Einsprungadresse Label E ab Programmzeile 111 dem 
Hauptprogramm „Differentialgleichung 1. Ordnung‘ aus Band 1 (Tabelle 4.2.1 bis Zeile 110) hinzu- 
gefügt werden. 

Hierzu ist es zweckmäßig, zunächst die Eingabespeicher zu definieren: 


1. Vom Hauptprogramm vorgegebene Speicher: 


STO A: Anfangszeit t, in Sekunden 
STO B: Anfangswert des Stromes i(to) in Ampere 
STOC: Schrittweite At in Sekunden 
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2. Im Unterprogramm festgelegte Speicher: 


STOD: 
STO E: 
STO 8: 
STO9: 


Spannungsamplitude U in Volt 


Frequenz f der Spannungspulse in Hertz 


Widerstand R in Ohm 
Induktivität L in Henry 


Tabelle 4.2.1 Anweisungsliste „Induktiver Lastkreis mit Freilaufdiode‘ 


Start 


@8l_ wLBli 847 
062  RÜLRA 846 
863  5TC1 849 
064  RCLB 058 
085  STO2 851 
B06 xLBLk 852 
#7 RCL2 853 
Bes Reli 854 
882 CFö 855 
eia 6SBE 855 
all  KÜLE 857 
slz x 856 
aiıs 5705 859 
814 5T06 868 
615 & B6l 
öls = 862 
al? ReLE 663 
818 + 864 
e1l3 RC 855 
26 Z 866 
821 = B657 

22 Reli BEE 
623 + 8659 
a24 $TO> 878 
025  65BE 871 
BZEe  RÜLE ?2 
027 x 673 
BEE 5Tüs B74 
623 2 UP 875 
038 Ä 876 
631  ST+ä 877 
632 4 876 
433 > 673 
834 RÜLE 886 
835 + @81 
036  RÜLS 862 
a3r  ESBE 883 
036 RELC 054 
833 x 885 
ea STOr 886 
a4 Fi? 887 
042 6708 885 
043 65Ba 883 
644 F&? 898 
645 E£TOb [®: 
646 #LBLO 892 


RCL? 895 
2 i 894 

x —895 
$T+8 896 
RCL? @97 
RCL2 896 
+ 094 
RCLI 168 
RCLG 181 
+ 182 
$T03 183 
steil 184 
6SBE 105 
RÜLC 186 
x 107 
$5T+8 166 
[3 189 
$T:8 118 
RCL& UP 111 
ST+2 112 
Rci1 113 
Spt 114 
FSE 115 
FRI =Dt 1l6 
RÄLZ2 117 
FSE 118 
PRIs=Di 113 
6T0k 128 
&LBLo 121 
RÜLF 122 
RCLE 23 
x=r? 124 
£&T02 125 
- 126 
RCL6ö 127 
RCLS 12€ 
- 129 

= 138 
ABS 131 
r 132 

1 133 

5 134 
xY9 135 
gTol 13€ 
P 137 

5 138 
133 


KEY? 
6702 
2 
sF6 
“LBL2 
RELC 
Fo? 
x 
PSE=D> At 
sToc 
RTN 
#LBLE 
xer 
STOr 
Key 
RCLE 
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Das in Tabelle 4.2.1 angegebene Programm beginnt mit der Einsprungadresse Label E in Zeile 111. 
Im X-Register steht die aktuelle Zeit t und im X-Register steht der aktuelle Strom it) zur Ver- 
fügung. In Zeile 113 wird der Strom i{t) nach STO 7 zwischengespeichert. Mit den Anweisungen 
115 bis 120 wird das Argument 360 f t gebildet, so daß nach dem Unterprogrammaufruf Label 3 
die Zeitfunktion g(t) gemäß Gl. (4.2.5) im X-Register zur Verfügung steht. Mit den folgenden 
Anweisungen der Programmzeilen 122 bis 131 wird die Gl. (4.2.2) ausgewertet und mit der Rück- 
sprunganweisung in Zeile 132 das Differential des Stromes di/dt im X-Register an das Hauptpro- 
gramm übergeben. 


4.2.2 Test- und Anwendungsbeispiel 
Der Stromverlauf eines pulsförmig angesteuerten Stromkreises gemäß Bild 4.1.3 soll in Abhängig- 


keit von der Zeit für die ersten fünf Perioden der Spannung u bestimmt werden. Es gelten folgende 
Zahlenwerte: 


U=200 V R=100% 
f = 1000 Hz L=0,1H 
Eingaben: 0 STOA (Anfangszeit ty= 0) 


0 STOB _(Anfangsstrom i(0) =0) 
EEX-4 STOC (Schrittweite At =0,1ms) 
200 STOD (Spannung U = 200 V) 
1000 STOE (Frequenz f = 1000 Hz) 
100 STO8 (Widerstand R = 100.22) 

0,1 STO9 _(Induktivität L=0,1H) 


Programmvorbereitung: Flag 1 setzen, damit die automatische Schrittweitensteuerung unterdrückt 
wird. 


Start [a]: 


Ergebnisse: Es werden die Wertepaare t, i mit der Schrittweite 0,1 ms ausgegeben 


t[s] 
i[A] 


0,0002 0,0003 0,0010 


0,4448 


0,5627 0,6994 


Die Ergebnisse sind in Bild 4.2.2 für die ersten fünf Perioden der Spannung aufgetragen. 


Bild 4.2.2 


Spannungs- und Stromverlauf bei 
pulsförmiger Ansteuerung eines 
induktiven Lastkreises mit Frei- 
laufdiode 

a) Spannung u {t), 

b) Strom it) 


0 1 2 3 4 5 ms 
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3 Programmiertechnik 


3.1 
3.2 
3.3 
3.4 
3.5 
3.6 
3.7 
3.8 
3.9 


Speichertechnik 
Speicherarten 

Indirekte Adressierung 
Speicherregister Arithmetik 
Programmaufrufe 
Dateneingabe 
Ergebnisausgabe 
Programmverzweigungen 
Unterprogrammtechnik 


3.10 Fehlersuche 
3.11 Dokumentation 


4 Angewandte Mathematik 


4.1 Quadratische und kubische Gleichungen 
4.2 Lineare Gleichungssysteme bis vierter Ordnung 
4.3  Newtonsches Iterationsverfahren 
4.4 Numerische Integration 
4.5 Lösung von Differentialgleichungen nach dem Runge-Kutta-Verfahren 
4.6 Komplexe Rechnung 
4.7 Harmonische Analyse 
4.8 Lösung nichtlinearer Gleichungssysteme mit Hilfe der Newton-Raphson-Methode 
5 Finanzmathematik 
5.1 Rentenberechnung 
5.2  Finanzierungsberechnung 
5.3  Zinseszinsberechnung für Jahres- und Monatszyklen 
5.4  Wirtschaftlichkeitsberechnung von Investitionen 
5.5  Einkommensteuerberechnung 
6 Statistik 
6.1 Gaußsche Normalverteilung 
6.2 Binominalverteilung 
6.3 Klassifizierung durch Stichproben 
6.4  Regressionsanalyse 
6.5  D’Hondtsches Verteilungsverfahren 
7 Informatik 
7.1 Konvertierung zwischen Zahlensystemen 
7.2  Code-Umwandlungen 
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Die in diesem Band vorliegende Programmsammlung wurde auf den 


Rechnern HP-67 und HP-97 


123456189061 


ee] 


HEWLETT Ku PACKARD 


erstellt. Aus den verschiedensten Fachbereichen ausgewählt und von durchweg 
hohem Anforderungsniveau, dokumentieren diese Programme beispielhaft die viel- 
seitigen Einsatzmöglichkeiten und den hohen technischen Standart der Rechner 
HP-67 und HP-97. Zusammen mit der bewährten Computer-Logik UPN, der ein- 
fachen Programmierung und dem anspruchsvollen Qualitätsmaßstab, den HEWLETT- 
PACKARD an alle seine Produkte anlegt, stehen so dem Anwender Rechner zur 
Verfügung, die im besten Sinne benutzerfreundlich sind. 

Aber HEWLETT-PACKARD bietet noch mehr. Ergänzend zu den Rechnern liegt 
ein umfassendes Software-Angebot von mehr als 3000 weiteren Programmen aus 
Wirtschaft, Wissenschaft und Technik in praxisgerechter Form vor. Damit eröffnet 
sich den Benutzern von HEWLETT-PACKARD-Rechnern die Möglichkeit, die 
individuellesten Problemstellungen gezielt anzugehen und einer schnellen und 
exakten Lösung zuzuführen. 


ES 


Software-Information für den Benutzer von HP-Rechnern 


Baustatik: 
Umfangreiche Programmsammlung in 


DIN-Norm von Ing. (grad.}) W. Mücke: 


Berechnung von Einfeld- und Durch- 
laufträger, Beton-Stahlbeton-Bemes- 
sung nach DIN 1045, Fundamente, 
Rahmen, Ing.-Holzbau, Stahl- und 
Grundbau, Spannbeton- und Brücken 
bau. 


Bauphysik: 

Zahlreiche Programme für Heizungs-/ 
Klimatechnik, Akustik und Sanitär 
von Ing. (grad.) H. Markert 


Maschinenbau: 


Programme zur Berechnung von 
zylindrischen Schraubenfedern, 
Schrumpfen, Verzahnungen, Rohr- 
leitungen, mehrfach abgesetzte 
Maschinenwellen von Prof. Dr.-Ing. 
M. Wink 


Geodäsie: 


Programme aus dem Bereich Straßen- 
planung und Absteckung, Ingenieur- 

Geodäsie, Landesvermessung, Photo- 
grammetrie von Prof. Dipl.-Ing. 

H. Oberste-Brink-Bockholt 


Navigation: 

Programme zur Berechnung der 
astronomischen Standordbestimmung, 
Gestirnsuche, Koppelrechnung, 
Orthodrome, Luxodrome von 

R. Wurster 


Finanzmathematik: 

Berechnung der Einkommens-, Ge- 
werbe-, Körperschafts- und Lohnsteuer, 
Zeitwertberechnung, Betriebsrentenan- 
passung, Kontokorrenzrechnung von 
Dr. A. Braun 

Braess/Fangmeyer „Elektronische 
Effektivzinsberechnung‘, Annuitäten- 
darlehen, Darlehen und Wertpapiere 
mit Patenttilgung, Festdarlehen und 
gesamtfällige Wertpapiere von Dipl. 
Mathematiker H. Fangmeyer und 
Bankkaufmann A. Beer 


Forst- u. Waldwirtschaft: 


Berechnung des Aufhöhenfaktors, 
Aufhöhe, Mittendurchmesser, gesamte 
Erntefestmeter ohne Rinde. Inhalt des 
Einzelstammes von R. Düser 


Sämtliche Programme mit ausführlicher 
Dokumentation incl. Magnetkarten. 


Umfangreiche Software aus unserem Hause: 
Statistik — Maschinenbau — Finanzmathematik — Elektronik — Spiele — 
Vermessung — Navigation — Nuklearmedizin. 


„Users-Club” für Programmaustausch: 


Der HP Users-C/ub zählt über 5000 Mitglieder. Es gibt über 3000 Programme. 
Zusätzlich stehen weitere Programmsammlungen mit jeweils 10 bis 15 Program- 
men, ohne Magnetkarten, unterschiedlicher Sachgebiete zur Verfügung. 


Adressenliste: Weitere Programmsammlungen von Anbietern der verschiedensten 


Fachrichtungen auf Anfrage. 


Fordern Sie weitere Informationen beim Fachhandel an oder direkt bei: 


Hewlett-Packard GmbH/Vertrieb, Berner Straße 117, 
6000 Frankfurt 56, Tel. (0611) 5004-1 


Die unentbehrliche Orientierungs- 
und Anwendungshilfe 


Etwa zehn Jahre sind vergangen, seit die ersten Mikroprozessoren 
und Taschenrechner auf den Markt kamen. Seit etwa acht Jahren 
werden diese „Kleinstcomputer’' mit den vorher kaum für möglich 
gehaltenen Fähigkeiten in großen Stückzahlen benutzt. Und nur 
rund fünf Jahre ist es her, seit programmierbare Taschenrechner und 
„Jedermann-Computer‘ (Personal Computer) zur Verfügung stehen. 
Das „Taschenrechner + Mikrocomputer Jahrbuch’ wendet sich an 
den großen Leserkreis, der sich über diese rasante technische Ent- 
wicklung und deren Anwendung in den verschiedensten Lebens- 
bereichen systematisch und gezielt informieren will. 

Das Jahrbuch-Konzept geht von einer Zweiteilung aus: Der Fachteil 
enthält aktuelle Beiträge zu den Bereichen Taschenrechner, Mikro- 
computer und Speicher. Der Datenteil besteht im wesentlichen aus 
Produktübersichten, Bezugsquellen, Adressen, Literatur- und Zeit- 
schriftenlisten. Ein detailliertes Sachwortverzeichnis erleichtert das 
Auffinden besonderer Textstellen. 

Um auch weniger erfahrenen Lesern entgegenzukommen, haben die 
Fachbeiträge ein Glossar. Darin sind in knapper Form Fachaus- 
drücke erklärt, die im jeweiligen Aufsatz ohne Kommentar benutzt 
werden, 

Die Produktübersichten (Datentabellen) sind so vollständig wie 
irgend möglich gehalten. Das bedeutet, es sind auch Taschenrechner, 
Tischrechner, Mikroprozessoren sowie Mikrocomputer aufgenom- 
men, die nicht mehr hergestellt und verkauft werden. Dadurch ent- 
steht aber gerade der beabsichtigte systematischen Überblick. So 
werden beispielsweise Besitzer älterer Geräte in die Lage versetzt, 
eine Entscheidung darüber zu treffen, ob eine Neuanschaffung lohnt 
oder gar notwendig ist. 

Das gewählte Konzept stellt sicher, daß ältere Jahrbücher nicht von 
vornherein zum Bücherschrank- oder Papierkorbfüller abgewertet 
sind. Vielmehr werden sie als oft genutzte Nachschlagewerke und als 
Chronik der „Mikrorechnerentwicklung’ über Jahre hin nützlich 
sein. 


Taschenrechner + 
Mikrocomputer 
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Taschenrechner + 
Mikrocomputer 
Jahrbuch 


Anwendungsbereiche — Produkt- 
übersichten — Programmierung — 
Entwicklungstendenzen — 
Tabellen — Adressen 


Herausgegeben von 

Harald Schumny 

Mit zahlr. Abb., Tabellen und 
Adressen. Ca. V1, 254 S. 

18,5 X 24 cm. Kart. 


Info-Gutschein 


Bitte informieren Sie mich (uns) ständig über Ihre 
Neuerscheinungen auf dem Gebiet: 


UL Taschenrechner U) Mikrocomputer 
Ich {wir) besitze (n) folgendes Gerät: 
TR: 


uC: 
Hauptanwendungsgebiete des TR bzw. uC: 


Diese Karte entnahm(en) ich (wir) dem Buch: 
Alt, Anwendung programmierbarer Taschenrechner, Band 1 


Meine (unsere) Buchhandlung: 


Gleichzeitig bestelle (n) ich (wir) folgende Bücher: 


Anzahl | Autor und Titel 


| Schumny, TR + a„C-Jahrbuch 
| 
I 
t 
| 
l 


Anschrift: 


Beruf /Branche 


Datum Unterschrift 


Friedr. Vieweg & Sohn Verlagsgesellschaft mbH : Wiesbaden 


Lieber Leser! 


Nenn Sie Interesse haben, aktiv an der 
Neiterentwicklung unseres Literaturprogramms 
zum Bereich TR + uC mitzuarbeiten, z.B. durch 
Veröffentlichung ausgetesteter Programme zu 
Jestimmten Anwendungsgebieten, dann schreiben 


Sie uns. 


Wir freuen uns über Ihre Nachricht und werden 
uns umgehend mit Ihnen in Verbindung setzen. 


Mit freundlichem Gruß 
Lektorat Fachbuch 


gez. Niclas 


Antwort 


Friedr. Vieweg & Sohn 
Verlagsgesellschaft mbH 


Postfach 5829 


D-6200 Wiesbaden 1 


Bitte 
mit 


50 Pf. 
freimachen! 


Anwendung programmierbarer Taschenrechner 


Diese Reihe bietet den Benutzern programmierbarer Taschenrechner eine reich- 
haltige Palette von Aufgabenstellungen aus den Anwendungsgebieten der Natur- 
und Wirtschaftswissenschaften an, für die Programme zur numerischen Lösung 
entwickelt werden. 


Jeder Band behandelt ein in sich abgeschlossenes Themengebiet: Nach einer kurzen 
Einführung in die Theorie der jeweiligen Problemstellung wird der Lösungsalgorith- 
mus entwickelt, das Programm dargestellt und kommentiert. 


Neben der direkten Nutzung der hier veröffentlichten Programme unterstützt diese 
Reihe den Leser wirkungsvoll bei der Ausarbeitung eigener Programmvarianten. 


Band 2: Allgemeine Elektrotechnik, Nachrichtentechnik, Impulstechnik 
für UPN-Rechner 


von Helmut Alt 


Band 2 enthält zahlreiche Programme für die Rechnertypen HP-67/97 und HP-41C 
zu folgenden Problemstellungen: Netzwerkberechnung /Elektrische Übergangsvor- 
gänge / Freie Schwingungen in elektrischen Stromkreisen / Netzwerkumwandlung / 
Harmonische Analyse /Schwingkreisberechnung /Wellen- und Betriebsparameter 
passiver Filter /Aktive Filter / Übertragungsverhalten elektrischer Netzwerke /Lei- 
tungsgleichungen /Zeitverhalten von Übertragungsgliedern /Induktiver Lastkreis 
mit Freilaufdiode. 


Dr.-Ing. Helmut Alt ist Mitarbeiter der Rheinisch-Westfälisches Elektrizitätswerk AG 
(RWE), Betriebsverwaltung Düren. 


ISBN 3-528-04151-X 


